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I.  INTRODUCTION 


This  report  presents  the  sixth  annual  summary  of  research  at  Ohio 
State  sponsored  by  the  Joint  Services  Electronics  Program  (JSEP).  The 
research  is  In  the  area  of  electromagnetics  and  the  specific  topics  are: 
(1)  Diffraction  Studies;  (2)  Hybrid  Techniques;  (3)  Integral  Equation 
Studies;  (4)  Time  Domain;  and  (5)  Transient  Signature  Measurements  of 
Radar  Targets  for  Inverse  Scattering  Research. 

The  following  sections  summarize  the  significant  accomplishments  of 
the  program  (Section  II)  and  the  research  by  work  unit  (Section  III). 
Researchers  and  their  publications  are  listed  under  each  work  unit.  A 
listing  of  the  present  research  programs  at  the  Laboratory  and  all 
reports  and  papers  published  by  the  Laboratory  during  the  past  year  are 
given  in  the  appendices. 
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II.  SIGNIFICANT  ACCOMPLISHMENTS 


Our  major  research  area  continues  to  be  the  analysis  of 
electromagnetic  radiation  and  scattering.  Associated  with  this  primary 
goal  Is  a  substantial  program  for  Improving  our  experimental 
capabilities.  This  was  supported  in  part  by  the  JSEP  program  and  has 
expanded  recently  with  other  sponsorship  so  that  our  anechoic  chamber 
compact  scattering  range  is  probably  the  best  of  Its  kind  at  the  present 
time.  This  Is  emphasized  by  the  fact  that  we  have  provided  information 
to  Scientific  Atlanta  (the  reflector  manufacturer)  for  Improving  these 
facilities.  We  have  also  been  contacted  by  several  other  manufacturers 
requesting  assistance  In  building  the  next  generation  of  compact  ranges. 
It  Is  noted  that  the  stepped  frequency  concept  originated  under  JSEP  was 
then  developed  using  the  compact  range  under  separate  funding  by  ONR. 
Additional  Improvements  were  made  under  funding  from  NASA.  This  has 
Indeed  been  a  most  rewarding  area  of  research. 

The  analysis  is  being  directed  more  and  more  to  the  study  of  the 
Interaction  of  electromagnetic  fields  and  penetrable  materials.  This  is 
being  pursued  from  (a)  a  high  frequency  asymptotic  viewpoint,  (b)  a 
surface  Impedance  model,  (c)  a  moment  method  approach  and  (d)  a  hybrid 
approach. 

However,  the  studies  of  diffraction  by  conducting  bodies  using  the 
Uniform  Geometrical  Theory  of  Diffraction  (UTD)  continue  to  be  a  very 
substantial  part  of  our  research  activities.  The  general  problems 
consist  of  the  radiation  from  antennas  mounted  close  to  perfectly 


conducting  edges  (and  also  the  scattering  from  objects  positioned  at 
edges),  as  well  as  paraxial  edge  diffraction  and  slope  diffraction  by 
perfectly-conducting  convex  cylinders.  The  radiation  studies  have 
helped  us  to  develop  more  general  computer  codes  on  other  programs  for 
the  analysis  of  antennas  on  or  near  complex  structures  (e.g.,  aircraft 
and  ships).  The  scattering  work  is  helping  on  programs  involving  radar 
cross  section  studies  and  target  Identification.  This  work  Is  basic  to 
the  development  of  computer  codes  for  calculating  the  patterns  of 
reflector  antennas  and  antennas  on  aircraft,  missiles,  satellites,  ships 
and  in  other  environments,  and  for  calculating  the  scattered  fields  of  a 
wide  range  of  complex  objects. 

As  we  have  observed,  some  of  the  above  efforts  have  been  supported 
by  other  sources  of  funding  but  are  a  direct  offspring  of  the  Joint 
Services  Electronics  Program.  Such  topics  Include  but  are  not  limited 
to: 

a)  Improved  RCS  pattern  measurement  techniques, 

b)  analysis  of  the  scattered  fields  of  aerospace  vehicles,  and 

c)  use  of  transient  techniques  for  low  frequency  target 
center  Isolation. 

Computer  codes  based  on  the  high  frequency  asymptotic  techniques 
for  antennas  mounted  on  various  types  of  vehicles  are  now  used  through 
out  the  aerospace  Industry.  A  long  term  goal  In  the  electromagnetic 
scattering  studies  Is  to  evolve  computer  codes  comparable  to  those 


for  antenna  patterns  now  In  general  use  by  the  aerospace  Industry.  One 
special  code  has  been  developed  for  this  purpose  but  requires 
substantial  knowledge  of  the  6TD  and  of  the  pertinent  scattering 
mechanism  by  the  user. 

The  time  domain  studies  In  electromagnetics  are  also  continued. 
Research  In  both  theoretical  and  experimental  methods  are  providing 
substantial  rewards.  The  theoretical  program  Involves: 

1.  Complex  Natural  Resonance  and  Geometrical  Procedures 

2.  Evaluation  of  Complex  Natural  Resonances  via 
Rational  Function  Approxlmants 

3.  Transient  Current  Density  on  a  Conducting  Sphere 

4.  Cavity  Studies,  and 

5.  K-Pulse  Studies. 

The  time  domain  experimental  program  Is  directed  toward  assembling 
complete  complex  scattering  data  on  several  canonical  target  shapes. 

The  complete  target  signature  data  Increases  substantially  In  value  when 
the  waveform  can  be  Interpreted  and  scattering  mechanisms  can  be 
Identified.  These  experimental  results  are  of  sufficient  quality  that 
this  has  become  practical.  Indeed,  they  have  been  used  several  times 
this  past  year  on  related  projects  when  unexpected  high  scattered  fields 
were  obtained  for  aircraft  targets  to  Identify  and  reduce  these  high 
scattering  mechanisms. 


Continued  research  In  time  domain  reflection  will  produce  a  more 
complete  understanding  of  the  Interaction  of  electromagnetic  fields 
with  complex  targets  and  will  also  Improve  our  target  Identification 
capabilities. 

Nine  graduate  students  have  been  Involved  in  the  above  research 
over  the  past  year.  Over  the  past  four  years,  with  the  support  of  the 
JSEP  at  Ohio  State,  there  have  been  10  students  granted  the  M.Sc.  degree 
In  Electrical  Engineering  and  9  students  granted  the  Ph.D.  degree  In 
Electrical  Engineering. 


III.  RESEARCH  SUMMARY 


A.  Diffraction  Studies 

Researchers:  R.G.  Kouyoumjlan,  Professor  (Phone:  (614)  422-7302) 
P.H.  Pathak,  Assistant  Professor 
N.  Wang,  Senior  Research  Associate 
R.  Tlberlo,  Visiting  Professor  and  Consultant 
R.  Paknys,  Graduate  Research  Associate 
M.  Buyukdura,  Graduate  Research  Associate 


Accompl 1 shments 


During  the  present  contract  period,  the  work  accomplished  In 
extending  the  uniform  geometrical  theory  of  diffraction  (UTD)  has  been 
substantial.  This  research,  and  the  technical  papers  based  on  this 
research  which  have  recently  appeared  (or  have  been  accepted  for 
publication)  are  described  below. 


1.  Diffraction  by  Non-Conducting  Surfaces 

a.  Smooth  dielectric  covered  and  Impedance  convex  surface 

A  study  of  the  electromagnetic  radiation  and  scattering  from  a 
conducting  surface  with  dielectric  coating  Is  of  great  Interest  In  that 
It  provides  an  understandl ng  of  the  effects  of  the  loading  on  the 
scattered  fields.  An  Interesting  application  of  coating  Is  to  control 
the  electromagnetic  scattering  characteristics  from  conducting  bodies 
such  as  an  aircraft,  missile,  satellite,  etc.  Also,  It  Is  useful  In 
predicting  the  radar  cross-section  of  structures  made  of  composite 
materials  or  conducting  bodies  coated  with  dielectric  materials. 

In  the  case  of  a  conducting  surface  coated  with  a  thin  layer  of 
dielectric  with  a  uniform  thickness,  the  surface  can  be  conveniently 
viewed  as  an  Impedance  surface  with  a  constant  surface  Impedance.  We 
have  developed  a  high  frequency  solution  for  the  problem  and  were  able 
to  predict  the  resonance  phenomena  of  the  radar  cross-section  of  a 
circular  cylinder  with  a  constant  surface  Impedance.  It  was  found  that 
surface  waves  with  almost  pure  Imaginary  propagation  constant  traverse 
around  the  cylinder  surface  with  diminishing  attenuation,  and  Interfere 
with  each  other  constructively  such  that  they  add  In  phase  to  give  the 
distinctive  resonance  phenomena  In  the  radar  cross-section.  Numerical 
values  for  the  propagation  constant  of  the  surface  wave,  which  are 
related  to  the  Regge  poles  of  the  Impedance  cylinder,  have  been  found. 
Also,  a  criterion  for  predicting  resonance  has  been  established  and 


the  correlations  between  the  resonance,  the  Regge  poles,  and  the  natural 
frequencies  of  the  Impedance  cylinder  have  been  demonstrated. 

It  should  be  noted  that  the  surface-impedance  model  Is  valid  only 
for  low-loss,  thin  dielectric  sheets  over  the  conducting  surfaces.  For 
a  coating  with  moderate  thickness,  a  more  accurate  approach  Is  being 
pursued.  More  recently,  work  has  been  carried  out  on  the  Investigation 
of  the  radar  cross-section  of  a  conducting  cylinder  coated  by  a 
dielectric  layer  of  uniform  thickness. 

Using  the  standard  Watson's  transformation  technique,  the  rigorous 
eigenfunction  solution  for  the  coated  circular  cylinder  Is  cast  Into  a 
ray  solution.  The  backscattered  field  from  the  coated  cylinder. 
Illuminated  with  an  Incident  plane  wave.  Is  obtained  by  summing  the 
geometrical -optics  contribution  and  the  surface-wave  contributions  which 
Include  all  the  multiply-encircled  surface  waves.  The  parameter  of 
fundamental  Importance  associated  with  surface-wave  fields  Is  the 
propagation  constant  of  the  surface  wave.  The  propagation  constant  Is 
related  to  the  Regge  poles  of  the  coated  cylinder.  By  a  combination  of 
numerical  procedures,  these  Regge  poles  were  obtained  by  solving  the 
high  frequency  approximation  of  the  resonance  equation  for  the 
dielectric  coated  conducting  cylinder. 

For  thick  enough  coatings,  the  surface  waves  moving  around  the 
cylinder  resemble  those  trapped  In  a  planar,  grounded,  dielectric  slab. 
It  Is  found  that  the  dominant  surface  waves  have  a  small  attenuation 
rate  so  that  all  the  multiply-encircled  waves  must  be  taken  Into 


consideration  In  a  sel f-consl stent  fashion.  As  a  consequence, 
resonances  appear  In  the  scattering  cross  section.  The  surface  wave 
Interpretation  of  the  resonance  phenomenon  and  the  condition  for 
predicting  the  resonance  frequencies  are  now  understood.  The  normalized 
radar  cross  section  of  the  coated  cylinder  obtained  from  the  high 
frequency  ray  solution  has  been  shown  to  compare  well  with  the 
eigenfunction  result. 

As  a  continuation  and  extension  of  the  work  Involving  the 
backscatterlng  from  the  coated  cylinder,  efforts  have  been  directed  to 
the  problems  of  radiation  and  bl static  scattering  from  smooth, 
dielectric  covered  and/or  Impedance  convex  surfaces.  The  ultimate  goal 
here  Is  to  develop  a  uniform  high-frequency  solution  for  the 
electromagnetic  fields  due  to  an  antenna  radiating  In  the  presence  of  a 
loaded  convex  surface.  This  solution  could  be  expressed  In  terms  of  an 
Integral  representation  similar  to  the  Fock-type  Integrals  which 
appeared  In  the  solution  for  the  perfectly-conducting,  smooth,  convex 
surfaces  [1,2,3]  except  that  the  Fock  Integral  now  contains  a 
denominator  In  the  Fock  Integrand  which  yields  a  resonant  condition  In 
terms  of  the  Impedance  or  dielectric  loading.  Far  from  the  geometrical 
shadow  boundaries,  the  Integrals  can  be  well  approximated  by  the 
geometrical  optics  (GO)  solution  In  the  Illuminated  region  and  by  the 
surface  (creeping)  or  ray  modal  residue  series  solution  in  the  shadow 
region,  respectively.  However,  within  the  transition  region,  l.e..  In 
the  neighborhood  of  the  shadow  boundaries,  only  the  Integral 
representation  provides  the  correct  solution.  Furthermore,  this 


transition  Integral  joins  smoothly  with  the  GO  solution  In  the  lit 
region  and  reduces  to  the  surface  ray  modal  or  residue  series  solution 
In  the  shadow  region.  Thus,  the  essential  task  Is  to  develop  an 
efficient  method  for  evaluating  this  Fock-type  transition  Integral 
associated  with  the  smooth,  dielectric-covered  and  Impedance  type  convex 
surface. 

The  Fock  (transition)  Inteqral  g(x,q),  associated  with  the 
canonical  problem  of  the  radiation  by  a  two  dimensional  (2-D)  magnetic 
line  source  on  a  circular  cylinder  with  a  constant  surface  Impedance 
boundary  condition,  has  been  derived  and  evaluated.  Some  numerical 
results  for  g(x,q)  are  presented  In  Figure  A-l.  Referring  to  Figure 
A-l,  our  results  for  a  lossy  curved  surface  (A  *  0.5,  1.0)  compare  well 
with  those  obtained  by  Walt  and  Conda  [4].  Note  that  g(x,q)  represents 
the  field  radiated  by  the  line  source  on  the  cylinder,  or  It  also 
represents  the  field  Induced  on  the  surface  of  the  cylinder  by  a  distant 
line  source  (located  off  the  surface)  via  reciprocity.  For  a  reactive 
surface  that  supports  a  surface  wave,  our  results  (q>0)  confirm  the  fact 
that  a  reactive  surface  Is  more  efficient  than  the  conducting  surface 
for  guiding  the  electromagnetic  energy.  Notice  that  In  Figure  A-l,  the 
curve  denoted  by  q*0  Is  for  the  case  of  a  perfectly  conducting  curved 
surface.  Here,  q  Is  a  parameter  related  to  the  value  of  the  surface 
Impedance  and  the  cylinder  radius.  The  shadow  boundary  Is  located  at 
x=0  In  Figure  A-l.  The  above  material  taken  from  [5]  has  now  been 
extended  In  the  current  time  frame  as  Is  discussed  In  the  following 
paragraphs. 


The  newly  developed  transition  Integral  g(x,z)  has  been  employed  to 


calculate  the  high  frequency  radiation  patterns  due  to  a  two 
dimensional  magnetic  line  source  mounted  on  a  circular  cylinder  with 
surface  Impedance  Zs.  Figures  A-2  through  A-5  compare  the  patterns 
obtained  by  the  eigenfunction  solution  (solid  line)  and  the  high 
frequency  GTD  solution  (dash  line).  It  should  be  noted  that  the  high 
frequency  GTD  solution  contains  the  newly  developed  transition  Integral 
g(x,q),  where 


q  =  -1 


1/3 


Zs 

*0 


Notice  the  agreement  between  the  two  solutions  Is  satisfactory  for 

Zs  *  Rs  +  j  Xs  up  to  10  +  j  200.  However,  for  Zs  =  10  +  j  300,  the  high 

frequency  GTD  solution  does  not  yield  accurate  results.  The  dlscrepency 

(2) 

Is  Introduced  by  the  approximation  used  for  the  Hankel  function  Hy  (x), 
namely 

( o\  1  pl/3 

and 
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where  Wj(t)  and  Wj ( t)  are  the  Airy  Integrals  defined  by  [4] 
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EIGENFUNCTION  SOLUTION 


- GTD  SOLUTION 

Figure  A-2.  Pattern  due  to  a  magnetic  line  source  radiating  from  a 
circular  cylinder  with  surface  impedance  Zs=Rs+jX$. 
ka  =  10,  Zs  =  lOn. 
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Figure  A-3. 


- EIGENFUNCTION  SOLUTION 
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Pattern  due  to  a  magnetic  line  source  radiating  from  a 
circular  cylinder  with  surface  Impedance  Zs*Rs+jXs. 
ka  *  10,  Zs  =  10  +  j  100. 


Figure  A-4.  Pattern  due  to  a  magnetic  line  source  radiating  from  a 


in  the  notation  of  Fock  where  the  quantity  t  is  defined  by 

2  1/3 

t  =  (“)  (v-x)  . 

It  turns  out  that  for  large  reactance  Xs,  the  curved  surface  supports 

surface  waves  which  are  tightly  trapped.  As  a  consequence,  a  better 

(21 

approximation  for  the  function  fr  (x). 


should  be  employed  [5]  in  the  analysis. 

During  the  previous  year,  the  high  frequency  GTD  solution  has  been 
modified  to  include  the  higher-order  terms  resulted  from  the 
approximations  employed  for  the  Hankel  functions.  Figures  A-6  through 
A-9  compare  the  new  GTD  solution  and  the  eigenfunction  solution.  It  is 
seen  that  the  agreement  between  the  two  solutions  is  excellent. 


EIGENFUNCTION  SOLUTION 
GTO  SOLUTION 


Figure  A-6.  Pattern  due  to  a  magnetic  line  source  radiating  from  a 
circular  cylinder  with  surface  Impedance  Zs=Rs+jXs. 
ka  *  10,  Zs  =»  10. 
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-  EIGENFUNCTION  SOLUTION 

- GTO  SOLUTION 


Pattern  due  to  a  magnetic  line  source  radiating  from 
circular  cylinder  with  surface  Impedance  Zs=R$+v1Xs. 
ka  *  10,  Zs  =  10  +  j  100. 


-  EIGENFUNCTION  SOLUTION 

- GTD  SOLUTION 


Figure  A-8.  Pattern  due  to  a  magnetic  line  source  radiating  from  a 
circular  cylinder  with  surface  Impedance  Zs*Rs+jXs. 
lea  =  10,  Zs  =  10  +  j  200. 


Another  topic  of  Interest  Is  the  mutual  coupling  between  antennas 
mounted  on  a  conducting  surface  covered  with  a  layer  of  dielectric 
material.  Work  has  been  In  progress  to  study  the  surface  Impedance 
model  of  the  dielectric-coated  surface.  It  Is  known  that  the  surface 
Impedance  model  provides  a  convenient  tool  to  study  the  electromagnetic 
characteristics  of  the  surface  waves  associated  with  the 
dielectric-coated  surface.  Some  preliminary  results  for  the  surface 
fields  due  to  a  magnetic  line  source  located  on  a  circular  cylinder  with 
reactive  surface  Impedance  have  been  obtained.  Figure  A-lOa  presents 
the  surface  fields  as  a  function  of  angular  separation  from  the  source, 
using  the  normalized  surface  Impedance,  ns  =  Zs/Z0,  as  a  parameter.  It 
can  be  seen  from  Figure  A-lOa  that  an  increase  In  the  value  of  the 
Inductive  impedance  gives  rise  to  surface  waves  which  are  more  tightly 
trapped  on  the  surface.  As  a  consequence,  the  field  Intensity  Increases 
In  the  shadow  regions  with  ?n  Increase  In  the  Inductive  reactance  of  the 
surface  and  also  the  standing-wave  feature  Is  more  profound.  Work  also 
has  been  carried  out  to  Isolate  the  dominant  contribution  to  the  surface 
field  on  the  cylinder.  Figure  A-lOb  shows  the  "one-ray"  solution  of  the 
problem  given  In  Figure  A-lOa .  Work  Is  now  in  progress  to  study  the 
high  frequency  characteristics  of  the  one-ray  solution.  It  should  be 
noted  that  this  one-ray  solution  Is  related  to  the  universal  transition 
1 ntegral 

W2(t)e-jXt 
Jt  W^(t)-qW2(t)  dt 
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Figure  A-lOa.  Surface  field  produced  by  a  magnetic  line  source  M 
on  the  circular  cylinder  with  a  surface  impedance  ns 
(eigenfunction  solution). 


Figure  A-lOb.  Dominant  (one-ray)  contribution  to  the  surface  field 
produced  by  a  magnetic  line  source  M  on  a  circular 
cylinder  with  a  surface  impedance  nc* 


b.  Discontinuity  of  surface  Immpedance 

The  problems  of  EM  scattering  by  discontinuities  In  surface 
Impedance  arise  not  only  In  the  control  of  the  RCS  of  perfectly 
conducting  structures  by  Impedance  loading,  but  they  also  arise,  for 
example.  In  the  design  of  flush  mounted  alrcraft/mlsslle  antennas  with  a 
thin  coating  (or  a  radome),  of  finite  extent,  and  In  the  studies  for 
Improving  the  Isolation  among  a  pair  of  flush  mounted  antennas  by 
placing  an  absorber  (which  can  be  modeled  by  a  lossy  Impedance  surface) 
In  between  those  two  antennas. 

A  problem  of  Interest  In  the  area  of  scattering  Is  that  of  plane 
wave  diffraction  by  a  strip  with  two  face  Impedance  (l.e.,  the  Impedance 
on  the  two  sides  of  the  strip  can.  In  general,  be  different).  This 
problem  has  been  analyzed  for  the  special  case  of  grazing  incidence;  the 
analysis  for  Illumination  at  aspects  away  from  grazing  on  the  strip  Is 
actually  simpler  than  the  case  of  grazing  Incidence.  It  Is  noted  that 
the  analysis  In  the  forward  scatter  direction  for  grazing  Incidence  Is 
complicated  because  the  diffraction  from  the  leading  edge  produces  a 
non-ray  optical  field  at  the  trailing  edge;  consequently,  the 
diffraction  of  this  non-ray  optical  field  from  the  trailing  edge  which 
contributes  to  the  field  diffracted  in  the  forward  direction  must  be 
handled  carefully.  Here,  the  diffraction  of  the  non-ray  optical  field 
Is  analyzed  via  a  spectral  extension  of  the  geometrical  theory  of 
diffraction.  A  paper  has  been  written  recently  which  describes  this 
work;  namely: 


"Scattering  by  a  Strip  with  Two  Face  Impedances  at  Edge-On 
Incidence",  by  R.  Tlberlo,  F.  Bessl,  6.  Manara,  and  G.  Pelosi, 
accepted  for  publication  In  J.  Radio  Science. 

The  geometrical  theory  of  diffraction  solution  employed  In  the  above 
analysis  (which  Is  formulated  In  the  spectral  domain)  Is  based  on  the 
work  of  Maliuzhlnets  [6].  The  result  for  the  diffraction  coefficient 
available  from  Maliuzhlnets  work,  dealing  with  the  problem  of  the 
diffraction  of  a  plane  wave  by  a  wedge  with  two  face  impedances.  Is 
non-uniform.  Consequently,  Maliuzhlnets  results  had  to  be  modified  here 
so  that  It  would  remain  uniformly  valid  across  the  geometrical  optics 
Incident  and  reflection  shadow  boundaries  formed  by  the  edge  of  the 
wedge . 

A  second  configuration  of  Interest  which  has  also  been  Investigated 
recently,  and  which  Is  somewhat  related  to  the  previous  one,  deals  with 
the  problem  of  diffraction  of  an  electromagnetic  plane  wave  by  the  edges 
of  an  Impedance  surface  patch  on  a  planar  perfectly-conducting  boundary 
of  Infinite  extent.  In  this  Investigation,  an  asymptotic  high  frequency 
solution  for  the  surface  field  Is  obtained  again  using  the  spectral 
theory  In  conjunction  with  a  uniform  version  of  Maliuzhlnets  result  [6]. 
Both  the  TE  and  TM  polarizations  are  considered  in  this  two-dimensional 
problem,  and  the  contribution  from  the  bound  surface  waves  (on  the 
Impedance  surface)  Is  Included  in  the  solution.  Furthermore,  this 
analysis  is  further  extended  to  treat  the  diffraction  by  a  planar 
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three-part  Impedance  surface  In  which  diffracted  rays  are  produced  by 
the  Illumination  at  the  edges  of  the  Impedance  surface  patch  that  lies 
In  the  middle  of  that  planar  three-part  impedance  surface.  A  paper 
describing  this  work  has  been  written  recently,  namely: 

"High  Frequency  Scattering  From  the  Edges  of  Impedance 
Discontinuities  on  a  Flat  Plane",  by  R.  Tlbereo  and  6.  Pelosl, 

IEEE  Transactions  on  Antennas  and  Propagation,  Vol .  AP-31, 

No.  4,  pp.  5910-596,  July  1983. 

In  addition  to  the  above  problems,  work  has  been  completed  on  the 
analysis  of  the  scattering  of  an  EM  plane  and  surface  waves  by  a  planar, 
two  part  Impedance  surface  In  which  one  part  Is  a  perfect  electric  or 
magnetic  conductor  and  the  other  part  is  characterized  by  a  non-zero 
surface  Impedance.  This  two-dimensional  analysis  Is  based  on  the 
Wlener-Hopf  procedure  with  Weinstein's  method  of  factorization  [7]. 

These  Wiener-Hopf  based  solutions,  which  have  been  obtained  for  both  the 
TE  and  TM  polarizations,  provide  uniform  (GTD)  diffraction  coefficients 
for  the  discontinuity  In  surface  Impedance,  and  besides  being  solutions 
to  Important  canonical  problems,  they  are  also  crucial  to  the 
development  of  uniform  GTD  solutions  for  the  Important  problems  of  the 
diffraction  by  a  thin  dielectric  or  ferrite  half  plane,  and  by  a  thin 
dielectric  or  ferrite  half  plane  on  a  perfectly  conducting  surface.  The 
problem  of  the  diffraction  by  a  thin  dielectric  or  ferrite  half  plane  Is 
discussed  later  In  part  c. 


An  analysis  of  the  high  frequency  EM  radiation  by  a  magnetic  line 
(or  line  dipole)  source  on  a  uniform  impedance  surface  which  partially 
covers  a  smooth  perfectly-conducting  convex  surface  has  been  performed 
in  the  past.  A  manuscript  describing  this  work, 

"An  Approximate  Asymptotic  Analysis  of  the  Radiation  from 
Sources  on  Perfectly-Conducting  Convex  Cylinders  with  an 
Impedance  Surface  patch",  by  L.  Ersoy  and  P.H.  Pathak, 

has  nearly  been  completed.  The  preceding  work  is  of  Interest,  for 
example,  in  the  study  of  fuselage  mounted  airborne  slot  antennas  for  the 
purposes  of  controlling  the  field  radiated  near  the  horizon  (or  shadow 
boundary)  over  that  In  the  absence  of  the  patch. 

The  work  on  the  diffraction  problems  described  above,  dealing  with 
structures  possessing  a  discontinuity  in  surface  impedance,  will  be 
continued  in  the  future  phases  of  this  study  in  order  to  extend  these 
solutions  to  treat  more  general  situations  (e.g.,  to  treat  3-D 
problems). 


c.  Diffraction  by  the  edge  of  a  thin  dielectric  or  ferrite 
half  plane 


The  diffraction  by  a  thin  dielectric  or  ferrite  half  plane  is  an 
important  canonical  problem  in  the  study  of  the  diffraction  of 
electromagnetic  waves  by  penetrable  bodies  with  edges.  The  excitation 
chosen  for  studying  this  problem  is  either  an  electromagnetic  plane 
wave,  or  a  surface  wave  incident  along  the  dielectric  surface;  both 
types  of  excitation  are  considered.  For  sufficiently  thin  dielectric  or 
ferrite  half  planes,  solutions  based  on  the  Wiener-Hopf  technique  have 
been  obtained  where  one  initially  approximates  the  effect  of  the  thin 
dielectric  or  ferrite  slab  by  an  equivalent  impedance  boundary 
condition.  The  Wiener  Hopf  solutions  to  the  latter  canonical  problems 
dealing  with  the  two  part  Impedance  surface  were  briefly  mentioned  in 
the  previous  part  (b)  dealing  with  discontinuities  In  surface  Impedance. 
The  present  analysis  begins  by  bisecting  the  semi -Infinite  dielectric  or 
ferrite  half  plane  by  an  electric  wall  In  the  first  case,  and  by  a 
magnetic  wall  in  the  second  case.  The  problem  of  plane  (or  surface) 
wave  diffraction  by  the  dielectric  or  ferrite  half  plane  Is  then 
constructed  by  appropriately  superimposing  the  corresponding  solutions 
for  the  electric  and  magnetic  wall  bisections,  respectl vely.  This 
procedure  Is  expected  to  yield  a  dielectric  half-plane  diffraction 
coefficient  which  Is  more  accurate  than  the  one  obtained  by  Anderson 
for  the  case  of  an  incident  electric  plane  wave  field  which  is  parallel 
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to  the  edge  of  the  thin  dielectric  half  plane  [8],  because  the  latter 
analysis  employs  an  approximate  "equivalent”  polarization  current 
sheet  model  for  the  thin  dielectric  half  plane.  The  approximation  in 
[8]  contains  only  a  part  of  the  information  present  in  the  more  general 
approach  being  employed  in  our  work;  consequently,  it  is  found  that  the 
previous  analysis  in  [8]  yields  a  diffraction  coefficient  which  is  valid 
only  for  an  extremely  thin  dielectric  half  plane.  Furthermore,  the 
equivalent  polarization  current  approximation  leads  to  a  rather 
complicated  Wiener-Hopf  analysis  when  the  magnetic  field  is  parallel  to 
the  edge;  the  latter  case  has  not  been  treated  by  Anderson  [8].  It  is 
also  noted  that  the  Wiener-Hopf  factors  for  the  case  treated  by  Anderson 
[8]  do  not  appear  to  be  well  behaved  for  near  edge-on  plane-wave 
incidence.  In  contrast,  the  Wiener-Hopf  factors  being  employed  in  our 
work  are  based  on  Weinstein's  factorization  procedure  [7]  which  appears 
to  overcome  the  difficulties  present  in  [8],  Furthermore,  the  present 
work  includes  both  dielectric  as  well  as  ferrite  half  planes. 

At  the  present  time,  the  diffraction  coefficients  for  the 
two-dimensional  case  of  both  TE  and  TM  plane  and  surface  wave  excitation 
of  the  thin  dielectric  and  ferrite  half-planes  have  been  obtained, 
and  they  have  been  tested  for  accuracy  for  the  dielectric  case.  Also, 
these  results  have  been  extended  to  treat  a  moderately  thick  dielectric 
or  ferrite  half  plane  which  is  excited  not  only  by  plane  and  surface 
wave  fields,  but  also  by  a  cylindrical  wave  (or  line  source  excited) 
field.  The  thickness  is  allowed  to  be  such  that  only  the  first  TE  and 


TM  surface  wave  can  exist.  The  Uniform  GTD  (UTD)  solution  obtained  for 
the  thin  dielectric  or  ferrite  half  plane  suggests  an  ansatz  upon  which 
the  aforementioned  extensions  to  the  moderately  thick  dielectric  or 
ferrite  half  plane,  and  also  to  cylindrical  wave  illumination  are  based. 
In  addition  to  extending  the  2-D  analysis  to  treat  moderately  thick 
dielectric  and  ferrite  half  planes,  work  has  also  been  completed  to 
include  loss  In  the  dielectric  and  ferrite  half  planes.  The  final 
result  for  the  lossy  case  Is  almost  as  simple  to  use  as  the  lossless 
case. 

A  typical  numerical  result  for  the  diffraction  of  a  cylindrical 
wave  by  a  moderately  thick  dielectric  strip  that  Is  based  on  the  present 
UTD  solution  is  shown  in  Figure  A-li.  The  UTD  based  result  in  Figure 
A-ll  is  compared  with  an  independent  numerical  moment-method  solution  of 
an  integral  equation  pertaining  to  this  line  source  excited  lossy 
dielectric  strip  geometry.  Since  the  line  source  illuminates  the  strip 
at  near  grazing  angles  of  incidence  in  Figure  A-ll,  it  is  especially 
important  in  this  near  grazing  incidence  case  to  include  surface  wave 
diffraction  effects  as  well.  The  contribution  of  all  the  dominant 
singly  and  multiply  diffracted  rays  is  included  for  calculating  the 
pattern  in  Figure  A-ll.  It  is  noted  from  the  pattern  in  Figure  A-ll 
that  the  total  UTD  field  is  continuous  at  the  reflection  and 
transmission  shadow  boundaries  as  it  should  be;  furthermore,  the  very 
close  agreement  between  the  totally  independent  UTD  and  moment-method 
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solution  In  Figure  A-ll  (and  also  In  other  cases  which  are  not  shown 
here)  Is  Indeed  gratifying.  The  effects  of  surface  wave  diffraction 
become  more  pronounced  In  the  lossless  case. 

Work  Is  currently  In  progress  to  treat  the  three-dimensional  (3-D) 
problems  of  the  diffraction  of  obliquely  Incident  plane,  conical,  and 
spherical  waves,  and  obliquely  Incident  surface  waves  by  a  thin  (or  at 
most  moderately  thick)  dielectric  half-plane.  These  Important 
generalizations  are  not  trivial  because  an  obliquely  Incident  plane  wave 
is  expected  to  launch  both  TE  and  TM  type  surface  waves  at  the  edge  of 
the  dielectric  half  plane.  In  other  words,  mode  coupling  between  TE  and 
TM  surface  waves  obliquely  Incident  on  the  edge  of  a  dielectric  half 
plane  Is  possible.  Consequently,  the  generalization  of  the 
two-dimensional  (2-D)  solutions  to  treat  the  corresponding 
three-dimensional  (3-D)  problems  Is  not  straightforward.  It  Is  also 
noted  that  even  though  a  relationship  exists  between  the  2-D  and  3-D 
solutions,  this  Is  not  a  direct  relationship  of  the  type  which  would 
allow  the  3-D  solutions  to  be  constructed  In  a  simple  way  from  the  2-D 
solutions  as  Is  possible  for  the  case  of  a  perfectly  conducting  half 
plane.  The  3-D  dielectric  half  plane  diffraction  solution  may  also  be 
of  Interest  In  the  area  of  Integrated  optics  and  millimeter  wave 
Integrated  circuits  where  open  dielectric  guiding  structures  are 
employed. 


2.  Diffraction  by  Perfectly-Conducting  Surfaces 
a.  Vertex  diffraction 

The  analysis  of  vertex  diffraction  Is  Important  because  one 
frequently  encounters  situations  where  an  antenna  radiates  In  the 
presence  of  planar  structures  with  edges  which  terminate  In  a  vertex  (or 
corner).  Also,  flat  plates  with  edges  are  used  to  model  aircraft  wings 
and  vertical  or  horizontal  staballzers  for  calculating  the  patterns  of 
antennas  on  aircraft.  In  these  problems,  the  antenna  pattern  Is 
affected  both  by  edge  diffraction  and  by  the  diffraction  at  the 
vertices. 

A  typical  vertex  In  a  planar,  perfectly-conducting  surface  Is  shown 
In  Figure  A-12.  A  more  general  vertex  In  a  planar  surface  Is  formed  by 
the  Intersection  of  two  otherwise  smooth,  curved  edges.  The  angle  a  Is 
the  Internal  angle  formed  by  the  tangents  to  each  of  the  curved  edges  at 
the  vertex. 

The  asymptotic  high-frequency  analysis  of  electromagnetic  vertex 
diffraction  Is  rather  complicated.  Vertices  not  only  shadow  the 
Incident  field,  but  they  also  shadow  the  edge  diffracted  fields.  The 
shadow  boundary  of  an  edge  diffracted  field  Is  a  conical  surface  whose 
tip  coincides  with  the  vertex  and  whose  axis  Is  an  extension  of  the 
shadowed  edge.  The  vertex  Introduces  a  diffracted  ray  which  penetrates 
the  shadow  regions;  moreover,  the  vertex  diffracted  field  must  also 
compensate  the  discontlnutles  In  the  incident  and  edge  diffracted 
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Figure  A-12.  Various  rays  associated  with  the  reflection  and 
diffraction  by  a  plane  angular  sector. 

fields  at  their  shadow  boundaries.  At  these  boundaries  the  vertex 
diffracted  field  assumes  Its  largest  magnitude  and,  hence.  Its  greatest 
Importance.  If  the  vertex  diffracted  field  is  omitted  in  the  GTD 
solution,  then  substantial  discontinuities  connected  with  the  shadowing 
of  the  Incident  and  edge  diffracted  fields  may  occur  In  the  calculated 
radiation  pattern. 

A  simple,  approximate  vertex  diffraction  coefficient  which  appears 
to  work  reasonably  well  In  certain  cases  has  been  obtained  at  the 
ElectroSclence  Laboratory  [9,10].  However,  this  result  has  been  found 
heurlstlcally  using  a  combination  of  theory  and  experiment;  It  therefore 
needs  to  be  Improved  In  order  to  be  useful  In  the  general  situations 


encountered  In  practice.  Nevertheless,  this  diffraction  coefficient 
offers  some  clues  for  constructing  the  more  refined  and  useful  vertex 
diffraction  coefficient,  which  we  expect  to  obtain  from  asymptotic 
analysis. 

The  objective  of  this  research  is  a  uniform  asymptotic  solution  for 
the  electromagnetic  field  scattered  by  the  vertex,  which  can  be 
Interpreted  ray  optically  In  terms  of  contributions  from  the  critical 
points  Qr,  Qi,  Q2  and  Qy.  The  approach  Is  to  approximate  the  current 
away  from  the  edges  by  the  geometrical  optics  current.  Near  the  edges, 
but  away  from  the  vertex,  the  local  half -plane  current  Is  used  and  near 
the  vertex  the  first  few  terms  of  the  eigenfunction  solution  [11]  are 
employed.  The  Integral  representation  of  the  field  radiated  by  these 
currents  Is  then  approximated  by  the  stationary  phase  method  to  obtain 
the  fields  diffracted  and  reflected  from  the  critical  points. 

Recently  a  method  for  asymptotically  evaluating  the  field  of  the 
current  near  an  edge  was  Investigated;  an  eigenfunction  representation 
of  the  edge  current  was  used.  The  preliminary  results  are  encouraging 
and  If  the  method  Is  successful.  It  will  be  applied  to  the  eigenfunction 
representation  of  the  current  near  the  vertex  to  obtain  the  vertex 
diffraction  coefficient. 

During  the  past  year,  the  confluence  of  the  critical  points  was 
examined  using  the  geometrical  optics  current  everywhere  on  the  vertex. 
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When  all  the  critical  points  are  close  to  each  other,  then  the 
expression  obtained  for  the  transition  function  Is  so  complicated  It 
does  not  appear  to  be  useful.  An  effort  Is  being  made  to  simplify  this 
expression. 

As  was  mentioned  earlier,  a  convergent  (eigenfunction)  solution  Is 
needed  to  represent  the  current  near  the  vertex.  Also,  a  convergent 
solution  Is  of  great  value  In  numerically  checking  the  approximate 
high-frequency  solutions  obtained  by  asymptotic  methods.  Therefore,  an 
effort  has  been  made  to  accurately  determine  the  dyadic  Green's  function 
for  the  plane  angular  sector. 

To  find  the  dyadic  Green's  function  we  begin  by  expanding  it  In 
terms  of  a  complete  set  of  vector  wave  functions  which  are  solutions  of 
the  vector  wave  equation  along  with  the  radiation  condition  and  the 
boundary  conditions  on  the  surface  containing  the  vertex.  The  vector 
wave  functions.  In  turn,  are  expressed  in  terms  of  scalar  wave  functions 
which  are  solutions  to  the  scalar  wave  equation  with  the  appropriate 
boundary  conditions.  Both  Neumann  and  Dlrlchlet  type  boundary 
conditions  must  be  satisfied  to  yield  a  complete  set  of  vector  wave 
functions.  The  final  step  of  the  solution  Involves  separating  the 
scalar  wave  equation  In  the  sphero-conal  coordinate  system.  The 
resulting  separated  equations  Include  the  spherical  Bessel  equation  and 
two  Lame*  equations  (one  with  periodic  boundary  conditions  and  the  other 
with  nonperiodic  ones)  which  are  coupled  through  the  two  eigenvalues 
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which  are  actually  the  separation  constants.  It  Is  precisely  these 
eigenvalue  pairs  which  serve  as  the  summation  Index  of  the 
dyadic  Green's  function. 

The  solution  Is  thus  ultimately  reduced  to  solving  for  the 
eigen-values  and  eigenfunctions  of  the  separated  Lame1  equations.  It  Is 
then  a  straightforward  procedure  to  construct  the  vector  wave  functions 
and  hence  the  dyadic  Green's  function.  Once  this  Is  found  one  can 
proceed  to  Investigate  a  wide  variety  of  problems  because  of  the 
versatility  and  general  nature  of  the  Green's  function  solutions. 

There  Is  no  known  closed  form  solution  to  the  Lame'  equations  and 
one  Is  therefore  attracted  to  an  Infinite  series  solution.  Earlier  work 
[11]  used  Fourier  sine  and  cosine  series  representations,  but  these 
resulted  In  the  need  to  solve  two  simultaneous.  Infinite 
continued-fraction  equations  for  the  eigenvalues  and  eigenvectors.  This 
solution  proves  to  be  numerically  formidable,  and  indeed,  almost 
Impossible  for  higher  order  eigenvalues  because  of  the  rapidly  varying 
nature  of  the  continued  fractions.  However,  Sahalos  has  used  these 
representations  [12]  and  reports  that  he  has  developed  methods  of 
calculation  which  are  more  efficient  and  accurate.  We  plan  to  study  his 
methods  before  doing  further  work  on  this  problem;  It  may  be  adequate 
for  purpose  of  calculating  the  dyadic  Green's  function  for  the  plane 
angular  sector. 


During  the  past  two  years  we  have  worked  on  the  special  case  of  the 
dyadic  Green's  function  where  o  =  *  and  the  plane  angular  sector  reduces 
to  a  half  plane.  This  Is  discussed  in  the  following  section. 
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b.  Paraxial  diffraction 


1)  At  edges 

As  was  mentioned  in  the  last  section,  the  dyadic  Green's  function 
for  the  plane  angular  sector  can  be  applied  to  the  half-plane  by  letting 
the  sector  angle  go  to  180°.  The  Lame'  differential  equations  decouple 
and  the  eigenvalues  and  eigenfunctions  are  easily  found;  furthermore, 
this  solution  can  be  extended  to  the  wedge  without  difficulty.  The 
advantage  of  this  spherical  wave  solution  for  the  dyadic  Green's 
function  of  the  wedge  over  that  given  In  [13]  is  that  It  can  be  used 
when  the  field  and  the  source  point  are  close  to  the  edge.  This  should 
make  It  possible  to  extend  the  UTD  solution  to  the  paraxial  region, 
l.e.,  to  grazing  Incidence  on  the  edge,  and  to  study  the  propagation  of 
edge  waves.  These  waves  are  Important  in  analyzing  the  scattering  from 
plates  Illuminated  at  or  near  grazing  Incidence  or  In  treating  the 
coupling  between  monopoles  or  slots  near  the  edge  of  a  wedge. 

During  the  past  year,  a  paper  was  written  on  the  spherical  wave 
representation  of  the  dyadic  Green's  function  for  the 
perfectly-conducting  wedge  [14];  Its  use  to  determine  the  fields  In  the 
paraxial  region  was  also  described.  A  summary  of  this  paper  was 
presented  at  the  1983  International  URSI  Symposium  [15]. 

The  use  of  this  dyadic  Green's  function  to  treat  the  scattering 
from  an  object  at  the  edge  of  the  wedge  was  described  In  the  preceding 
annual  report;  again  the  paraxial  region  was  of  special  Interest.  A 


T-Matrlx  method  was  proposed  for  non-spherical  objects.  In  the  period 
covered  by  this  report,  the  method  was  tested  by  choosing  a  sphere  whose 
center  is  displaced  along  the  edge  from  the  origin  of  the  coordinate 
system.  This  simulates  the  difficulties  encountered  in  treating 
a  non-spherical  scatterer  by  this  method. 

Calculations  of  the  echo  area  a  (normalized  with  respect  to 
wavelength  X)  are  shown  In  Figures  A-13a,b.  The  direction  of  incidence 
makes  an  angle  of  30°  with  respect  to  the  edge  of  the  wedge  (half-plane 
case),  the  radius  of  the  sphere  Is  0.25\,  and  its  center  is  displaced 
0 . 125 X  from  the  origin  of  the  coordinate  system.  The  incident  and 
scattered  electric  fields  are  polarized  in  the  n-di recti  on  in  Figure 
A-13a  and  in  the  ^-direction  In  Figure  A-13b.  The  echo  area  curves  are 
seen  to  be  in  excellent  agreement  with  those  calculated  from  an  (exact) 
eigenfunction  solution.  This  establishes  confidence  in  our  T-Matrix 
solution. 

Recently  some  work  has  been  done  on  the  problem  of  multiple 
scattering  where  there  are  two  or  more  objects  arrayed  along  the  edge  of 
the  wedge. 


Figure  A-13a.  Backscatter  from  a  spherical  boss  at  the  edge  of  a  half 
plane.  Incident  field  polarized  in  the  e-directlon. 


Figure  A-13b.  Backscatter  from  a  spherical  boss  at  the  edge  of  a  half 
plane.  Incident  field  polarized  In  the  ^-direction. 


11)  At  smooth,  quasl-cyllndrlcal  surfaces 


Several  papers  were  written  (under  the  JSEP  program)  and  published 
[1,2,3]  which  describe  Uniform  GTO  (UTD)  solutions  for  the  diffraction 
by  perfectly-conducting  convex  surfaces.  In  particular,  efficient  UTD 
solutions  for  the  problems  of  the  radiation  from  sources  both  off  and  on 
a  convex  surface  and  the  mutual  coupling  between  sources  on  a  convex 
surface  were  presented  in  [1,2,3].  These  UTD  solutions,  for  the 
problems  of  scattering,  radiation,  and  mutual  coupling,  which  are 
associated  with  the  radiation  by  antennas  in  the  presence  of  an 
arbitrary,  smooth  perfectly-conducting  convex  surface  that  Is 
electrically  large  represent  an  important  and  useful  contribution  to  the 
area  of  ray  methods  for  analyzing  the  EM  radiation  and  scattering  from 
complex  structures.  It  is  noted  that  the  effects  of  surface  ray  torsion 
on  the  diffracted  fields  are  explicitly  identified  in  these  solutions. 
Here,  the  diffracted  fields  are  associated  with  surface  rays  as  well  as 
with  rays  shed  from  the  surface  rays.  It  is  noted  that  these  surface 
rays  on  a  convex  surface  traverse  geodesic  paths  which  in  general  are 
torsional;  i.e.,  the  surface  ray  paths  are  twisted  (or  they  do  not  lie 
in  a  plane) . 

While  the  above  mentioned  UTD  solutions  for  sources  on  or  off  a 
smooth  perfectly-conducting  convex  surface  are  valid  under  very  general 
conditions,  they  must  be  modified  within  the  paraxial  regions.  For 
example,  these  solutions  must  be  modified  for  an  observation  point 


(either  on  or  off  the  surface)  which  lies  In  the  paraxial  region  of  an 
elongated  or  cigar  shaped  (quasl-cyl indrical )  convex  surface  whenever 
the  rays  from  the  source  to  that  observation  point  traverse  paths  which 
lie  within  the  paraxial  zone.  At  the  present  time,  the  solution  for  the 
surface  fields  of  a  source  on  the  convex  surface  has  been  extended  so  as 
to  be  valid  within  the  paraxial  region  of  quasi -cyl indrical ,  or 
elongated  convex  surfaces.  The  surface  field  (or  mutual  coupling 
related)  solution  provides  the  electric  current  density  which  is  induced 
by  a  source  on  the  same  surface.  Such  a  surface  current  density  can  be 
incorporated  into  the  usual  radiation  integral  to  find  the  field 
radiated  by  this  current;  that  step  has  also  been  performed  presently. 
Thus,  the  field  radiated  by  the  source  on  the  surface  has  been  obtained 
via  an  asymptotic  high  frequency  evaluation  of  the  radiation  Integral 
containing  this  surface  current  density  which  remains  valid  in  the 
paraxial  zone.  The  asymptotic  evaluation  of  this  radiation  integral 
associated  with  a  general  convex  surface  has  been  done  carefully  in 
order  to  obtain  a  simple  and  useful  uniform  solution  for  the  radiated 
field.  It  may  be  remarked  that  previous  attempts  by  others,  which  have 
been  reported  earlier,  were  not  very  successful  in  asymptotically 
evaluating  such  radiation  integrals  associated  with  a  general  convex 
surface.  We  have  found  that  such  previous  attempts  did  not  lead  to  very 
useful  results  because  of  a  simplifying  assumption  which  is  inherent  in 
their  work.  Our  present  analysis  reveals  that  such  an  asymptotic 
analysis  must  be  done  by  a  different  approach  which  circumvents  the 


Inherent  approximation  present  In  the  previously  reported  approaches  to 
obtain  meaningful  results.  The  solution  to  the  problem  of  EM  radiation 
in  the  paraxial  zone  by  sources  on  a  convex  surface  is  therefore 
currently  nearing  completion.  When  the  work  on  the  radiation  problem  is 
completed,  It  will  be  extended  subsequently  to  deal  with  the  problem  of 
EM  scattering  within  the  paraxial  regions  In  almost  the  same  manner  as 
the  radiation  solution  was  obtained  via  an  extension  of  the  surface 
field  (mutual  coupling  related)  solution  through  the  use  of  the 
radiation  Integral. 

c.  Slope  diffraction 

If  the  field  incident  at  the  edge  of  a  wedge  or  the  shadow  boundary 
of  a  convex  surface  has  a  rapid  spatial  variation,  a  slope  diffraction 
term  must  be  added  to  the  ordinary  UTD.  The  slope  diffraction  term  is 
proportional  to  the  spatial  derivatives  of  the  incident  field.  It 
ensures  that  the  spatial  derivatives  of  the  resulting  radiation  pattern 
are  continuous  at  the  shadow  and  reflection  boundaries  so  that  there  are 
no  “kinks"  In  the  calculated  pattern. 

The  derivation  of  slope  diffraction  terms  for  ordinary  wedges 
formed  by  plane  surfaces  and  for  convex  cylinders  has  been  described  in 
the  two  preceding  annual  reports,  and  a  dissertation  based  on  this  work 
has  been  written  [16].  During  the  present  period,  the  rough  draft  of  a 
paper  on  this  subject  has  been  written  and  there  has  been  an  effort  to 
extend  the  solution  to  a  curved  wedge  formed  by  curved  surfaces. 


d.  Diffraction  by  a  discontinuity  in  surface  curvature 


The  diffraction  by  an  edge  in  a  perfectly  conducting  surface 
formed  by  a  discontinuity  In  surface,  such  as  occurs  in  the  cone  sphere 
and  the  hemi spherically  capped  cylinder,  contributes  strongly  to  the 
radiation  in  the  vicinity  of  the  reflection  boundary  where  it 
compensates  for  the  discontinuity  in  the  reflected  field.  It  may  also 
contribute  strongly  to  the  axial  backscatter  from  a  body  of  revolution 
with  a  discontinuity  in  surface  curvature.  The  geometry  of  the  problem 
is  depicted  in  Figure  A-14. 

A  diffraction  coefficient  for  the  discontinuity  in  surface 
curvature  has  been  described  in  the  two  preceding  annual  reports.  The 
solution  is  restricted  to  the  case  where  the  direction  of  incidence  is 
normal  to  the  edge.  The  purpose  of  the  present  work  to  obtain  a  simpler 
expression  for  this  diffraction  coefficient  and  to  generalize  it  to  the 
oblique  incidence  case. 

It  has  been  noted  by  Chu  [17]  that  the  curved  wedge  diffraction 
coefficient  obtained  by  Kouyoumjian  and  Pathak  [18]  can  be  used  to 
calculate  the  diffraction  by  a  discontinuity  in  surface  curvature  in  the 
vicinity  of  the  reflection  boundary.  Although  the  first  two  terms  In 
the  curved  wedge  diffraction  coefficient  vanish  in  this  case,  the  second 
pair  of  terms  associated  with  the  reflected  field  do  not.  They 
correctly  compensate  the  discontinuity  in  the  reflected  field  at  the 
reflection  boundary.  However,  outside  the  transition  region  of  the 
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reflection  boundary,  the  diffracted  field  obtained  by  this  method  does 
not  reduce  to  the  correct,  non-uniform  asymptotic  form  given  by  Senior 

[19] .  Starting  with  the  expression  for  the  curved  wedge  diffraction 
coefficient,  a  diffraction  coefficient  for  the  discontinuity  In 
curvature  was  derived,  which  reduces  to  the  solution  given  by  Senior 
away  from  the  reflection  boundary  and  exactly  compensates  the 
discontinuity  In  the  reflected  field  at  the  reflection  boundary.  It  can 
be  used  at  oblique  Incidence  and  is  computationally  simpler  than  the 
diffraction  coefficient  described  In  the  1982  annual  report.  Away  from 
grazing  Incidence  and  diffraction,  the  two  diffraction  coefficients  give 
virtually  the  same  numerical  results.  This  Is  demonstrated  in  Figures 
A-15a  and  A-15b  where  the  soft  and  hard  diffraction  coefficients, 
respectively,  are  plotted  as  a  function  of  4».  In  these  figures,  the 
direction  of  Incidence  Is  normal  to  the  edge,  $’=60°,  ao=10A,  atfp5A, 
and  the  field  point  Is  In  the  far  zone. 

Since  the  new  diffraction  coefficient  Is  based  on  the  curved  edge 
diffraction  coefficient.  It  can  not  be  used  when  the  direction  of 
Incidence  or  diffraction  grazes  the  surface  unless  the  surface  In 
question  Is  planar.  It  Is  expected  that  this  restriction  can  be  removed 
by  Incorporating  the  generalized  reflection  coefficient  given  by  Pathak 

[20]  Into  the  solution. 


3.  Diffraction  of  Non-Ray  Optical  Fields  (Source  Close  to 
An  Edge) 

In  the  conventional  form  of  the  Uniform  GTD,  it  Is  assumed  that  the 
Incident  field  Is  a  ray-optical  field,  which  Implies  that  It  Is 
polarized  in  a  direction  perpendicular  to  the  Incident  ray.  In  general, 
this  requires  that  the  source  of  the  Incident  field  be  sufficiently  far 
from  the  point  of  diffraction  so  that  the  component  of  the  incident 
field  parallel  to  Its  ray  path  (the  component  in  the  radial  direction 
from  the  source)  Is  negligible  at  the  diffraction  point.  However,  In 
some  applications  this  Is  not  the  case,  e.g.,  a  monopole  antenna  may  be 
mounted  at  or  very  close  to  the  edge  of  a  ship  or  the  edges  of  wings  and 
stabilizers.  This  case  Is  also  of  Interest  In  the  development  of  the 
Hybrid  GTD /Moment  Method  solution,  where  It  Is  desired  to  calculate  the 
Input  Impedance  of  a  wire  antenna  close  to  an  edge. 

An  asymptotic  solution  for  the  diffraction  of  the  fields  of 
electric  and  magnetic  dipoles  close  to  the  edge  of  a  wedge  has  been 
obtained.  The  basic  approach  is  to  determine  the  field  close  to  the 
edge  of  the  wedge;  the  far  field  of  the  dipole  Is  then  found  using 
reciprocity.  The  asymptotic  analysis  of  the  diffracted  field  proceeds 
as  done  earlier  [18]  except  that  higher  order  terms  are  retained.  The 
diffraction  coefficient  has  the  form  of  a  3x2  matrix  which  contains 
radial  and  cross-polarized  components.  This  solution  has  been  used  to 
accurately  calculate  the  fields  of  dipoles  which  are  only  a  few  tenths 
of  a  wavelength  from  the  edge  [21].  In  the  present  period,  this  result 


has  been  Improved  further  and  the  components  of  the  diffraction 


coefficient  have  been  expressed  more  compactly.  This  work  has  just 


recently  been  extended  to  the  case  of  spherical  wave  illumination  of  the 


wedge. 


The  dyadic  Green's  function  described  in  section  2.b  converges 


rapidly  when  the  source  point  is  close  to  the  edge  and  so  it  is 


convenient  to  employ  it  to  numerically  check  the  solutions  described  in 


the  preceding  paragraph;  extensive  comparisons  of  this  type  have  been 


carried  out.  It  is  hoped  that  the  two  types  of  solutions  for  edge 


diffraction  can  be  combined  so  that  we  will  have  a  more  useful 


computational  algorithm  for  edge  diffraction.  This  is  particularly 


important  in  the  paraxial  region  where  the  high  frequency  solutions 


described  in  the  last  paragraph  lose  accuracy. 
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4.  Caustic  Field  Analysis 


The  GTO  Is  a  very  convenient  and  accurate  procedure  for  analyzing 
high  frequency  radiation,  scattering,  and  diffraction  problems. 

However,  the  GTD  suf fe  s  from  a  limitation  Inherent  in  ray  methods; 
namely.  It  cannot  be  employed  directly  to  evaluate  fields  at  and  near 
focal  points  or  caustics  of  ray  systems.  The  fields  at  caustics  must, 
therefore,  be  found  from  separate  considerations  [22,23]. 

In  certain  problems  such  as  In  the  diffraction  by  smooth,  closed 
convex  surfaces  or  by  surfaces  with  a  ring-type  edge  discontinuity,  it 
Is  possible  to  employ  the  GTO  indirectly  to  evaluate  the  fields  In  the 
caustic  regions  via  the  equivalent  ring  current  method  [24,25]. 

However,  even  the  equivalent  ring  current  method  falls  If  the  incident 
or  reflection  shadow  boundaries  are  near  or  on  a  caustic. 

The  recently  developed  uniform  GTD  (or  UTD)  solution  for  the 
scattering  and  diffraction  of  waves  by  a  convex  surface  [1,2,3]  offers 
clues  as  to  how  It  may  be  employed  indirectly  to  obtain  the  far  zone 
fields  in  caustic  regions  where  the  surface  is  Illuminated  by  a  distant 
source.  In  the  latter  case,  the  shadow  boundary  and  caustic  transition 
regions  tend  to  overlap.  The  far  zone  fields  In  the  near  axial 
direction  of  a  closed  surface  of  revolution  Illuminated  by  an  axially 
directed  plane  wave  can  be  expressed  in  terms  of  an  equivalent  ring 
current  contribution  plus  a  dominant  term  which  may  be  Interpreted  as  an 
"effective  aperture  Integral".  The  latter  Integral  can  be  evaluated  in 
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closed  form.  In  the  near  zone,  where  the  shadow  boundary  and  caustic 
directions  are  sufficiently  far  apart,  only  the  equivalent  ring  current 
contribution  must  remain  significant.  Presently,  such  a  solution  has 
been  developed  for  the  special  case  of  axial  Incidence  on  a  rotationally 
symmetric  convex  surface  (or  revolution).  The  generalization  of  that 
solution  to  treat  non-axlal  Incidence  and  also  closed  convex  surfaces 
which  are  not  necessarily  surfaces  of  revolution  forms  the  subject  of 
future  Investigation. 

Another  related  and  interesting  problem  which  is  presently  under 
study  involves  the  analysis  of  high  frequency  electromagnetic 
diffraction  by  perfectly  conducting  planar  surfaces  (or  plates)  with  a 
smooth  convex  boundary.  In  this  problem,  the  planar  face  of  the  plate 
gives  rise  to  caustic  effects.  For  certain  aspects,  the  caustic  and 
specular  reflection  directions  can  coincide  so  that  a  direct  use  of  UTD 
becomes  Invalid.  The  present  study  Is  aimed  at  determining  a  proper 
modification  of  the  UTD  so  that  a  new  and  efficient  equivalent  current 
method  can  be  used  In  these  situations.  Away  from  the  caustic  and 
confluence  of  caustic  and  specular  directions,  the  modified  UTD  solution 
must  reduce  to  the  conventional  GTD  solution.  More  related  problems 
Involving  a  confluence  of  diffracted  ray  caustic  and  geometrical  optical 
shadow  boundary  transition  regions  will  be  Investigated  in  the  future 
phases  of  this  work. 
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5.  Electromagnetic  Dyadic  Green's  Function 

Electromagnetic  dyadic  Green's  functions  are  useful  in  that  they 
allow  one  to  calculate  the  fields  of  an  arbitrary  source  distribution  in 
a  systematic  fashion.  In  the  case  of  special  Green's  dyadics  which 
account  for  the  source  radiation  in  the  presence  of  any  boundaries, 
i.e.,  structures,  they  form  the  starting  point  for  the  asymptotic  high 
frequency  analysis  of  canonical  problems  in  the  uniform  geometrical 
theory  of  diffraction  (or  UTD).  In  the  latter  application,  it  is 
implied  that  a  Green's  dyadic  for  the  canonical  geometry  is  known  in  the 
first  step;  this  Greens'  dyadic  is  then  generally  transformed  into  a 
useful  form  in  terms  of  appropriate  integrals  if  it  is  not  originally 
available  in  that  form  (since  it  could  originally  be  given  as  a  series 
expansion),  and  finally,  the  intgrals  are  evaluated  asymptotically  to 
extract  the  high  frequency  UTD  solution  of  interest. 

A  relatively  simple  method  has  been  developed  recently  for 
constructing  the  complete  eigenfunction  expansion  of  time-harmonic 
electric  (E)  and  magnetic  (H)  fields  within  exterior  or  -'nterior  regions 
containing  an  arbitrarily  oriented  electric  current  point  source.  In 
particular,  since  these  fields  E  and  H  pertain  to  the  case  of  point 
source  excitation,  they  yield  directly  the  complete  eigenfunction 
expansion  of  the  corresponding  electric  and  magnetic  dyadic  Green's 
functions  Ge  and  Gm  that  are  associated  with  E  and  H,  respectively. 


The  expansion  of  Ge  and  obtained  here  Is  expressed  in  terms  of 
only  the  solenoidal  type  eigenfunctions  which  satisfy  the  required 
boundary  and  radiation  conditions.  In  addition,  the  expansion  of  Gg 
also  contains  an  explicit  dyadic  delta  function  term  which  is  required 
for  making  that  expansion  complete  at  the  source  point.  This  explicit 
dyadic  delta  function  term  in  ^  is  obtained  readily  in  the  present 
approach  from  a  simple  condition  that  the  eigenfunction  expansion  must 
satisfy  at  the  source  point,  provided  one  interprets  that  condition  in 
the  sense  of  distribution  theory. 

A  paper  has  been  wtitten  on  the  above  topic,  namely 

"On  the  Eigenfunction  Expansion  of  Electromagnetic  Dyadic 
Green's  Functions",  by  P.H.  Pathak,  IEEE  Transactions  on 
Antennas  and  Propagation,  Af P-31,  No.  6,  pp.  837-846, 

November  1983. 
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A  General  Description  of  Hybrid  Techniques 

The  method  of  moments  (MM)  provides  a  means  of  solving 
electromagnetic  boundary  value  problems  in  terms  of  a  set  of 
simultaneous  linear  equations.  In  general,  the  electromagnetic  boundary 
value  problem  is  formulated  as  an  integral  equation  for  the  unknown 
surface  fields  on  the  antenna  or  scatterer  and  the  integral  equation  is 
then  reduced  to  a  system  of  equations  by  expanding  the  unknown  in  terms 
of  a  basis  set  and  by  enforcing  this  expansion  to  satisfy  the  boundary 
conditions  in  some  average  sense  through  the  use  of  testing  functions. 
However,  the  MM  procedure  can  become  inefficient  and  cumbersome  if  the 
number  of  unknowns  (coefficients  of  the  expansion  or  basis  functions) 
becomes  large  as  is  the  case  for  antennas  or  scatterers  which  are  not 
small  in  terms  of  wavelength.  On  the  other  hand,  the  geometrical  theory 
of  diffraction  (GTO)  exploits  the  local  nature  of  high  frequency  wave 
propagation,  diffraction  and  radiation,  thereby  reducing  the  antenna 
radiation  or  scattering  problem  to  calculating  the  fields  associated 
with  just  a  few  rays  emanating  from  edges,  tips,  and  shadow  boundaries 
(of  smooth  convex  surfaces),  and  also  from  other  discontinuities  in  the 


geometrical  and  electrical  properties  of  the  antenna  or  scatterer. 
Although  the  GTD  Is  a  high  frequency  technique.  It  works  rather  well, 
even  for  structures  which  are  only  moderately  large  In  terms  of  the 
wavelength.  However,  the  use  of  the  GTD  Is  limited  by  the  number  of 
available  diffraction  coefficients  for  characterizing  a  particular  type 
of  electrical  and/or  geometrical  discontinuity.  It  Is  obvious  that  a 
procedure  Is  desirable  which  would  overcome  the  limitations  of  the 
Individual  MM  and  GTD  approaches.  Such  a  procedure,  referred  to  as  the 
"hybrid"  GTD-MM  procedure,  can  Indeed  overcome  the  limitations  of  the 
Individual  W  and  GTD  approaches  by  actually  combining  the  best  features 
of  both  methods.  In  particular,  GTD  provides  the  form  of  the  local 
field  over  any  part  of  the  antenna  or  scattering  structure,  which  Is  at 
least  moderately  large  In  terms  of  the  wavelength;  hence,  the  form  of 
the  GTD  field  could  be  viewed  as  a  set  of  basis  functions  for  the 
expansion  of  the  unknowns  In  the  MM  formulation.  The  unknown 
coefficient  associated  with  this  type  of  GTD  basis  or  expansion 
functions  Is  then  the  diffraction  coefficient  for  the  surface  field 
calculations  If  the  unknown  In  the  Integral  equation  happens  to  be  the 
surface  field.  Thus,  by  using  the  local  GTD  field  form  outside  the 
region  where  the  structure  Is  small  In  terms  of.  the  wavelength,  the 
number  of  unknowns  Is  thereby  vastly  reduced  in  the  MM  procedure.  The 
expansion  for  the  unknown  within  regions  (of  the  structure)  which  are 
small  In  terms  of  wavelength.  Is  of  course,  done  according  to  the 
conventional  MM  approach  (perhaps  using  a  subsectional  basis  set  such  as 
rectangular  pulses,  etc).  Clearly,  the  hybrid  GTD-MM  procedure  can 


solve  problems  far  more  efficiently  than  the  MM  procedure  as  the 
frequency  increases.  Also,  it  can  provide  a  useful  check  on  future 
diffraction  coefficients  as  and  when  they  become  available. 

While  the  hybrid  6T0-MM  procedure  will  in  general  be  employed  to 
obtain  diffraction  coefficients,  other  hybrid  techniques  which  combine 
ray  methods  and  numerical  methods  different  from  the  MM  procedure  will 
also  be  studied.  Thus,  in  a  broader  sense,  the  area  of  hybrid 
techniques  will  emphasize  useful  combinations  of  high  frequency  or  ray 
techniques  with  numerical  methods  for  solving  a  variety  of  Interesting 
and  useful  electromagnetic  radiation  and  scattering  problems. 

Presently,  two  classes  of  diffraction  problems  are  being 
investigated;  namely,  the  diffraction  by  perfectly-conducting  structures 
Involving  special  edge  or  tip  type  discontinuities,  and  the  diffraction 
by  discontinuities  associated  with  the  junction  of  dielectric  and 
perfectly-conducting  structures  which  are  connected  in  special  ways. 

In  particular,  the  problems  in  the  first  catagory  are  the  problems  of 
diffraction  by  a  discontinuity  in  surface  curvature  (e.g.,  the  smooth 
join  of  a  half  plane  and  a  cylinder);  by  an  edge-type  discontinuity  in 
an  otherwise  smooth  convex  surface;  by  the  tip  of  a  cone-type  structure; 
by  the  tip  of  a  semi -infinite  wire;  and  by  an  aperture  in  a  thick 
screen  (and  related  configurations).  The  problems  in  the  second 
catagory  are  the  problems  of  diffraction  by  a  truncated  dielectric  slab 
recessed  in  a  perfectly-conducting  structure  as  in  Figure  B-l(a),  by  the 
junction  of  a  dielectric  half  plane  which  is  joined  to  a  thick 
perfectly-conducting  half  plane  as  in  Figure  B-l(b);  and  the  problem 
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of  the  radiation  of  a  single  element  or  an  array  of  electrically  small 
microstrip  patch  antennas  which  are  placed  conformally  on  an 
electrically  large  perfectly-conducting  smooth  convex  surface  (such  as 
the  fuselage  of  an  aircraft  or  a  missile).  The  latter  problem,  as 
depicted  in  Figure  B-2,  requires  a  combination  of  low  frequency 
techniques  (for  handling  the  electrically  small  microstrip  patches)  with 
the  high  frequency  ray  technique  (for  handling  the  electrically  large 
convex  surface). 

Currently,  work  is  In  progress  to  obtain  an  efficient,  hybrid 
GTD-MM  solution  to  the  problem  of  the  diffraction  by  the  structure  shown 
In  Figure  B-l(a).  The  dielectric  Is  made  Inhomogeneous  near  the 
junction  of  the  dielectric  conducting  body  in  order  to  reduce  the 
surface  wave  reflection  from  that  junction.  A  computer  program  searches 
for  the  selection  of  dielectric  material  to  accomplish  this  purpose. 

Some  results  of  the  surface  wave  reflection  coefficient  and  the 
radiation  pattern  for  the  case  where  the  dielectric  Is  homogeneous  are 
shown  in  Figures  B-3  and  B-4.  Also  shown  In  these  figures  are  solutions 
of  the  same  problem  obtained  by  Pathak  [1]  using  a  Wlener-Hopf  type 
technique.  The  agreement  between  the  two  methods  Is  generally  very 
good.  One  should  note  that  the  technique  employed  by  Pathak  In  [1]  Is 
not  feasible  In  dealing  with  the  problem  where  the  dielectric  Is 
Inhomogeneous.  The  Inhomogeneous  dielectric  problem  Is  now  being 
Investigated  using  the  hybrid  GTD-MM  technique.  Both  surface  wave  and 
plane  wave  Incidences,  TE  and  TM  polarizations,  when  the  boundary 
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gure  B-3.  Surface  wave  reflection  coefficient  as  a  function  of  the 
dielectric  thickness.  Solid  curve  Is  reproduced  from 
reference  [1]. 
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Figure  B-4.  Radiation  pattern  of  the  structure  In  Figure  B-l(a)  with 
surface  wave  Incidence.  Solid  curve  Is  reproduced  from 
reference  [1]. 


beneath  the  dielectric  is  either  an  electric  conductor  or  a  magnetic 
conductor,  are  considered.  Combinations  of  these  general  solutions 
enable  one  to  solve  the  problem  posed  in  Figure  B-l(b).  These  solutions 
will  be  reported  and  published  in  the  future. 

A  hybrid  GTD-MM  solution  to  the  problem  of  an  exterior  line  source 
excitation  of  an  aperture  in  a  parallel  plate  wave  guide  as  seen  in 
Figure  B-5  was  also  initiated  during  this  past  year.  It  is  of  interest 
to  know  the  energy  scattered  from  and  penetrated  through  the  aperture. 

A  proper  choice  of  special  Green's  functions  in  the  integral  equation 
formulation  of  this  problem  requires  one  to  find  only  the  unknown 
electric  field  in  the  aperture.  The  aperture  field  can  be  obtained 
using  the  moment  method  technique  for  small  apertures  and  the  hybrid 
technique  for  large  apertures.  An  example  of  the  aperture  field  is 
shown  in  Figure  B-6  which  was  obtained  using  the  moment  method 
technique.  Extension  to  the  case  of  a  very  large  aperture  using  hybrid 
technique  will  be  evaluated  in  the  future.  The  fields  scattered  by  the 
aperture  Into  the  interior  and  exterior  regions  can  be  found  easily  once 
the  aperture  field  is  determined.  An  analyses  of  two  related  problems 
involving  the  EM  diffraction  by  an  aperture  in  a  thick  perfectly 
conducting  screen,  and  by  a  groove  In  a  planar  conducting  surface, 
respectively,  which  employed  a  combination  of  moment  method  and  multiple 
scattering  techniques,  was  mentioned  in  the  last  annual  report  [2],  A 
technical  report  and  a  paper  dealing  with  these  two  related  problems  is 
currently  in  preparation;  namely. 


Figure  B-5.  Exterior  line  source  excitation  of  a  parallel  plate 
waveguide  through  an  aperture  In  the  wall. 


"EM  Scattering  by  Slits  and  Grooves  in  Thick  Perfectly- 
Conducting  Planar  Surfaces",  by  R.  Kautz,  P.H.  Pathak,  and 
L.  Peters,  Jr., 

A  third  problem  of  interest  involving  the  EM  diffraction  by  a 
perfectly  conducting  half  {circular)  cylinder,  which  has  also  been 
investigated  during  the  past  year  is  illustrated  in  Figure  B-7.  The 
geometry  chosen  for  this  problem  is  one  of  the  simplest  examples  of  an 
edge  in  an  otherwise  smooth,  perfectly  conducting  finite  size  convex 
body.  A  hybrid  GTD-MM  technique  has  been  applied  to  analyze  the 
scattering  by  this  half-cylinder  configuration.  At  the  present  time, 
only  the  hybrid  GTD-MM  analysis  of  the  TE  (or  H  parallel  to  the  cylinder 
edge)  case  has  been  completed.  A  corresponding  analysis  of  the  TM  (or  E 
parallel  to  the  cylinder  edge)  case  is  currently  being  pursued.  This 
hybrid  solution  for  the  TE  case  is  very  efficient  since,  in  general, 
only  8-12  unknowns  are  required  to  represent  the  induced  surface  current 
in  the  moment  method  procedure  which  employs  the  mixed  (hybrid)  basis 
representation  for  the  current  in  terms  of  a  few  pulse  basis  functions 
for  the  current  around  the  edges,  and  the  uniform  GTD  type  basis 
functions  for  the  current  elsewhere  on  the  surface  of  the  half 
cylinder.  It  is  noted  that  the  above  number  of  unknowns  stays  the  same 
in  this  hybrid  GTD-MM  approach  even  if  the  electrical  size  of  the 
cylinder  is  increased;  in  contrast,  the  conventional  MM  approach  always 
requires  an  increase  in  the  number  of  unknowns  with  an  increase  in  the 
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Figure  B-7.  Diffraction  by  a  half  cylinder  of  radius  =  a. 

cylinder  size.  Some  typical  results  for  the  TE  radiation  patterns  of  an 
electric  current  line  source  in  the  presence  of  the  half  cylinder  of 
Figure  B-7  that  are  based  on  the  hybrid  GTD-MM  procedure  are  illustrated 
in  Figures  B-8  and  B-9  which  also  show  a  comparison  with  the 
corresponding  conventional  MM  solutions.  For  a  cylinder  whose  perimeter 
is  approximately  10A  (here  A=wavelength)  in  these  figures,  the  number  of 
unknowns  (using  a  pulse  basis  set)  in  the  conventional  MM  solution  is 
100;  whereas,  in  the  present  hybrid  GTD-MM  solution,  the  number  of 
unknowns  in  only  12.  A  paper  describing  this  work  is  planned  in  the 
near  future. 
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Figure  1 


-8.  Pattern  of  an  electric  line  source  in  the  presence  of  a 
perfectly  conducting  half  cylinder. 


rf  v-W  '  .  . 


The  half  cylinder  configuration  of  Figure  B-7  highlights  the 
effects  of  edge  excited  surface  rays.  It  is  planned  in  the  future 
phases  of  this  study  to  numerically  extract  the  uniform  GTD  diffraction 
coefficient  for  edge  excited  surface  rays  using  the  hybrid  GTD-MM 
procedure.  This  work  dealing  with  the  diffraction  by  an  edge  in  an 
otherwise  smooth  convex  surface  is  a  natural  and  proposed  extension  of 
the  work  reported  in  the  previous  years  on  the  development  of  a  hybrid 
GTD-MM  solution,  and  an  approximate  uniform  asymptotic  solution  to  the 
problem  of  the  diffraction  by  a  perfectly  conducting  convex  surface  with 
a  discontinuity  in  curvature.  A  paper  dealing  with  the  latter 
approximate  uniform  asymptotic  solution  is  currently  in  preparation: 

"Diffraction  by  a  Perfectly  Conducting  Surface  with  a 
Discontinuity  in  Surface  Curvature",  by  C.W.  Chuang. 
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Introduction 

The  goal  of  our  research  is  to  perform  basic  work  in  the  area  of 
integral  equation  solutions  of  general  or  "arbitrarily"  shaped 
structures.  It  is  anticipated  that  based  upon  this  work  general  purpose 
method  of  moments  (MM)  computer  codes  will  be  developed  (although  not 
under  JSEP  sponsorship)  for  the  analysis  of  "arbitrary"  structures.  In 
the  past,  we  have  studied  such  basic  problems  as  a  perfectly  conducting 
rectangular  plate,  the  junction  of  a  wire  to  a  plate  near  the  edge  of 
the  plate,  the  junction  of  several  plates,  and  the  polygonal  plate 
[1-6].  Much  of  this  basic  work  has  already  been  implemented  into 
general  purpose  MM  computer  codes  [7].  The  philosophy  of  these  codes  is 
to  model  the  structure  as  an  interconnection  of  wires  and  plates.  For 
example,  we  have  riodeled  the  Concorde  aircraft  as  an  interconnection  of 
14  polygonal  plates.  This  model  was  used  to  accurately  predict  the 
magnitude  (RCS)  and  phase  of  the  backscattering  from  the  Concorde  at  a 
frequency  where  it’s  overall  length  was  about  1.9  wavelengths  [8,9]. 

One  can  easily  envisage  how  the  polygonal  plates  could  be  interconnected 


to  model  other  aircraft,  ships,  buildings,  etc.  Thus,  the  computer  code 
Is  capable  of  modeling  a  very  general  structure.  Antennas  on  the 
structure  are  modeled  using  thin-wires. 

All  of  the  above  work  has  dealt  with  perfectly  conducting  plates. 
Thus,  It  was  suitable  only  for  the  modeling  of  (essentially)  perfectly 
conducting  structures.  However,  recently  there  has  been  Increasing 
Interest  in  the  use  of  non-metalllc  materials  for  the  construction  of 
aircraft.  These  non-metal  lie  materials,  often  called  composite 
materials,  are  used  to  modify  the  RCS  of  the  aircraft,  or  to  make  the 
aircraft  frame  strong  and  yet  lightweight.  It  is  our  future  goal  to  be 
able  to  model  a  general  structure  composed  of  metallic  and/or  composite 
materials.  While  the  modeling  of  the  perfectly  conducting  surfaces  Is 
well  advanced,  the  following  steps  are  required  to  realize  this  goal: 

1.  Oevelop  a  model  for  an  Isolated  material  plate.  A 
material  plate  is  a  thin  planar  three  dimensional  slab 
composed  of  some  dielectric/ferrite  material.  Also,  it 
can  be  a  perfectly  conducting  plate  coated  on  one  or  both 
sides  by  a  dielectric/ferrite  material.  This  can  be  used 
to  model  absorber  coated  plates. 

2.  Develop  a  model  for  the  junction  of  a  material  and  a 
perfectly  conducting  plate. 

3.  Extend  the  above  solutions  to  the  non-planar  case. 


4.  Implement  the  above  in  the  form  of  a  general  purpose 
computer  code  capable  of  treating  an  interconnection  of 
wires,  perfectly  conducting  plates,  and  material  plates. 

Step  1  is  now  essentially  complete  [20,21],  and  will  be  described 
in  detail  in  the  next  section.  This  will  include  a  description  of  the 
theoretical  solution  in  terms  of  for**  coupled  integral  equations,  the 
simplification  of  these  equations  using  a  surface  impedance 
approximation,  and  the  generation  of  numerical  results  using  the  method 
of  moments.  Steps  2  and  3  have  been  proposed  for  future  JSEP  research. 
Step  4  is  in  the  area  of  applied  research  or  development  and  must  be 
done  under  support  other  than  JSEP. 

1.  Theoretical  Solution 

a.  Integral  equation  development 

Figure  C-la  shows  the  side  view  of  a  material  plate  with 
constitutive  parameters  (u,e),  and  of  thickness  t,  immersed  in  free 
space,  i.e.,  a  homogeneous  medium  with  constitutive  parameters  (uo,^o). 
The  plate  is  excited  by  the  known  impressed  currents  (Jj.Mf)  which  in 
free  space  radiate  the  incident  fields  (Ej.Hj),  and  the  unknown  fields 
(£,HJ  in  the  presence  of  the  plate.  Let  S  denote  the  surface  enclosing 
the  plate.  Using  the  surface  equivalence  principle.  Figure  t^-la  can  be 
decomposed  into  the  two  cases  shown  in  Figures  C-lb,c  [10].  In  Case  A, 
the  constitutive  parameters  of  all  space  are  (uo»eo),  the  fields  are 
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zero  inside  the  plate,  the  unknown  fields  exist  outside  the  plate 

and  the  surface  currents  (£,M)  flow  on  the  surface  S.  The  currents 
U,M)  radiate  the  scattered  fields  (£S,HS)  in  free  space.  In  Case  B, 
the  constitutive  parameters  of  all  space  are  (u,e),  the  fields  are  zero 
outside  the  plate,  the  unknown  fields  (E^»H^)  exist  inside  the  plate,  and 
the  surface  currents  (J.1  ,M* )  =  (-J,-M)  flow  on  the  surface  S.  We  will 
consider  that  the  plates  are  sufficiently  thin  so  that  the  surface 
currents  on  the  sides  of  the  plate  will  be  ignored.  Note  that  use  of 
the  surface  equivalence  principle  allows  one  to  evaluate  the  fields 
either  inside  or  outside  a  material  body,  although  one  must  use  the 
appropriate  equivalent  currents  (l.e..  Case  A  outside  or  Case  B  inside) 
Further,  this  statement  is  not  changed  by  the  use  of  an  abc  such  as 
surface  impedance  relationships. 

Referring  to  Figure  C-l,  define: 

Sj  =  the  top  of  surface  S,  l.e.,  that  part  of  S  at  x  =  t/2 

Sp  =  the  bottom  of  surface  S,  l.e.,  that  part  of  S  at 
x  =  -t/2 

£(x,y,z)  =  the  electric  current  flowing  on  the  closed 

surface  S 

Jj(y.z)  =  the  top  electric  current,  l.e.,  that  part  of 
£  on  Sj 

Jp(y,z)  =  the  bottom  electric  current,  l.e.,  that  part 
of  on  Sfl 

J5(y,z)  =  the  sum  electric  current,  l.e.,  a  single  surface 

current  equal  to  the  vector  sum  Jy  +  and 
arbitrarily  considered  t:o  be  on  5y 


v  '_  -.  -V-^ 


Jp(y,z)  =  the  difference  electric  current,  I.e.,  a  single 
surface  current  equal  to  the  vector  difference 
Jj  -  Jp  and  arbitrarily  considered  to  be  on  Sy 

M,  Wf,  Mg,  M$,  and  Mp  have  identical  meanings,  except  that 

they  refer  to  magnetic  currents. 


Now  consider  the  problem  of  Figure  C-lb.  Since  a  null  field  exists 
Interior  to  S,  there  will  be  zero  reaction  between  the  currents  (J^M) 
and  (Jj.Mf),  and  a  test  source  interior  to  S.  If  the  test  source  is 
electric  with  current  density  J$»  this  leads  to  the  electric  reaction 


integral  equation  (ERIE) 


-  <Jt,(J,M)>  =  VJ 


where 


VJ  =  <Jt,(Ji,Mi)> 


By  definition  [11],  the  reaction  between  sources  a  and  b  is 


(C.l) 


(C.2) 


<a,b>  =  jj  (Ea  •  Jk  -  Ha  •  Mb)  ds  =  <b,a>  (C.3) 

b  ”  ” 

where  the  integral  is  over  the  region  occupied  by  the  source  b. 

Figure  C-2  shows  the  side  view  of  the  plate  with  an  electric  test 

source  £t  the  Interior.  considered  to  be  a  surface  current  in 

UJ 

the  yz  plane,  and  its  fields  in  free  space  will  be  denoted  (<j-,  j^) . 

The  plate  currents  (<^,M)  are  explicitly  shown  as  ( Jj ,Mj )  on  Sy  and 
(J_BtM]B)  on  Sr.  Using  the  linearity  of  the  reaction  operator.  Equation 
(C.l)  becomes 


-  <it.(iT,MT)>  -  <Jt,(J£,M3)>  =  VJ 


or 


(C.4) 


-  J  j  (EJ  •  Jt  -  HJ  •  My)  ds  -  J  J  (E°  •  Jr  -  H°  •  MB)ds  =  , 

ST  “t  "  “t  - 

f  C.  5) 

Since  the  tangential  electric  and  magnetic  fields  of  Jy  are  even  and  odd 
functions  of  x,  respectively 

£t(V  x  nT  =  it(SQ]  X  nT  (C.6a) 

Ht(Sj)  x  nj  =  x  nj  =  Jy/2  •  (C.6b) 

“  —  lim  t+o  — 

Using  the  symmetry  relations  of  Equations  (C.6),  we  move  the  bottom 
currents,  (£b»Mb)»  from  SB  to  Sy,  and  Equation  (C.5)  becomes 

J  J 

-  J  J  [E  •  (Jy  +  JB)  -  H  •  (Mr  -  MB)1  ds  =  VJ  ,  or  (C.7) 

Sy  -t 

-  <Jt.(J5,Mp)>  =  VJ.  (C.8) 

Note  that  Equation  (C.8)  applies  on  the  open  surface  Sy. 

Next  consider  the  case  where  the  test  source  in  Equation  (C.3)  is 
the  magnetic  current  My.  This  leads  to  the  magnetic  reaction  integral 
equation  ( MR  IE) 

-  <Mt,(J,M)>  =  VM  (C.9) 

where 
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VM  =  <Mt,(Ji,M1)>  . 


(C.10) 


Figure  C-2  shows  the  side  view  of  the  plate  with  a  magnetic  test  source 

Mt  in  the  interior.  Mt  is  considered  to  be  a  surface  current  in  the  yz 

_M  M 

plane,  and  its  fields  in  free  space  will  be  denoted  (£t»  _t)* 

Again,  using  the  linearity  of  the  reaction  operator.  Equation  (C.9) 
becomes 


-  <Mt,(Jj,Mj)>  -  <Mt,(JB,MB)>  »  VM  , 
or 


(C.ll) 


MM  MM 

-  J  J  (E  •  JT  -  H  •  MT)ds  -  j  j  (E  •  Jp  -  H 
S  j  ^  ^  Sg  ^  ^ 


•  Mg)ds  =  VM  . 

(C. 12 ) 


Since  the  tangential  electric  and  magnetic  fields  of  M-t  are  0,1(1  an(1  even 
functions  of  x,  respectively 


£t(V  x  nT  =  “  £t(V  x  nT  =  £t/Z  •  (C.  13a ) 

lim  t+o 

H^(Sj)  x  n^  =  H^(Sg)  x  nT  •  (C.13b) 

Using  the  symmetry  relations  of  Equations  (C.13),  we  again  move  (Jb.^b) 
from  Sr  to  Sj,  and  Equation  (C.12)  becomes 

-JJ  l£t  ‘  -  V  “  •  (£r  +  V1  ds  =  v”  ■ 


(C.14) 


-  <  Mt,  >  »  VM  ,  (C.15) 

valid  on  the  open  surface  Sj. 

Equations  (C.8)  and  (C.15)  constitute  two  of  the  four  equations 
needed  to  solve  for  the  four  unknowns,  an^  In 

obtaining  these  equations,  we  did  not  require  that  the  plate  thickness 
go  to  zero.  The  thickness  does  have  to  be  sufficiently  small  that  the 
side  currents  can  be  ignored.  To  obtain  the  remaining  two  equations, 
the  zero  reaction  concept  with  test  sources  outside  S  is  applied  to  the 
equivalent  problem  of  Figure  C-lc.  An  arbitrary  test  source  (^t»^t) 
exterior  to  S  in  Figure  C-lc  has  zero  reaction  with  the  currents 
(J',M')  =  (-J,-M),  or 

<(it,”t),(i,M)>  =0-  (C. 16) 

Explicitly  showing  the  integrations  over  the  top  and  bottom  surfaces. 
Equation  (C.16)  becomes 

ii  (Et  *  JT  -  Ht  *  +  JJ(Et  *  Jr  "  Ht  *  MB)ds  =  0  • 

St  -  “  ~  -  SB  " 

(C.17) 

where  (£t»Ht)  are  the  fields  of  (^ti^t)  in  the  homogeneous  medium  (u,e). 
In  practice,  two  independent  equations  can  be  obtained  from  Equation 
(C.16)  by  first  placing  (£t»!!t)  i°  the  plane  x  =  (t  +  At)/2  and  then  in 
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the  plane  x  =  -(t+At)/2  where  At  Is  a  positive  value.  Thus,  Equations 
(C.8),  (C.15),  and  (C. 16 )  twice  constitute  four  coupled  integral 
equations  which  could,  in  principal,  be  solved  for  and 

In  order  to  properly  account  for  plate  thickness,  one  should  convert  the 
sum  and  difference  currents  to  top  and  bottom  currents  prior  to 
computing  the  scattered  fields  (Interior  or  exterior  to  the  plate). 
However,  for  thin  plates  and  exterior  field  points  whose  x  distance  from 
the  plate  Is  large  In  comparison  to  plate  thickness,  one  may  compute  the 
scattered  fields  from  (£5,^3)  only. 

2.  Impedance  Boundary  Condition 

In  this  section  we  will  obtain  the  surface  impedance  relationships 
for  a  penetrable  dielectric/ferrite  plane  sheet  or  slab.  We  will  also 
obtain  the  surface  impedance  of  an  impenetrable  perfect  conductor  backed 
on  one  or  both  sides  by  a  dielectric/ferrite  slab.  This  will  be  done  by 
considering  the  simple  one  dimensional  problem  of  a  plane  wave 
normally  incident  upon  a  dielectric/ferrite  slab.  Although  the  solution 
to  this  problem,  in  terms  of  plane  wave  reflections,  appears  In  many 
texts  on  electromagnetic  theory  [12],  we  will  Investigate  the  solution 
using  the  surface  equivalence  principle  and  the  reaction  integral 
equations  obtained  in  the  previous  section.  For  the  case  of  the  slab, 
Equation  (C.16)  has  simple  solutions  from  which  we  can  deduce 
relationships  between  the  equivalent  electric  and  magnetic  surface 


H 

N 
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currents,  l.e.,  the  slab  surface  impedance  relationships.  In  the  next 
section  we  will  employ  the  approximation  that  these  relationships  {which 
are  exact  for  the  slab)  apply  also  to  the  finite  or  three  dimensional 
material  plate.  This  will  eliminate  two  of  the  four  Integral  equations 
for  the  material  plate. 

Figure  C-3a  shows  a  plane  wave  of  the  form 


£l  =  2  Eo  e 


(C.lBa) 


ii  -  i  Eo  eY°x  \ 


( C . 18b) 


where  y<>  =  fuoeo  and  n0  »  Ajo  Incident  upon  a  slab  of  thickness  t 

V  e0 

and  with  parameters  (u,e).  The  Illuminated  surface  of  the  slab  at 
x  =  t/2  will  be  referred  t r  as  Sj,  the  top  surface.  The  surface  at 
x  =  -t/2  will  be  referred  to  as  Sb,  the  bottom  surface.  In  direct 
analogy  to  Figure  C-l,  the  equivalent  currents  will  be  constant  sheets 

A  Ak  A  A 

denoted  ( «MTy ^  on  ST  and  (JBz»"MBy^  on  SB*  Referring  to  Figure  C-lc 
and  Equation  (C.16),  there  will  be  zero  reaction  between  a  constant 
sheet  current  test  source,  located  outside  the  slab,  and  the  equivalent 
currents.  If  the  test  source  Is  given  by  the  electric  current 


it  *  -2z 


(C.19) 


and  Is  located  in  the  plane  x  *  xt  (Utl>t/2),  then  Its  fields  In  the 
medium  (w,e)  are 
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Figure  C-3. 


(a)  Side  view  of  material  slab  Immersed  In  free  space, 

(b)  Side  view  of  a  coated,  perfectly  conducting  sheet 
Immersed  In  free  space, 

(c)  Transmission  line  equivalent  of  (a),  and 

(d)  Transmission  line  equivalent  of  (b). 
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E  t  -  z  n  e 


-Ylx-Xtl 


Ht  *  -y  sgn  (x-xt)  e 


-Ylx-xtl 


(C.20a) 


(C.20b) 


where  n  *  JL  ,  y  *  jw  ,  and  sgn(x)  «  1  If  x  >  0  and  -1  If  x  <  0. 
e 

Inserting  Equations  (C.19)  and  (C.20)  Into  Equation  (C.16)  yields 


-Yt/?  yt/2  -yt/?  yt/y 

JTne  6  +  JBne  c  +  f^e  ,£  -  M0e  =  0  for  xt  <  -t/2  , 

(C.21a) 


,  Tt/2  ,  -yt/2  M  Yt/2  -Yt/2  .  „ 

Jjne  +  J0ne  -  Kj.e  +  M^e  =  0  for  x^  >  t/2  . 


Solving  Equations  (C.21)  for  Hr  and  Mb  yields 


Mr  -  zttjt  ♦  ztbjb  , 


MB  “  ZBTJT  +  ZBBJB 


(C.21b) 


(C.22a) 


( C . 22b ) 


where 


ZTT  *  ZBB  *  T,/tanhTt  » 


( C . 23a ) 


ZTB  *  ZBT  *  r'/s1nh^t 


(C . 23b) 


These  Impedances  are  equal  to  the  two  port  Impedance  parameters  of  the 
transmission  line  of  length  t,  characteristic  Impedance  n,  and 


n.lYlVt’ 


propagation  constant  y,  shown  In  Figure  C-3c.  They  are  comparable  to 
equations  by  Mltzner  for  an  absorbing  shell  [25]. 

If  the  slab  had  been  Impenetrable,  such  as  the  perfectly  conducting 
ground  plane  coated  on  one  or  both  sides  by  a  dielectric/ferrite 
slab  of  Figure  C-3b,  then  the  two  port  Impedances  would  be  those  of 
the  transmission  line  model  of  Figure  C-3d.  In  this  case  we  have 
ZTB  =  ZBT  *  and 

Z-jt  =  njtanhyjtj  ,  (C.24a) 

ZBB  "  •  (C.24b) 

3.  Surface  Impedance  Approximation 

In  the  last  section,  two  Impedance  boundary  relationships 
(Equations  (C.22))  for  a  one  dimensional  dielectric/ferrite  slab  were 
obtained.  In  this  section.  Equations  (C.22)  are  used  as  an 
approximation  for  Equation  (C.16).  For  a  one  dimensional  slab  at  normal 
incidence  these  relationships  are  exact.  For  a  three  dimensional 
material  plate  they  are  not  exact  since: 

(1)  the  Incident  wave  may  not  have  normal  Incidence, 

(2)  the  edges  and  corners  of  the  plate  will  modify  the 
relationships. 
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Despite  this,  we  will  use  the  one  dimensional  slab  surface  Impedance 
relationships  for  the  three  dimensional  material  plate.  This  Is  an 
application  of  what  Senior  [13]  termed  an  approximate  boundary  condition 
(abc).  As  pointed  out  by  Senior,  there  usually  Is  no  rigorous 
justification  for  the  use  of  an  abc.  The  justification  Is  In  their 
usefulness  In  predicting  results  obtained  by  measurements  or  by  more 
rigorous  techniques.  Of  course,  the  motivation  for  their  use,  as 
compared  to  more  rigorous  techniques.  Is  that  they  simplify  the 
analysis.  It  Is  possible  that  better  results  would  be  obtained  by  using 
a  surface  Impedance  which  Is  dependent  upon  the  angle  of  Incidence, 
however,  for  simplicity  this  was  not  done  here.  Writing  Equation  (C. 22 ) 
In  terms  of  sum  and  difference  currents  gives 


!!s  *  zssi  x  "t  +  zsd£p  x  "t 


(C.25a) 


!!d  a  zds£s  x  "T  +  ZD0i)  x  nT  . 


where 


ZSS  *  (ZTT  +  ZTB  "  ZBT  '  ZBB^/2  » 


( C . 25b) 


(C.26a) 


ZSD  =  (ZTT  “  ZTB  “  ZBT  +  ZBB^/2  » 


(C . 26b) 


Combining  Equations  (C.8),  (C.15),  and  (C.25)  yields 
"  =  -<it’(i>’ZDsi>x  nT  +  ZDD^)X  nT^>  *  yJ  *  (C.27a) 

~  <^V^)>  =  -^.(Jp.Zssi51  nT  +  ZSd£px  s  vM  •  (C.27b) 

Equations  (C.27)  represent  two  coupled  Integral  equations  on  the  open 
surface  $t  which  can  be  solved  for  the  electric  currents  J$  and  Jp. 
Having  solved  for  J$  and  Jp,  the  magnetic  currents  and  Mp  can  b® 

a 

found  from  Equation  (C.25).  We  reemphasize  that  the  ZDSl!sxnT  +  ZDDi)xnT 

A  A 

and  zss£sxnT  +  zSpfpxnT  terms  In  Equation  (C.27)  are  not  electric 
currents,  but  rather  are  the  magnetic  currents,  ^p  and  respectively. 

In  the  remaining  part  of  this  section.  Equation  (C.27)  Is  examined 
for  various  special  cases. 

CASE  1  -  PERFECTLY  CONDUCTING  PLATES 

If  the  plate  Is  a  perfect  conductor,  then  Z$s  =  Zqo  =  Zds  =  Z$d  *  0 

and  Equation  (C.27)  reduces  to 

-<Jt.(Js.°)>  ■  v°  ,  (C.28a) 

-«*.($>.<>)>  *  V"  .  (C .  28b) 

Equation  (C.28)  was  previously  presented  and  solved  by  the 
authors  [14]. 


CASE  2  -  UNIFORM  MATERIAL  PLATE 

If  the  plate  is  of  uniform  u  and  e,  then  from  Equations  (C.23)  and 
(C.26),  Z$$  *  Zoo  *  0  and 

ZS0  =  ZTT  ‘  ZTB  =  n(coshYt-l)/sinhYt  ,  (C.29a) 

Z0S  =  ZTT  +  ZTB  =  n(coshYt+l)/sinhYt  .  (C.29b) 

In  this  case.  Equation  (C.27)  reduces  to 

“  =  "<Jt,^!ls*zDs£sx  nT^  =  * 

( C . 30a ) 

"  <Mt»^i)*MS^  =  "<Mt*^JD,ZSDi)x  nT^  *  V 

”  (C.30b) 

An  Important  point  is  that  the  equations  uncouple  for  this  case. 

However,  note  that  both  equations  must  be  solved  to  find  (£5,^5)  which 

radiate  the  far-zone  scattered  field.  In  general,  for  any  symmetric 

slab  (i.e.,  Zjt-Zbb  and  Zjb'Zbt)  t,,e  nations  uncouple, 

CASE  3  -  VERY  THIN  UNIFORM  MATERIAL  PLATE 

Now  consider  Equation  (C.30a)  in  the  limit  as  the  thickness  Yt  goes 
to  zero.  By  reciprocity.  Equation  (C.30a)  can  be  written  as 


where  £(J$)  and  £{Mp)  are  the  free-space  electric  fields  of  J$  and  Mp  = 

* 

Zqs  J3  xnTt  respectively,  and  the  Integral  Is  over  the  surface  of  the 
arbitrary  electric  test  source  J$.  Note  that  Is  located  in  the 

Interior  of  the  closed  surface  S  and  is  directed  tangential  to  the  broad 
surfaces  of  S. 

Since  Equation  (C.31)  must  apply  for  an  arbitrary  Jj»  it  is  clear 

that 

itan'^S1  *  £tan(V  '  '  Li  tan  ,nter,or  t0  s-  (C.32) 

As  yt  0,  the  tangential  electric  field  of  Mp  will  be  essentially  {see 
Equation  (C.13a)) 

£tan{V  =  ‘  nT  x  V2  ‘  (C.33) 

Also  from  Equation  (C.29b). 

lfm  Z  *  2n  .  ( C . 34 ) 

Yt+o  DS  Yt 

Combining  Equations  (C.25b),  (C.32),  (C.33),  and  (C.34)  gives 

"Itan^1  *  ^  =  II  tan  •  (C.35) 

Equation  (C.35)  gives  the  current  on  a  surface  of  sheet  Impedance  n/yt. 
It  Is  equivalent  to  that  used  by  Harrington  to  study  thin  dielectric 
shells  [15]  and  by  Senior  to  study  resistive  strips  [16].  Note  that 


when  treating  Impedance  sheets,  one  simply  adds  the  term  n  Jc  to  the 

yT 

equation  for  perfectly  conducting  surfaces.  In  a  MM  solution  of 
Equation  (C.35)  this  would  result  In  an  Impedance  matrix  of  the  form 
[Z]  +  [aZ]  where  [Z]  Is  the  Impedance  matrix  for  the  perfectly 
conducting  surface,  and  a  typical  term  of  [aZ]  would  be 

*  h  U  in  •  d»  '  ,C-36> 

Here  ^  represents  the  mth  current  expansion  mode  and  the  n^  test 
mode  or  weighting  function.  Equation  (C.36)  can  be  evaluated  in  closed 
form  for  any  of  the  usual  choices  for  the  modes.  This  permits  the  thin 
symmetric  slab  to  be  treated  as  a  very  simple  modification  of  the 
perfect  conductor.  Note  that  AZmn  a  0  unless  modes  and  overlap. 
Also,  for  extremely  small  yt,  AZ^o  will  be  large,  and  this  will  result 
In  small  £5. 

CASE  4  -  COATED  CONDUCTING  PLATE 

If  the  plate  Is  a  perfect  conductor  coated  on  one  or  both  sides  by 
a  dielectric/ferrite,  then  from  Equations  (C . 25 )  and  (C.26) 

ZSS  a  ZDD  *  (nitanhYjtj  -  n2tanhy2t2)/2  .  (C.37a) 


ZSD  =  ^DS  *  (nj^an^Yjt^  +  n2tanhy2t2)/2  . 


( C . 37b ) 


and  Equations  (C.27)  reduce  to  the  coupled  equations 


~  *  -<£.t*^^5*^SD2sx  nT  +  zss£px  nT^>  =  ^  *  (C.38a) 

"  <Mt.CJp.M5)>  55  "^t’C £p*ZSs£sX  nT  +  ZSDi)X  nT^>  *  ^  •  (C. 38b) 


4.  Moment  Method  Solution 

The  current  distributions  on  the  material  plates  can  be  found  by 
solving  the  Integral  equations  of  the  previous  section  by  the  numerical 
technique  known  as  the  method  of  moments  (MM)  [18].  Here  we  will 
present  the  MM  solution  to  Equations  (C.27)  which  represent  the  most 
general  case,  rather  than  one  of  the  special  cases  which  follow.  The 
solution  Is  similar  to  that  for  perfectly  conducting  surfaces  [1,2,7], 
and  thus  our  description  will  be  brief,  emphasizing  only  those  parts 
which  are  new. 

The  first  step  In  the  MM  solution  Is  to  expand  the  unknown 
currents,  and  £p  In  a  finite  set  of  N  basis  functions: 


Then  using  the  abc's  of  Equations  (C.25),  and  are  the  dependent 
unknowns 

3s  -  l  (ZSS  !Sn  +  ZSD  !Dn^  fjix  nT  »  (C.40a) 

n=l 

N 

4s!  (ZDS  !Sn  +  ZD0  W  fnx  nT  •  (C.40b) 

n=l 

Choosing  our  test  modes  Identical  to  the  expansion  modes  (except  that 
they  are  In  the  yz  plane),  Equations  (C.27)  reduce  to  the  set  of  2N 
simultaneous  linear  equations  written  in  block  matrix  form  as 


'lz!E!  1  £°s  Lzf"l’  _  _  Ido  _[£]_  I  I  isj  .  I  voj  , 

|  zss  [z*T]  (  [z®]  *  z$D  [znm]  j  j  To  |  j  *»n 

(C.41) 


where 


Z™  *  -  II  ^  •  (C.42a) 

n 

Z™  ’  -  IJ  !!iJ  •  (Ip  *  "t)  ds  ■  (C.42M 

C  -  -  JJ  iS  •  In  ds  •  (C.42C) 
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,MM 

4mn 


(C.42d) 


■  ij  •  (in *  "t>  ds  • 

Vm  =  JIJ  (if  •  i  •  s!>  .  <c.42e) 

V"  =  UJ  (if  •  ill  -  Hi  •  $  dv  •  (C.42f ) 

The  integrations  in  Equations  (C.42a-C.42d)  are  over  the  surface  of  the 

j  j  MM 

nth  expansion  mode  located  on  Sj.  (E^,  and  (J^,,  £,„)  are  the  free 
space  fields  of  the  mth  electric  and  magnetic  test  modes,  respectively. 
Since  the  expansion  and  test  modes  used  in  this  solution  are  identical 
(except  for  a  shift  of  t/2),  the  method  is  essentially  a  Galerkin 

solution,  and  each  of  the  matrices  [ZEE]  and  [ZMM]  are  symmetric.  Also 
[ZEM]  >  .  [ZME]T. 

A  typical  rectangular  plate  is  shown  in  Figure  C-4.  Two  orthgonal 
and  overlapping  sets  of  surface-patch  modes  cover  the  surface  of  the 
plate  forming  a  two-dimensional  vector  surface  density  which  could  be 
used  to  represent  £5,  £p,  or  Mp.  The  normal  components  of  £5  an<l 
must  vanish  at  a  plate  edge,  however,  the  normal  components  of  £p  and 
do  not.  Thus,  we  include  a  series  of  edge  modes,  consisting  of  a  single 
monopole  with  terminals  at  the  plate  edge.  This  allows  for  a  non-zero 
normal  component  of  current  at  the  plate  edge.  The  plate  in  Figure  C-4 
is  made  up  of  ten  monopole  edge  modes  and  seven  dipole  interior  modes. 
The  monopole  edge  modes  will  be  identically  zero  for  £5  and  ^D»  but  not 

for  £p  and 


i, 
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5.  Numerical  Results 


In  this  section  we  will  present  numerical  results  for  electro¬ 
magnetic  scattering  by  material  plates.  Our  results,  computed  using 
Equations  (C.41)  and  (C.42)  will  be  compared  with  measurements. 

Figure  C-5  shows  the  radar  cross  section  (RCS)  versus  angle  of 
Incidence  for  a  0.813A  square  by  0.051A  thick  penetrable  dielectric 
plate  with  relative  dielectric  constant,  er  *  e/e0  =  4.  Both  TE  and  TM 
polarizations  are  shown  and  compared  with  measurements  by  Thai  and 
Finger  [18].  The  agreement  Is  seen  to  be  good  near  broadside,  however. 
It  worsens  as  one  approaches  grazing.  This  Is  to  be  expected,  since  the 
surface  Impedance  relationships  were  evaluated  for  a  normally  Incident 
wave.  Also  shown  In  Figure  C-5  Is  the  RCS  computed  by  the  Impedance 
sheet  approximation  of  Equation  (C.35).  The  sheet  impedance 
approximation  produces  a  less  accurate  result  since  It  requires  |yt|«l, 
while  In  this  case  I rtl =0.64.  Our  experience  Indicates  that  lytl  should 
be  less  than  about  0.2  to  use  the  Impedance  sheet  approximation. 

The  Insert  in  Figure  C-6  shows  a  1.5A  square  perfectly  conducting 
plate  coated  on  one  side  by  a  narrow  band  dielectric/ferrite  absorber 
with  er=7.8-jl.6,  yr=1.5-j0.7,  and  thickness  0.065A.  This  Is  an  example 
of  an  Impenetrable  plate  with  asymmetric  surface  Impedance 
relationships.  Figure  C-6  shows  the  RCS  of  this  plate  evaluated  with 
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.  (a)  TE  radar  cross-section  of  a  square  dielectric  plate 
(b)  TM  radar  cross-section  of  a  square  dielectric  plate 
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the  surface  Impedance  relationships  of  Equation  (C.37),  and  compared 
with  measurements.  Note  that  the  theory  correctly  predicts  about  a  20 
dB  drop  In  RCS  when  the  Incident  field  is  broadside  to  the  absorber, 
rather  than  the  perfect  conductor. 

Figure  C-7  shows  the  RCS  at  broadside  for  a  fiberglass  strip  dipole 
of  width  O.lx,  thickness  0.032a,  and  er=4.8  versus  the  dipole  length  L. 
Our  compuratlon  Is  slightly  below  the  previous  computations  (based  upon 
the  more  exact  equivalent  volume  current  representation)  and 
measurements  of  Richmond  [19].  Since  the  surface  impedance 
approximation  Is  felt  to  be  most  accurate  away  from  corners  and  edges, 
the  degree  of  error  In  Figure  9  Is  probably  less  than  one  would  expect, 
considering  that  It  Is  being  applied  to  an  extremely  small  plate. 

Figure  C-8  shows  the  RCS  at  broadside  for  an  aluminum  strip  dipole 
of  width  0.0056 A,  conductivity  35  megamho/meter  (l.e.,  skin  depth  6  = 
1.03  x  10-5)  versus  the  dipole  length  L.  Curves  are  shown  for 
thicknesses  of  t  a  6/1000,  6/100  and  6/10,  and  6.  The  agreement  between 
our  results  and  the  more  exact  results  of  Richmond  [19]  Is  seen  to  be 
good,  especially  for  the  larger  t  values. 
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Figure  C-8.  Radar  cross-section  of  rectangular  aluminum  strips. 


6.  Summary 


An  Integral  equation  formulation  for  the  electromagnetic  scattering 
from  a  thin  material  plate  has  been  presented.  The  material  plate  Is 
viewed  as  an  open  surface.  In  order  that  the  solution  appear  as  an 
extension  or  modification  of  that  for  perfectly  conducting  plates. 
Solutions  for  the  sum  and  difference  electric  and  magnetic  equivalent 
surface  currents  are  obtained  In  the  form  of  four  coupled  integral 
equations.  By  using  the  abc  or  surface  Impedance  relationships  two  of 
the  four  Integral  equations  are  eliminated.  The  surface  impedance 
relationships  are  shown  to  be  the  two-port  Impedance  parameters  of  a 
simple  transmission  line  model.  By  Including  the  off-diagonal  terms  In 
the  two-port  Impedance  matrix,  penetrable  as  well  as  impenetrable  plates 
can  be  treated.  Several  special  cases  of  the  two  coupled  integral 
equations  are  Investigated,  and  It  Is  shown  that  for  electrically  very 
thin  plates,  the  results  simplify  to  a  single  equation  equivalent  to  the 
Impedance  sheet  approximations  used  by  Harrington  [15]  and  Senior  [16]. 
The  equations  were  solved  by  the  method  of  moments,  and  accurate  results 
were  obtained.  Future  work  In  this  area  should  include  solutions  for 
non-planar  geometries,  and  the  Intersection  of  a  perfectly-conducting 
and  a  material  plate. 
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D.  Scattering  by  Absorber  Geometries 

Researchers:  J.H.  Richmond,  Professor  (Phone:  (614)  422-7601) 

The  increased  use  of  materials  including  absorbers  on  aerospace 
vehicles  leads  to  requirements  for  more  efficient  analytic  treatment  of 
such  configurations.  One  such  analysis  was  developed  under  another 
project  [1].  The  expansion  of  this  technique  requires  a  substantial 
basic  research  program  to  support  the  application  of  these  materials. 
Thus  this  research  is  to  be  expanded  on  the  Joint  Services  Program.  The 
original  results  presented  in  [1]  are  outlined  in  the  following  pages. 

Accompl i shments 

1.  Introduction 

An  efficient  moment-method  solution  (with  only  three  equations  and 
three  unknowns)  was  developed  for  plane-wave  scattering  by  a  thin  lossy 
dielectric  strip  with  infinite  length.  We  expanded  the  field  in  the 
dielectric  region  as  the  sum  of  three  plane  waves:  the  forced  wave  and 
two  traveling  surface  waves.  Galerkin's  method  was  applied  to  determine 
the  amplitudes  of  these  waves  and  finally,  the  far-zone  scattered  field 
was  determined  by  considering  the  electric  polarization  currents 
radiating  in  free  space.  This  work  is  reported  in  reference  [1], 


111 


2.  Scattering  by  Ferrite  and  Ferrite-Coated  Cylinders  and  Strips 

Background 

Although  Richmond  [2,3]  published  the  first  moment-method 
solutions  for  scattering  by  dielectric  cylinders  with  arbitrary 
cross-sectional  shape  In  1965,  the  technique  has  not  yet  been  extended 
to  ferrite  cylinders.  A  dielectric  cylinder  is  modeled  with  electric 
polarization  currents,  while  a  ferrite  cylinder  requires  magnetic  (as 
well  as  electric)  polarization  currents. 

In  1973,  Richmond  [4]  developed  a  user-oriented  moment  method  for 
scattering  by  a  perfectly  conducting  polygon  cylinder.  Later,  Wang  [5] 
extended  the  theory  and  computer  program  to  Include  a  thin  dielectric 
coating  on  the  polygon  cylinder.  While  promising,  Wang's  results  might 
be  Improved  by  Including  surface  waves  In  the  analysis. 

In  the  area  of  RCS  control ,  Professor  Leon  Peters  has  suggested 
that  a  thin  ferrite  coating  may  offer  unique  properties  that  cannot  be 
duplicated  with  a  thin  lossy  dielectric  coating.  This  idea  has  been 
reinforced  by  the  numerical  results  of  Hill  [6]  for  the  propagation 
constant  of  a  surface  wave  on  a  thin  lossy  ferrite  coating. 

To  take  advantage  of  these  unique  properties,  we  require  an 
analysis  of  scattering  by  a  thin  lossy  ferrite  slab,  or  a  conducting 
strip  with  a  thin  ferrite  coating.  An  analysis  for  the  lossless  case 
has  been  developed  recently  by  Burnside  and  Burgener  [7]  with  the  GTD 
format.  Their  results  agree  closely  with  moment  method  calculations 


except  for  those  situations  where  surface  wave  effects  are  dominant. 

The  surface  waves  are  not  Included  In  their  analysis. 

As  discussed  In  an  earlier  section,  Pathak  and  Rojas-Teran  have 
applied  the  Wiener-Hopf  technique  to  the  thin  dielectric  half-plane. 

They  have  obtained  excellent  results,  except  for  difficulties  In  the 
vicinity  of  grazing  Incidence. 

Richmond  [1]  has  developed  an  efficient  moment-method  solution 
(with  only  three  equations  and  three  unknowns)  for  plane-wave  scattering 
by  a  thin  lossy  dielectric  strip  (see  Figure  D-l)  with  the  incident 
electric  field  vector  parallel  with  the  edges  of  the  strip  as  In  Figure 
D-2.  The  electric  field  intensity  in  the  dielectric  region  was 
expressed  as  the  sum  of  three  traveling  waves  as  follows: 

fjx  fx  -fx 

Ez(x,y)  *  Cje  cosh(g^y)  +  C2e  cosh(gy)  +  C3e  cosh(gy)  . 

The  first  term  represents  the  "forced  wave",  or  the  field  induced  In  an 
inflnltely-wide  dielectric  slab.  The  last  two  terms  represent  surface 
waves  traveling  In  opposite  directions  across  the  thin  dielectric  strip. 
Gal erk In's  method  is  employed  to  determine  the  constants  Ci,  C2  and  C3. 
Finally,  the  far-zone  scattered  field  is  calculated  by  considering  the 
polarization  currents  radiating  In  free  space.  Figure  D-3  illustrates 
the  electric  field  distribution  Induced  in  a  thin  dielectric  strip  by  a 
plane  wave  with  grazing  Incidence.  The  results  obtained  with  the  new 
method  (plane-wave  expansion)  show  excellent  agreement  with  those  from 
the  old  method  (pulse  bases)  which  involved  a  system  of  200  simultaneous 
linear  equations. 
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Future  Research 


We  propose  to  Investigate  scattering  by  thin  dielectric  strips  with 
the  aid  of  moment  methods  and  plane  wave  basis  functions.  Having 
already  demonstrated  the  technique  with  the  electric  field  parallel  to 
the  edges  of  the  strip,  we  propose  to  extend  It  to  the  following 
situations: 

a.  thin  dielectric  strip  with  the  other  polarization 

b.  thin  ferrite  strip 

c.  perfectly  conducting  strip  with  thin  ferrite  coating 

d.  perfectly  conducting  polygon  cylinder  with  thin  ferrite 
coatl ng. 

In  each  phase  of  this  Investigation,  the  results  calculated  by  the  new 
method  will  have  to  be  tested  against  experimental  measurements  or 
Independent  calculations.  Since  measurements  or  Independent  data  are 
not  available  In  most  cases,  alternative  formulations  and  programs  will 
have  to  be  developed  for  test  purposes.  This  In  Itself  will  represent  a 
substantial  effort  and  a  very  useful  product. 
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Introduction 

The  optimum  exploitation  of  the  Interaction  of  electromagnetic 
waves  and  material  objects  can  only  be  achieved  with  a  clear 
understanding  of  how  the  physical  properties  of  the  object  (size,  shape 
and  composition)  In  situ  are  related  to  scattered  or  radiated  fields. 
For  whatever  purpose,  e.g.,  target  Identification,  radar  cross  section 
control,  specific  radiation  waveforms,  etc.,  a  basic  understanding  of 
the  electromagnetic  Interactions  Is  vital.  Our  research  has  shown  that 
a  time  domain  approach  provides  Ideal  Insight  Into  how  such  an 
understanding  can  be  realized.  Increasingly,  solutions  to  radiation  or 
scattering  problems  are  sought  over  broad  bandwldths.  The  canonical 
response  waveforms  (response  waveforms  to  Impulse,  step  or  ramp 
excitation)  are  Ideally  suited  for  this  purpose  because  the  response  to 
an  arbitrary  excitation  Is  simply  obtained  using  convolution. 


Much  of  our  recent  research  has  stressed  the  complex  natural 
resonances  of  a  scatterer  or  radiator,  that  is,  those  complex 
frequencies  which  ellicit  singularities  in  the  transfer  function 
relating  excitation  and  response  in  a  linear  system.  For  distributed 
parameter  systems  such  as  electromagnetic  scatterers  and  radiators,  the 
response  to  an  aperiodic  excitation  consists  of,  first,  a  forced 
response  as  the  wavefront  moves  across  the  body  and  then  a  free  or 
natural  response  as  the  wavefront  moves  beyond  the  body.  For  fixed 
time-independent  residues,  only  the  free  response  can  be  described  in 
terms  of  the  natural  oscillation  frequencies  of  the  object. 

The  general  goals  of  our  research  include: 

1)  Prediction  of  scattered  and  radiated  waveforms  from  objects  of 
Increasing  complexity  for  point  or  plane  wave  electromagnetic 
sources  with  arbitrary  waveforms. 

2)  Investigation  of  target-dependent  excitation  waveforms  and/or 
signal  processing  algorithms  to  identify  and  optimize 
electromagnetic  response  from  specific  targets. 

3)  Determination  of  the  effects  of  system  bandwidth  and  noise 

on  recommended  techniques  for  detection  and  identification  of 
radar  targets. 


A  very  real  measure  of  the  significance  of  the  basic  research 
performed  on  a  program  Is  the  transfer  of  the  developed  techniques  and 
methods  to  other  less  basic  programs.  The  extraction  of  the  complex 
natural  resonances  of  a  target  from  measured  mul tl frequency  scattering 
data  using  rational  function  approxlmants  and  the  recognition  of  a 
target  using  such  resonances  In  prediction-correlation  processing  (both 
concepts  developed  on  this  program)  are  now  being  extensively  tested  on 
another  program/  This  other  program  has  resulted  in  a  unique  new 
broadband  compact  range  reflectivity  measurement  facility  capable  of 
rapidly  measuring  the  complex  polarization  scattering  matrix  of  model 
targets  over  a  12:1  band  (present)  with  planned  extensions  to  30:1  band 
(future)  and  possibly  a  100:1  band  (ultimately).  This  other  program 
(ESI  Project  714190)  has  as  one  major  goal  the  establishment  of 
prediction-correlation  target  recognition  statistics  using  complex 
natural  resonances  and  a  realistic  data  base.  Accordingly,  no  future 
research  on  routine  extraction  of  complex  natural  resonances  from 
measured  scattering  data  nor  on  application  of  prediction-correlation 
processing  for  target  recognition  Is  anticipated  on  this  program.  Our 
research  Is  summarized  in  the  following  paragraphs.  Written 
publications  and  oral  presentations  are  listed  at  the  end  of  this 
section.  We  stress  here  the  fact  that  we  have  been  extremely  active  In 
disseminating  our  results  to  the  research  community. 

*Project  714190,  Contract  #  N00014-82-K-0037  between  The  Ohio  State 
University  Research  Foundation  and  the  Department  of  the  Navy,  Office 
of  Naval  Research,  Arlington,  Virginia. 


Accompl 1 shments 


1.  Complex  Natural  Resonances  and  Geometrical  Procedures 

The  K-pulse  provides  a  means  of  relating  surface  waves  on  a 
structure  to  the  complex  natural  resonances  of  the  structure  [1].  For 
solid  structures  such  as  spheres,  prolate  spheroids,  circular  and 
elliptical  cylinders,  etc.,  the  method  has  been  found  to  be  limited  only 
by  the  geometrical  theory  of  diffraction  (GTD)  estimate  of  the 
attenuation  and  phase  shift  of  surface  rays  on  the  structure  [1,2]. 

That  Is,  the  first  one  or  two  lowest  frequency  (natural  oscillations) 
natural  resonances  predicted  are  Incorrect  because  the  estimates  are 
asymptotic.  Higher  frequency  "pole"  estimates  are  excellent  for  spheres 
and  circular  cylinders  and  good  to  adequate  for  other  classical 
geometries.  Therefore,  for  sufficiently  simple  solid  target  geometries 
a  method  for  obtaining  estimates  of  the  complex  natural  resonances  has 
been  delineated.  The  method  Itself  (circumferential  paths)  also 
refutes  the  Idea  that  the  complex  natural  resonances  are  simply  a 
mathematical  tool  and  unrelated  to  the  target's  physical  properties. 

For  wire  structures,  results  were  somewhat  less  satisfactory  except 
for  cases  where  the  wire  was  very  thin  (length  to  diameter  ratios  (L/D) 
greater  than  1000 ) [23.  It  has  been  found  that  for  "fatter"  wires. 
Improved  estimates  of  the  attenuation  and  excess  phase  shift  of  current 
on  the  wire  are  now  available  [3].  For  a  straight  wire,  these  formulas 
yield  pole  locations  with  Improved  estimates  of  the  damping  but  somewhat 


poorer  estimates  of  the  oscillation  frequencies.  By  properly  combining 
the  new  and  old  analytical  estimates  of  attenuation  and  phase  shift, 
good  estimates  of  the  pole  locations  can  now  be  obtained  for  both  "fat" 
and  "thin"  wires.  These  results  are  presently  being  used  to  modify  the 
analytical  estimates  of  the  pole  locations  of  bent  wires,  curved  wire 
segments,  loops  and  other  wire  geometries  [2].  The  delayed  technical 
report  [4]  on  these  procedures  is  now  being  modified  and  the  report 
should  be  available  by  early  1984. 

Related  questions  concern  the  frequency-dependent  current 
reflection,  transmission  and  radiation  characteristics  at  wire  junctions 
and  the  radiative  coupling  at  end  points  of  the  wires.  The  improved 
estimate  of  current  attenuation  and  excess  phase  shift  will  permit  a 
brief  reexi mi  nation  of  these  parameters  but  it  is  not  intended  to 
further  delay  the  technical  report  on  these  topics. 

2.  Complex  Natural  Resonances  Via  Rational  Function  Approximants 

A  definitive  report  on  the  use  of  rational  function  approximants  to 
extract  complex  natural  resonances  of  a  target  from  measured  multiple 
frequency  scattering  data  Is  In  process*[5].  It  has  been  found  that 
a  combination  of  filtering,  windowing,  averaging  and  judicious  selection 
of  test  frequencies  (resonance  region)  permits  a  consistent  set  of 

*The  processing  methods  developed  were  completed  on  this  contract. 
Application  of  the  methods  and  techniques  to  actual  measured  scattering 
data  was  done  on  the  related  contract  (ESL  Project  714190)  mentioned 
earlier.  Target  recognition  algorithms  are  being  studied  on  the  latter 
contract. 
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complex  natural  resonances  (poles)  to  be  extracted  from  measured 
scattering  data.  There  are,  however,  several  limitations.  These 
limitations  and  some  Initial  conclusions  are  listed  below: 

a)  The  measured  data  must  have  an  Integrated  signal  to  noise 
ratio  of  at  least  15.0  dB  before  consistent  poles  can  be 
extracted.  This  assumes  that  the  noise  Is  additive. 
Multiplicative  noise  remains  a  moot  question. 

b)  The  limitation  In  a)  suggests  that  In  most  situations  the 
scattering  data  for  estimating  poles  must  be  taken  In  a 
controlled  environment.  This  means  that  target 
recognition  algorithms  should  not  use  direct  extraction  of 
poles  from  full  scale  data.  Fortunately,  prediction- 
correlation  target  recognition  does  not  extract  poles 
directly  and  Is  therefore  unaffected  by  this  limitation. 

c)  The  extracted  pole  set  for  a  target,  while  clearly 
consistent,  does  not  display  precise  excitation 
Invariance.  Therefore,  as  a  practical  matter  a 
geometrically  complex  target  will,  for  target  recognition 
purposes,  have  to  be  considered  as  several  targets,  with  a 
set  of  averaged  poles  related  to  a  given  span  of  aspects 
and/or  polarizations. 


d)  The  report  In  process  also  uses  GTD  estimates  to  obtain 
the  creeping  wave  poles  for  a  thin  circular  disk  geometry. 

From  the  results,  the  K-pulse  and  response  waveform  for 
edge-on  Incidence  for  the  disk  are  presented. 

No  further  studies  of  processing  techniques  for  extracting  complex 
natural  resonances  from  measured  scattering  or  radiation  data  are 
anticipated  on  this  contract.  We  shall,  however,  continue  to  explore 
geometrical  and  asymptotic  procedures,  particularly  with  regard  to 
future  K-pulse  studies. 

3.  Transient  Current  Density  Waveforms  on  a  Conducting  Sphere 

A  technical  report  presenting  the  transient  current  waveforms 
Induced  on  a  conducting  sphere  when  Illuminated  by  a  plane 
electromagnetic  wave  with  Impulsive  time  dependence  has  been  published 
[6].  On  the  Illuminated  side  of  the  sphere,  the  Inadequacies  of  the 
physical  optics  approximation  In  both  the  E  and  H  planes  are  Illustrated 
and  a  simple  correction  form  Is  suggested.  On  the  shadowed  side  of  the 
sphere,  the  changing  form  of  the  current  density  waveforms  Is  discussed 
and  demonstrated.  This  report  demonstrates  that  first  order  corrections 
of  the  physical  optics  backscatter  approximation  are  possible  even  for 
an  object  which  does  not  have  a  depolarized  scattering  component.  This 
report,  with  some  additional  references  to  research  on  time  domain 
scattering  and  radiation  performed  at  the  ElectroSclence  Laboratory, 
will  also  be  published  In  a  book  of  NATO  Conference  Proceedings  [7]. 
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4.  Cavity  Structures 


Research  on  the  electromagnetic  scattering  from  finite  and  open 
circular  waveguide  structures  is  continuing.  The  ultimate  goal  is  a 
realistic  modeling  of  the  intake  and  exhaust  cavities  of  jet  engine 
configurations.  At  the  moment,  particular  emphasis  is  on  the  complex 
natural  resonances  of  such  structures  and  on  the  K-pulses  to  be 
constructed  from  such  complex  natural  resonances.  It  should  be 
emphasized  that  some  portions  of  the  research  on  cavity  scattering  have 
been  done  on  other  projects.* 

While  completely  consistent  results  for  all  aspects  of  the  problem 
have  not  yet  been  obtained,  some  preliminary  results  have  proved 
publishable.  A  paper  [8]  which  demonstrates  a  technique  for  using 
rational  function  approximations  to  combine  low  frequency  and  high 
frequency  scattering  estimates  has  been  accepted  for  publication  by  WAVE 
MOTION.  In  this  case,  the  low  frequency  results  come  from  the 
Wiener-Hopf  solution  for  an  open  circular  waveguide  [9]  and  the  high 
frequency  estimates  from  the  asymptotic  unified  theory  of  diffraction 
(UTD)  [10].  The  targets  in  this  case  are  finite  circular  waveguides 
with  various  Internal  (transverse  to  guide  axis)  loads.  The  rational 
function  approximations  permit  a  smooth  spanning  and  transition  of  the 
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University  Research  Foundation  and  the  Department  of  the  Navy,  Office 
of  Naval  Research,  Arlington,  Virginia  and  project  712661,  Contract 
#  F19628-80-C-0056  between  The  Ohio  State  University  Research  Foundation 
and  the  Department  of  the  Air  Force,  Electronic  Systems  Division,  Air 
Force  Systems  Command,  Hanscom  Air  Force  Base,  Massachusetts. 


spectral  range  between  the  lowest  cut  off  frequency  for  the  waveguide 
(TEn  mode)  and  a  higher  frequency  where  the  UTD  results  are  apparently 
valid.  Another  publication  [11]  stresses  a  rational  function 
approximant  for  extracting  complex  natural  resonances  from  measured  or 
calculated  multiple  frequency  scattering  data.  Included  here  also  are 
estimates  of  certain  of  the  complex  natural  resonances  of  finite  and 
open  circular  waveguides. 

The  problem  of  a  finite  circular  waveguide  loaded  with  a  shorted 
conical  structure  is  currently  being  studied  using  a  method  of  moments 
program  for  rotationally  symmetric  scatterers.  At  present,  a  portion  of 
the  scattered  field  spectrum  is  being  generated  and  the  poles  extracted 
using  a  rational  function  fit  approximation.  A  numerical  search 
procedure  will  ultimately  be  added.  The  new  compact  range  reflectivity 
facility  mentioned  in  the  introduction  will  materially  aid  in  providing 
verification  data  for  new  cavity  scattering  results  obtained  during  the 
next  interim. 

5.  K-pulse  Studies 

The  K-pulse  waveform  as  Intepreted  in  the  time  domain  is  a  unique 
waveform  of  minimal  duration  which,  as  input  to  a  specified  linear 
system  (scatterer  or  radiator)  produces  response  waveforms  of  finite 
duration  at  all  points  of  interest.  Derivations  of  the  input  waveform 
would  also  possess  this  property,  thus  uniqueness  requires  that  the 
highest  order  discontinuity  be  of  minimal  order  as  well.  It  must  be 
recognized  that  while  it  would  appear  that  any  linear  combination  of 


derivatives  of  finite  orders  would  appear  to  satisfy  the  minimal 
duration  test,  this  process  also  would  Introduce  new  zeros.  The  further 
condition  that  the  K-pulse  transform  zeros  match  only  the  system  poles 
must  be  Impressed.  The  K-pulse  for  a  distributed  parameter  system  can 
be  similarly  modified  to  be  the  K-pulse  for  this  distributed  parameter 
system  and  a  ’specified  lumped  content  system  without  further 
lengthening. 

In  a  recent  paper  [12],  Pearson  purports  to  show  that  while  the 
singularity  expansion  Is  an  Ill-suited  means  for  modeling  a  scatterer's 
response  to  a  time-limited  excitation  (since  it  cannot  be  made 
practically  useful  through  truncation),  one  is  mistaken  to  conclude  that 
it  is  Incorrect  or  Incomplete  through  this  fact.  We  disagree  with  this 
conclusion  primarily  because  of  his  simplified  definition  for  the 
Laplace  transform  of  a  K-pulse.  It  is  unlikely  that  this  disagreement 
will  be  resolved  in  the  near  future.  We  note,  however,  that  two  papers 
[13,14]  to  appear  In  the  near  future  both  demonstrate  that  there  is  a 
nonzero  region  In  the  space- time  domain  where  the  scattered  field  cannot 
be  expressed  in  terms  of  natural -mode  contributions  only  and  where  the 
conventional  SEM  representation  does  not  hold.  One  must  recognize 
however,  that  this  result  Is  for  a  linearly  polarized  electromagnetic 
pulse  of  finite  duration  incident  on  a  homogeneous  or  inhomogeneous 
lossy  dielectric  slab.  Thus,  here  the  target  is  neither  finite  nor 
perfectly  conducting. 
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The  question  raised  by  Kennaugh  [1]  regarding  the  completeness  of 
the  SEM  representation  has  no  bearing  on  the  future  usefulness  of  the 
K-pulse.  It  remains  a  unique  and  powerful  tool  in  identification  and 
imaging  and  also  appears  to  be  useful  in  antenna  studies.*  We  have  only 
recently  resumed  research  on  the  K-pulse  as  unpublished  notes  and 
memorandums  by  Professor  Kennaugh  became  available.  First  priority  is 
being  given  to  completion  of  the  report  "K-pulse  for  Transmission 
Lines". 

Recei:t  Publications  Under  Time  Domain  Studies 

D.L.  Moffatt,  C-Y  Lai  and  T.C.  Lee,  "Time  Domain  Electromagnetic 
Scattering  by  Open-Ended  Circular  Waveguide  and  Related  Structure", 
accepted  for  publication  in  WAVE  MOTION,  October,  1983. 

D.L.  Moffatt,  J.D.  Young,  A. A.  Ksienski,  E.  Lin  and  C.  Rhoads, 

"Transient  Response  Characteristics  in  Identification  and  Imaging"  in 
Microwave  Imaging  Techniques  edited  by  C.K.  Chan  and  T.L.  Lim,  a  series 
on  "Benchmark  Papers  in  Electrical  Engineering  and  Computer  Science", 
edited  by  John  B.  Thomas,  Hutchinson  Ross  Publishing  Company, 
Stroudsburg,  Pennsylvania,  1984. 


*See  JSEP  proposal.  Time  Domain  Studies,  for  1983-1984,  pp.  77-90. 


Electromagnetic  Imaging,  edited  by  W.M.  Boerner,  D.  Reidel  Publishing 
Company,  Dordrecht,  Holland,  1984. 

E.M.  Kennaugh  and  D.L.  Moffatt,  "Transient  Current  Density  Waveforms  on 
a  Spherical  Scatterer",  in  Proceedings  of  MATO  Advanced  Research 
Workshop  on  Inverse  Methods  in  Electromagnetic  Imaging,  edited  by  W.M. 
Boerner,  D.  Reidel  Publishing  Company,  Dordrecht,  Holland,  1984. 

T.W.  Johnson  and  D.L.  Moffatt,  "Electromagnetic  Scattering  by  an  Open 
Circular  Waveguide",  Radio  Science,  Vol .  17,  No.  6,  pp.  1547-1556, 
November-December,  1982. 

D. L.  Moffatt  and  C.M.  Rhoads,  "Radar  Identification  of  Naval  Vessels", 
IEEE  Transactions  on  Aerospace  and  Electronic  Systems,  Vol.  AES-18,  No. 
2,  March  1982. 

E. M.  Kennaugh  and  D.L.  Moffatt,  Comments  on  "Impulse  Response  of  a 
Conducting  Sphere  Based  on  Singularity  Expansion  Method",  Proceedings 
of  the  IEEE,  Vol.  70,  No.  3,  March  1982. 


L.C.  Chan,  D.L.  Moffatt,  and  L.  Peters,  Jr.,  "Subsurface  Radar  Target 
Imaging  Estimates",  IEEE  Trans,  on  Antennas  and  Propagation,  Vol.  AP-29, 
No.  2,  March  1981. 


E.M.  Kennaugh,  "Polarization  Dependence  of  RCS  -  a  Geometrical 
Interpretation",  IEEE  Trans,  on  Antennas  and  Propagation,  Vol .  AP-29, 

Mo.  2,  March  1981. 

L.C.  Chan,  D.L.  Moffatt,  and  L.  Peters,  Jr.,  "Improved  Performance  of  a 
Subsurface  Radar  Target  Identification  System  Through  Antenna  Design", 
IEEE  Trans,  on  Antennas  and  Propagation,  Vol.  AP-29,  No.  2,  March  1981. 
Invited. 

E.M.  Kennaugh,  "Opening  Remarks",  IEEE  Transactions  on  Antennas 
Propagation,  Vol.  AP-29,  No.  2,  March  1981. 

D. L.  Moffatt,  J.D.  Young,  A. A.  Kslenskl,  H.C.  Lin,  and  C.M.  Rhoads, 
"Transient  Response  Characteristics  in  Identification  and  Imaqlng",  IEEE 
Trans,  on  Antennas  and  Propagation,  Vol.  AP-29,  No.  2,  March  1981, 
Invited. 

E. M.  Kennaugh,  "The  K-pulse  Concept",  IEEE  Trans,  on  Antennas  and 
Propagation,  Vol.  AP-29,  No.  2,  March  1981. 

D.L.  Moffatt,  "Ramp  Response  Radar  Imagery  Spectral  Content",  IEEE 
Trans,  on  Antennas  and  Propagation,  Vol.  AP-29,  No.  2,  March  1981. 

D.B.  Hodge,  "Scattering  by  Circular  Metallic  Disc",  IEEE  Transactions  on 
Antennas  and  Propagation,  Vol.  AP-28,  No.  5,  September  1980. 


L.C.  Chan,  D.L.  Moffatt,  and  L.  Peters,  Jr.,  “Estimation  of  the  Complex 
Natural  Resonances  of  a  Class  of  Subsurface  Targets",  in  Acoustical , 
Electromagnetic  and  Elastic  Wave  Scattering  -  Focus  on  the  T-matrix 
Approach,  edited  by  V.K.  and  V.V.  Varadan,  Pergamon  Press,  1980. 

L.C.  Chan,  D.L.  Moffatt,  and  L.  Peters,  Jr.,  "A  Characterization  of 
Subsurface  Radar  Targets",  Proc.  IEEE,  Vol .  17,  No.  7,  July  1979. 

Recent  Oral  Presentations  Under  Time  Domain  Studies 

D.L.  Moffatt  and  T.C.  Lee,  "Time-Dependent  Radar  Target  Signatures 
Synthesis  and  Detection  of  Electromagntics  Authenticity  Features", 
presented  at  NATO  Advanced  Research  Workshop  on  Inverse  Methods  in 
Electromagnetic  Imaging,  Bad  Windshelm,  Federal  Republic  of  Germany, 
September  1983.  Invited. 

D.L.  Moffatt,  "Contributions  by  Professor  E.M.  Kennaugh  to  High 
Resolution  Broadband  Radar  Pol arlmetry",  presented  at  the  Second 
Workshop  on  Polarlmetric  Radar  Technology,  Huntsville,  Alabama,  May 
1983.  Invited. 

D.L.  Moffatt,  N.  Wang,  C.Y.  Lai  and  T.C.  Lee,  "Special  Response 
Waveforms  -  Cavity  and  Related  Structures",  presented  at  "Workshop  on 
Research  Techniques  in  Wave  Propagation  and  Scattering",  The  Ohio  State 
University,  October  18-0ctober  21,  1982.  Invited. 

D.L.  Moffatt,  "Time  Domain  Electromagnetics  -  An  Elementary  Review", 
presented  to  the  IEEE  Local  Chapter  on  Antennas  and  Propagation  at 
Kharagpur,  India  on  March  24,  1983.  Invited. 


E.M.  Kennaugh  and  D.L.  Moffatt,  "Improved  Response  Waveforms  for  the 
Circular  Disk",  presented  at  the  1982  Joint  International  IEEE/ 

APS  Symposium  and  National  Radio  Science  Meeting,  Albuquerque,  New 
Mexico,  May  1982. 

D.L.  Moffatt  and  Chun-Yue  Lai,  "Electromagnetic  Scattering  by  Loaded 
Open  Circular  Waveguides",  National  Radio  Science  Meeting,  University  of 
Colorado,  January  1982. 

D. L.  Moffatt  and  J.G.  Bagby,  "Approximate  Characteristic  Equations  for 
Simple  Structures",  National  Radio  Science  Meeting,  University  of 
Colorado,  January  1982. 

E. M.  Kennaugh,  "Canonical  Electromagnetic  Response  Waveforms  and  Their 
Estimation",  International  Union  of  Radio  Science,  XX*h  General 
Assembly,  Washington,  D.C.,  August  1981. 

T.W.  Johnson  and  D.L.  Moffatt,  "Electromagnetic  Scattering  by  Open 
Circular  Waveguides",  National  Radio  Science  Meeting,  Los  Angeles, 
California,  June  1981. 

E.M.  Kennaugh,  "Prediction  of  Cavity  and  Natural  Resonance  Frequencies 
by  GTD",  International  IEEE  AP-S/USNC  Symposium,  Quebec,  Canada,  June 
1980. 

E.M.  Kennaugh,  "Resonances  of  Thin-Walled  Penetrable  Scatterers", 
International  IEEE  AP-S/USNC  Symposium,  Quebec,  Canada,  June  1980. 
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O.L.  Moffatt  and  K.A.  Shubert,  "Pulse  Response  Waveforms  of  Aircraft", 
National  Radio  Science  Meeting,  Seattle,  Washington,  June  1979. 


D.L.  Moffatt  and  C.M.  Rhoads,  "An  Update  on  Naval  Vessel 
Identification",  National  Radio  Science  Meeting,  Seattle,  Washington, 
June  1979. 

D.L.  Moffatt,  "A  Chronological  History  of  Radar  Target  Imagery  at  The 
Ohio  State  University",  National  Radio  Science  Meeting,  Seattle, 
Washington,  June  1979. 

Recent  Technical  Reports  Under  Time  Domain  Studies 

Joint  Services  Electronics  Program  Annual  Report  710816-1,  December 

1978. 

"The  Calculation  of  Far  Field  Scattering  by  a  Circular  Metallic  Disk", 
Report  710816-2,  D.B.  Hodge,  February  1979. 

"Electromagnetic  Scattering  by  a  Metallic  Disk",  Report  710816-3, 

D.P.  Mithouard  and  D.B.  Hodge,  September  1979.  Thesis. 

Joint  Services  Electronics  Program  Annual  Report  710816-6,  December 

1979. 

"Electromagnetic  Scattering  by  Open  Circular  Waveguides",  Report 
710816-9,  T.W.  Johnson  and  D.L.  Moffatt,  December  1980.  Dissertation. 
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Joint  Services  Electronics  Program  Third  Annual  Report  710816-10, 
December  1980. 

Joint  Services  Electronics  Program  Fourth  Annaul  Report  710816-11, 
December  1981. 

"Transient  Current  Density  Waveforms  on  a  Spherical  Scatterer",  Report 
710816-13,  E.M.  Kennaugh  and  D.L.  Moffatt,  October  1983. 

Technical  Reports  in  Preparation  Linder  Time  Domain  Studies 

E.M.  Kennaugh,  D.L.  Moffatt  and  N.  Wang,  "K-pulse  for  Transmission 
Unes". 

D.L.  Moffatt,  J.G.  Bagby,  and  H.T.  Kim,  “A  Geometrical  Approach  to 
Characteristic  Equations". 

T.C.  Lee  and  D.L.  Moffatt,  "Approximate  Methods  for  Obtaining  the 
Complex  Natural  Electromagnetic  Oscillations  of  an  Object". 

Ph.D.  Dissertations  Under  time  Domain  Studies 

T.C.  Lee,  1983 

T.W.  Johnson,  1980 

L.C.  Chan,  1978 

M.Sc.  Thesis  Under  Time  Domain  Studies 
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F.  Transient  Signature  Measurements  of  Radar  Targets 
for  Inverse  Scattering  Research 

Researchers:  J.D.  Young,  Research  Scientist  (Phone:  (614) 

422-6657) 

E.K.  Walton,  Senior  Research  Associate 
W.  Leeper,  Graduate  Research  Associate 

Introduction 

The  objective  of  the  effort  is  to  assemble  "complete"  complex 
scattering  cross  section  data  on  several  canonical  targets.  The  data 
are  complete  in  the  following  respects: 

1.  Both  amplitude  and  phase  of  the  backscattered  return  are 
measured. 

2.  The  data  cover  a  96:1  bandwidth. 

3.  The  results  are  actual  backscatter,  i.e.,  precisely  at  0° 
bistatic  angle. 

4.  The  backscatter  is  measured  for  three  polarizations: 
vertical -vertical ,  horizontal -horizontal ,  vertical-h(  -’'rontal. 

5.  The  measurements  were  made  at  aspects  of  every  15°  from 
0°  -  180°. 

6.  The  vertical -horizontal  (cross)  polarized  measurements  were 
made  at  aspects  out  of  the  principal  planes. 


The  complete  target  signature  data  for  canonical  targets  is  much 
more  valuable  when  they  are  interpreted  and  scattering  mechanisms  are 
identified.  The  impulse  response  signature  is  being  used  as  the  vehicle 
for  this  study. 

This  section  describes  progress  in  several  aspects  of  this  study. 
The  experimental  system  has  been  significantly  improved,  and  its 
configuration  is  reviewed.  New  data  have  been  collected,  particularly 
on  the  ellipsoid  and  disk  targets,  and  a  few  illustrative  signatures  are 
presented.  The  study  of  interpretation  of  scattering  mechanisms  is 
represented  by  a  brief  discussion  of  the  ellipsoid  features.  Finally,  a 
study  of  optimum  approaches  to  sampling  the  signature  data  in  frequency 
and  aspect  angle  is  summarized. 

1.  Measurement  System 

The  system  which  has  been  used  to  acquire  the  data  is  shown  in 
Figure  F-l.  This  differs  from  the  block  diagram  of  prior  measurement 
systems  primarily  in  the  use  of  a  Scientific  Atlanta  compact  antenna 
range  reflector  to  obtain  far-zone  radar  data  for  targets  up  to  four 
feet  in  extent.  To  the  manufacturers  knowledge,  this  is  the  first 
system  in  the  world  to  use  their  reflector  for  acquiring  transient 
response  signature  data. 

The  system  has  been  described  in  two  oral  papers  [1,2]  and  a 
written  paper  which  has  been  submitted  to  IEEE/APS  [3],  The  performance 
of  this  facility  is  an  advance  In  the  state  of  the  art,  and  it  is  being 
directly  copied  by  several  other  organizations.  A  portion  of  the 
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gure  F-l.  Diagram  of  proposed  advance  target  measurement. 


signature  effort  was  devoted  to  improvement  of  this  range,  with  the 
result  that  at  the  close  of  this  effort  the  system  had  an  effective 
noise  floor  of  approximately  -50  dBsm  from  1  to  12  GHz,  and  -  70  dBsm  in 
the  range  of  8  to  12  GHz.  Other  JSEP  research  areas  have  contributed 
concepts  which  are  being  applied  to  enlarge  the  "sweet  spot"  area  of 
this  range  and  extend  its  usable  frequency  range  to  a  final  goal  of  1  to 
100  GHz. 

For  the  new  system,  a  comparison  of  measured  and  theoretical  sphere 
response  signatures  covering  a  48:1  frequency  band  are  shown  in  Figures 
F-2  and  F-3. 

2.  Target  Signature  Measurements 

The  two  targets  which  received  the  most  attention  during  the  year 
were  the  4:2:1  ellipsoid  and  the  circular  disk.  The  coverage  of  the 
frequency,  polarization,  and  aspect  angle  coverage  of  the  ellipsoid  is 
particularly  complete.  It  Is  summarized  in  Table  F-l.  The  data  base 
on  the  ellipsoid  has  been  supplied  to  two  other  organizations  at  the 
present  time,  and  is  available  to  others.  A  report  specifying  the  data 
base  format  for  the  data  which  are  supplied  on  a  floppy  disk  has  been 
written  [4], 

A  reference  frame  for  discussion  of  aspect  angle  data  on  the 
ellipsoid  is  shown  In  Figure  F-4.  Figures  F-5a  through  F-5d  show 
approximate  impulse  response  signatures  for  the  ellipsoid  at  several 
look  angles.  Figure  F-6a  through  F-6e  shows  cross-polarized  impulse 


TABLE  F-l 


DATA  BASE  SPECIFICATIONS  FOR  ELLIPSOID  SCATTERING 
SIGNATURE  MEASUREMENTS 

(a) 

1.5  ELLIPSOID 


TILT 

POL. 

FLAT 

EDGE -UP 

TIP-UP 

0° 

V 

2-12* 

2-12* 

2-12* 

0° 

H 

2-12* 

2-12* 

2-12* 

0° 

X 

1-12 

30° 

X 

1-12* 

1-12* 

45° 

V 

1-4.5 

1-6.3 

4-12 

6-12 

45° 

H 

1-4.5 

1-4.5 

4-12 

4-12 

45° 

X 

1-12* 

1-12* 

60° 

X 

1-12* 

1-12* 

(b) 

12.0  INCH  ELLIPSOID 


0° 

H 

6-6.3 

6-12 

30° 

V 

6-12 

45° 

V 

1.0-6. 3 
6-12 

1.3-6 

6-12 

0 

in 

H 

1-6.3 

6-12 

1-6.3 

6-12 

45° 

X 

1-6.3 

1-12* 

1-6.3 

6-12 

All  measurements  at  15°  aspect  intervals. 

TILT  is  the  elevation  of  the  major  axis  from  horizontal 
FLAT  means  with  the  2  axis  horizontal. 

EDGE -UP  means  with  the  1  axis  horizontal. 

TIP-UP  means  with  the  4  axis  vertical. 
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Figure  F-2.  Calibrated  waveforms. 

a.  sphere 

b.  second  sphere 

c.  background. 


AMPLITUDE  (dBSM) 


Figure  F-3.  Frequency  domain  waveforms. 

a.  sphere 

b.  second  sphere 

c.  background. 
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(a)  NORMALIZED  SPHERE  IMPULSE  RESPONSE 

(b)  ICFOHPI.  4:2:1  ELLI PSOID  ,  FLAT  ,  0*  (  NOSE -ON  ) ,  HP 

(c)  TCF30H2 .  4:2:1  ELLIPSOID, FLAT, 30*  OFF  NOSE-ON, HP 

(«J)  TCD901 .  4:2=1  ELLIPSOID,  F  LAT  ,  90°  (  BROADSIDE  ),  HP 


TIME  ,  NORMALIZED 


Figure  F-5.  Impulse  response  waveforms  for  4:2:1  ellipsoid,  compared  to 
the  sphere  impulse  response. 


■  ili  :ui'- 


0.0  0.0  1.0  2.0  3.0  4.0 


O. 

■n 


(a)  CDXK2.  XPOL  IMPULSE  RESPONSE, 45  TILT,  FLAT,  270° 


(b)  C0XK4.  XPOL  IMPULSE  RESPONSE  ,45  TILT,  FLAT,  300* 
(C)  C0XK5A.  XPOL  IMPULSE  RESPONSE  .45TILT  .FLAT, 315* 

( d )  CDXKS.  XPOL  IMPULSE  RESPONSE  ,  45TILT,  FLAT,  330* 

(e)  CDXLO.  XPOL  IMPULSE  RESPONSE  .45TILT,  FLAT, 360° 


TIME,  NORMALIZED 


Figure  F-6.  Cross-polarized  Impulse  response  vs  aspect  angle  for  a 
4:2:1  ellipsoid  (10:1  bandllmlted  approximation). 
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response  signatures  for  the  ellipsoid.  These  data  have  been  discussed 
in  a  separate  report  [5],  and  will  be  submitted  for  publication. 

A  pair  of  impulse  response  signatures  for  the  circular  disk  at  look 
angles  shown  on  the  figures  are  shown  in  Figures  F-7  and  F-8.  These 
data  are  being  used  by  Professor  Moffatt  and  colleagues  in  his  JSEP 
studies  in  order  to  clarify  and  assess  the  differences  in  the  predicted 
signatures  of  the  circular  disk  arrived  at  by  different  analytical 
techniques. 

3.  Interpretation  of  Scattering  Mechanisms  from 
Impulse  Signatures 

Following  up  on  last  year's  JSEP  studies  in  this  area,  reports  and 
papers  refining  and  interpreting  the  sphere-capped  cylinder  and  spheroid 
data  have  been  presented  and  submitted  for  publication  [6,7,8].  The 
results  show,  in  particular,  how  easily  UTD  mechanisms  can  be 
understood  in  the  time  domain,  and  how  the  scattering  mechanisms  are 
clearly  separable  in  a  time-history  point  of  view. 

The  same  techniques  may  be  applied  to  the  signatures  of  Figure  F-5. 
For  the  nose-on  ellipsoid  case,  the  specular  reflection  off  the  nose  and 
the  creeping  wave  around  the  rear  of  the  target  are  identifiable,  and 
their  amplitude  and  position  are  related  correctly  to  basically  similar 
returns  from  a  sphere.  As  the  ellipsoid  is  rotated  in  azimuth,  the 
specular  term  grows  as  the  local  radius  of  the  specular  point  increases, 
and  the  time  delay  between  the  specular  return  and  the  creeping  wave 
decreases  proportional  to  the  projection  of  the  target  along  the 
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line-of-sight.  At  broadside,  the  time  of  arrival  and  character  of  the 
creeping  wave  term  is  significantly  different  for  hh  and  vv 
polarizations,  as  is  expected  due  to  the  different  paths  of  the  creeping 
wave  propagation  around  the  body. 

4.  Optimum  Sampling  of  a  Scattering  Signature  in 
Aspect  Angle  and  Frequency 

Acquiring  comprehensive  target  scattering  signature  data  is 
equivalent  to  sampling  its  scattering  properties  in  aspect  angle  and 
frequency.  The  accuracy  of  the  data  is  improved  and  cost  of  data 
acquisition  reduced  if  optimal  sampling  is  done.  The  following 
paragraphs  summarize  present  progress  on  a  study  in  this  area. 

This  discussion  will  explore  applications  of  the  N-dimensional 
sampling  theorem  to  invese  scattering.  More  specifically,  the 
application  is  on  the  reconstruction  of  the  impulse  response  [9]  of  a 
target  using  one-dimensional  sampling  theorem. 

In  theory,  an  impulsive  electric  field  is  sent  to  the  target  at  a 
given  aspect  angle.  The  measurement  of  the  back  scattered  echo  is 
referred  to  as  the  impulse  response  of  the  target  at  that  particular 
aspect  angle.  Another  approach  to  obtain  the  Impulse  response  Is 
similar  to  deconvolution:  dividing  the  complex  frequency  response 
spectrum  of  the  target  by  its  input  signal  frequency  spectrum;  then 
the  result  Is  Inverse  Fourier  Transformed  to  give  the  impulse  response. 
In  mathematical  terms: 


151 


FU)  |  H(u)  |  C(oi) 

I _ I 

Let  i)  f(t)  be  the  input  in  time, 

F ( oj)  be  its  Fourier  transformed  pair, 
ii)  H(t)  be  the  impulse  response  of  the  target, 

H(m)  be  its  Fourier  transformed  pair, 
iii)  C(t)  be  the  output  in  time  (i.e.,  the  scattering) 
C(w)  be  its  Fourier  transformed  pair, 
then  by  the  convolution  theorem  in  Fourier  transform. 


According  to  Kennaugh  and  Moffatt  [10],  the  impulse  response  will 
decay  exponentially  for  large  values  of  t-r/c,  where 


t  =  time 

r  =  distance  between  observation  point  and  the  origin 
c  =  speed  of  light. 


Therefore,  one  can  define  a  settling  time  when  the  impulse  response  has 
its  amplitude  embedded  in  noise.  For  all  practical  purposes,  the 
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settling  time  will  be  the  end  of  the  impulse  response,  though  the 
impulse  response  may  linger  on  for  infinite  time  in  reality.  Defining 
the  end  of  the  impulse  response  is  another  way  of  saying  that  the 
impulse  response  is  time  limited.  From  the  sampling  theorem  [11]: 

A  time  limited  signal  can  be  reproduced  from  its  discrete 
frequency  values,  if  it  is  sampled  over  the  complete  frequency 
domain  using  the  Nyquist  rate.  As  a  consequence,  the  impulse 
response  may  be  reproduced  by  measuring  its  frequency  response 
values  at  the  Nyquist  sampling  frequencies. 

Next,  the  use  of  an  N-dimensional  sampling  theorem  is  considered. 
Any  physical  object  is  finite  in  dimensions.  If  it  is  located  in  a 
rectangular  coordinate  system  (X;  Xi,  X2,  X3),  then  it  can  be  said  to 
have  limited  dimensions  in  the  X  coordinates.  Namely,  the  object  is 
confined  in: 

Xj  <  Xj  <  xj 

x*  <  x2  <  x* 

X!  <  X3  5  X3 

where  X*’  X2,  X®,  X^’  X2,  X^  are  some  real  constants. 

Now,  consider  the  impulse  response  of  such  an  object.  An  impulsive 

magnetic  field  TTi  incident  on  the  body  will  set  up  current  (J)  on  the 
surface.  As  a  result,  J  =  n  x  H1  will  start  generating  a  scattered 


field  in  all  directions.  After  H’  has  washed  over  the  target,  the 
current  created  and  its  scattered  field  will  start  to  exponentially 
decay.  Therefore,  at  one  aspect  angle,  the  start  of  the  target  response 
is  defined  by  the  Initial  time  (T-j ) ;  the  time  (Tf)  denotes  the  beginning 
of  the  exponential  decay. 

Let  X  be  the  length  of  the  object  at  an  aspect  angle 

At  be  the  time  difference  between  the  critical  and  final  time 
at  that  angle 

C  be  the  speed  of  light 

At  1  Tf  —  Ti . 

Then, 

2X  =  C[At] 

-C(Tf-Ti)  . 

Putting  the  concepts  together,  the  confinement  of  an  object  in 
space  at  one  angle  is  equivalent  to  time  limiting  of  its  impulse 
response. 

While  there  is  only  one  sampling  theorem,  there  are  different  forms 
of  the  theorem,  depending  on  how  the  time  limited  or  frequency  limited 
signal  is  assumed  to  repeat  itself.  In  the  3-D  sense,  this  relates  to 
whether  the  target  is  confined  by  a  unit  of  sphere,  parallel  piped, 
ellipsoid,  etc.,  and  as  assumption  on  whether  the  unit  is  repeated  in 
the  3-dimensional  space. 


Since  the  target  is  of  arbitrary  shape  already,  one  would  not  want 


to  complicate  the  problem  by  assuming  an  arbitrarily  shaped  unit  to 
confine  the  target.  The  logical  approach  would  be  to  confine  the  object 
by  either  a  sphere  or  parallel  piped.  A  decision  rule  may  be  derived 
after  the  efficiency  has  been  defined. 

The  following  formula  is  a  modified  version  of  the  original  formula 
defined  by  Peterson  and  Middleton  [12] 


Given: 


Conservative  estimate  of  the  volume  of  the  object 
=  Volume  of  the  smallest  sphere  enclosing  the  object 

=  efficiency  using  the  sphere. 


n_  = 


Conservative  estimate  of  the  volume  of  the  object 
Volume  of  the  smallest  peral Tel  piped  enclosing  the  object 


=  efficiency  using  the  parallel  piped. 

Rule:  If  ns  >  np,  then  use  spherical  confinement 

If  np  >  ns,  then  use  parallelplped  confinement. 

N-dimensional  sampling  theorem  (from  Petersen  and  Middleton  [12]): 

A  function  fCx)  whose  Fourier  transform  F(u>)  vanishes  over 
all  but  a  finite  portion  of  wave-number  space  can  be  everywhere 
reproduced  from  its  sampled  values  taken  over  a  lattice  of 
points  {l\Y\  +  12^2  +  •••  +  *NVn}’»  *1**2* .-•*■»!  =  0,  ±1,  ±2,... 
provided  that  the  vectors  { Vj }  are  small  enough  to  ensure 
non-overlapping  of  the  spectrum  F(to)  with  its  images  on  a 
periodic  lattice  defined  by  the  vectors  {lT^  w1*^h  Vj  .  ^ 

2irojk. 
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One  can  measure  and  reproduce  the  original  8-space  signal.  The 
procedure  for  reconstruction  is 

1.  Sample  the  frequency  response  at  the  sampling  lattice 

2.  Interpolate  the  sampled  points 

3.  Fourier  transform  the  frequency  response  to  obtain  the 
original  signal. 

Consider  the  non-overlapping  condition  in  the  N-dimensional 
sampling  theorem.  The  requirement  is  'non-overlapping'  of  the  spectrum 
F(u>)  with  its  images  on  a  periodic  lattice.  In  other  words,  the 
periodic  lattice  is  not  uniquely  defined,  e.g.,  an  arbitrary  guard  band 
can  be  put  around  the  spectrum  F(u),  thus  changing  the  periodic  lattice 
Furthermore,  the  periodicity  defined  by  {u^}  is  arbitrary.  Figure  F-9 
shows  a  2-0  spectrum  with  different  periodicity.  However,  an  efficient 
sampling  lattice  may  be  defined.  From  Petersen  and  Middleton  [12]:  "Ai 
efficient  sampling  lattice  is  one  which  uses  a  minimun  number  of 
sampling  points  to  achieve  exact  reproduction  of  a  wave-number  limited 
function".  In  other  words,  the  closest  pack  of  the  spectrum  and  its 
images  without  overlapping  will  enable  efficient  sampling. 

With  the  ease  of  calculation  in  mind,  the  following  periodicity  of 
the  spectrum  and  its  images  are  defined  as  in  Figure  F-10.  From 
Petersen  and  Middleton: 


F(  o>)  =  l  F(ul>])G(u)-uU]) 


-7-  '  7  -  -~h 

(o)  ONE  TYPE  OF 
PERIODICITY 


GUARD 

BAND 


(b)  SAME  AS  (a) 

WITH  GUARD  BAND 


- ,  -r-  -  r 

i 


(C)  PERIODICITY 
2-D  spectra. 

(a)  One  type  of  periodicity. 

(b)  Same  as  (a)  with  guard  band. 

(c)  Periodicity  different  from  (a) 
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[A]  =  *i,  12,  *3  where  *1,  *2.  *3  =  0.  ±1,  *2... 
utfc]  =  ifi*!  +  U2*2  +  U3*3  . 


a 


■cvcv: 


G{i:;)  =  TRl^p^f) 


2a)icos{Ra)i//T)cos(Riu>2)  ~  2uicos(2Rui//7) 
-2/T  [s1n{Ri<jq//7)]  [s1n(Ru)2)]  (u>2) 


For  3-Dimensional  case: 


R  v?  =  | 


i  i  in 


R  *3  =  J-1//7 


'1//7  ”| 

I  1  - 

1//7  1  IT  u2  = 

I 

1//7  | 


"  0  ”1 
0  I  1  - 

I  *  u3 
i  in  i 
I 


n//7  -| 

i  i  i 

1-1  ht  i  k 


There  is  a  u)2  factor  missinq  in  the  third  term  of  this  expression  in 
Reference  [12]. 


%  - •  ••  ’»  *. v  ’  *.  *•  /».*.  •-  *.  A  J 


2(R//2)  [oij+u^+Wj  -  2 ( o>^+o>2+(a> +  u^^)] 


-2a)jSi n(/2Rcj^ )  -  2aj2Sln(/3Ra)2)  -  2u)2Sin(/2Ru)3) 


+ ( u>^  +aj2+tjJ3 )  si  nR/Z/^Wj+u^+Wj) 


+( u>l  -a>2+(^2 ) sf  nR//  2  (a)^  -ui2+<*>3 ) 
+( ail— a>2— ) S1  nR/^2(u)^-a)2“(i>3) 
+(wj+aj2-W2)si  nR//2(w^ +012-^3) 


Implementation 


Depending  on  the  shape  of  the  object,  the  object  may  be  confined  in 
a  sphere  or  a  parallel  piped  such  that  the  isotropic  or  parallel  piped 
form  of  the  sampling  theorem  [12]  may  be  employed  respectively  and 
effectively.  But,  in  general,  most  objects  can  easily  be  fitted  into  a 
rectangular  box;  therefore,  there  is  no  advantage  to  go  to  the  more 
general  case  of  parallel  piped  confinement.  Nevertheless,  programs  are 
developed  for  both  the  isotropic  and  parallel  piped  sampling. 

Thus  far,  the  development  in  the  area  has  been  in  computer 
programming.  After  specifying  isotropic  or  parallel  piped  sampling,  a 
file  containing  sampling  angles  and  radii  in  the  2-D  frequency  domain 
will  be  output.  Binary  search  is  employed  to  arrange  the  angles  in 
ascending  order.  The  output  file  is  then  used  to  make  actual 


measurement  or  storage  selection. 


Since  the  Isotropic  or  parallel  piped  sampled  data  are  governed  by 
the  sampling  theorem,  they  are  sufficient  to  define  the  object  in  the 
2-D  frequency  domain  plane.  Furthermore,  the  sampling  theorem  may  be 
viewed  as  an  interpolation  formula.  Interpolation  to  get  the  1  aspect 
angle  swept  over  a  frequency  band  data  (which  after  Fourier  transform 
can  provide  the  impulse  response  of  that  aspect  angle)  and  interpolation 
to  get  the  RCS  measurement  for  a  particular  frequency  inside  the  swept 
frequency  band  are  available.  The  interpolation  into  2-D  rectangular 
coordinates  in  the  frequency  domain  is  possible,  but  as  pointed  out  in 
the  Mensa,  et  al .  paper  [14],  It  is  very  time  consuminq.  The  natural 
approach  to  obtain  2-D  impulse  response  of  the  object  will  be  to  confine 
the  object  in  a  square  box  and  take  samples  at  the  sampling  points,  thus 
filling  the  2-D  frequency  plane.  The  rectangular  grid  frequency  domain 
data  may  then  be  transformed  into  the  time  domain. 

Since  the  FFT  program  was  written  to  operate  on  rectangular 
coordinates  and  not  on  polar  coordinates,  the  isotropic  sampled  points 
required  a  different  approach  to  obtain  the  2-D  time  response.  A  2-D 
Fourier  transform  operating  on  polar  coordinated  points  is  required.  In 
the  Mensa,  et  al  paper  [14],  they  reduced  the  double  integral  Into  a 
single  integral  involving  convolution  using  DOlar  coordinates.  Using  a 
packaged  subroutine  FFT  three  times,  the  2-D  time  domain  waveform  was 
constructed  for  each  frequency  ring.  A  similar  approach  was  used 
successfully  here. 


5.  Frequency  Domain  Decomposition  of  Scattering  Mechanisms 
Via  the  Time  Domain 

It  has  long  been  recognized  that  the  high  frequency  scattering  from 
a  body  can  be  considered  as  a  superposition  of  local  scattering 
mechanisms.  The  mechanisms  can  be  categorized  and  the  common  ones  are 
reflection  and  edge,  corner,  tip  and  surface  diffractions.  Most  of 
these  mechanisms  can  be  calculated  accurately  in  the  high  frequency 
regime  based  upon  asymptotic  solutions.  However,  as  the  wavelength  of 
the  incident  field  becomes  comparable  to  and  smaller  than  the  size  of 
the  body,  these  asymptotic  solutions  gracefully  fail.  A  new  approach  to 
obtain  the  numerical  bi static  values  of  these  mechanisms  for 
frequencies  in  the  resonant  region  of  the  scatterer  can  be  obtained 
through  a  time  domain  approach. 

The  time  domain  approach  is  simply  to  gate  out,  from  the  impulse 
representation  of  the  scatterer,  as  much  of  the  desired  mechanism  as 
possible  without  including  any  significant  amount  of  any  other 
mechanism.  The  gate  is  nothing  more  than  a  window  concept  that  is 
common  in  communication  applications  where  only  a  finite  record  of  data 
is  available  for  processing.  The  window  symmetrically  weights  the 
desired  portion  of  the  time  response  by  a  factor  that  is  unity  at  its 
center  and  smoothly  goes  to  zero  at  its  endpoints.  The  window  portion 
of  the  time  response  is  then  numerically  Fast  Fourier  Transformed  (FFT) 
to  obtain  the  frequency  domain  representation  of  the  scattering 
mechanism.  The  original  time  domain  representation  of  the  scatterer  may 
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be  obtained  directly  through  a  time  domain  solution,  time  domain 
measurements,  an  eigenfunction  solution  or  wideband  frequency 
measurements.  In  any  event,  the  original  time  domain  scattering 
response  is  an  approximation  to  the  true  impulse  response. 

This  technique  has  been  successfully  employed  on  canonical  shapes 
where  the  eigenfunction  solution  is  readily  available,  i.e.,  the  sphere 
and  the  infinitely  long  circular  and  elliptic  cylinders.  Here,  the  band 
limited  impulse  response  was  obtained  by  applying  the  FFT  to  unwindowed, 
lowpass  frequency  domain  information.  The  results  are  excellent 
throughout  most  of  the  data  record  available.  The  endpoints  of  the 
processed  record  become  distorted  due  to  the  standard  limitations 
with  data  windowing.  The  extent  of  this  distortion  is  small  and 
located  in  the  Rayleigh  region  and  high  frequency  PO  region  of  the 
scatterer  where  other  techniques  are  available  to  predict  the 
scattering. 

An  example  of  the  results  for  this  technique  are  shown  in  Figures 
F-12  through  F-16.  They  pertain  to  the  far  field  monostatic  scattering 
due  to  an  incident  magnetic  plane  wave  field  polarized  parallel  to  the 
cylinder  axis  from  a  perfectly  conducting  circular  cylinder  with  a  three 
inch  radius.  Comparisons  of  the  reflected  field  obtained  by  a 
Luneberg-Kl ine  expansion,  subtracting  the  creeping  wave  field  from  the 
eigenfunction  solution  and  the  time  domain  approach  are  shown  in  Figure 
F-12.  A  comparison  of  the  creeping  wave  field  as  predicted  by  GTD  and 
the  time  domain  approach  is  shown  in  Figure  F-13.  A  comparison  of 


scattered  field  based  upon  an  eigenfunction  solution  and  that  based  upon 


the  sum  of  the  the  reflected  and  creeping  wave  fields  as  predicted  by 
the  time  domain  approach  is  shown  in  Figure  F-14.  The  previously 
mentioned  end  effects  are  illustrated  in  these  figures.  The  approximate 
time  response  for  the  reflection  and  creeping  wave  mechanisms  are  shown 
in  Figure  F-15.  These  were  obtained  by  appropriately  windowing  the 
frequency  domain  representations  for  the  mechanisms  and  applying  the  FFT 
to  them.  A  comparison  between  the  directly  transformed  eigenfunction 
solution  and  the  sum  of  the  two  mechanisms  shown  in  Figure  F-15  is  shown 
in  Figure  F-16.  The  small  disagreement  in  late  time  is  attributed  to 
the  finite  length  of  the  time  window. 

The  main  benefit  from  this  approach  is  the  determination  of 
scattering  mechanisms  in  the  resonant  region  of  a  body. 
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Figure  F-12.  Reflected  field  comparison. 
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APPENDIX  I 

PROJECT  TITLES  AND  ABSTRACTS 

Project  784652  An  Advanced  Prototype  System  for  Locating  and  Mapping 
Underground  Obstacles 

The  objective  of  this  program  is  to  develop  a  portable  video  pulse 
radar  system  for  locating  and  mapping  underground  objects  to  a  depth  of 
10-15  feet.  The  emphasis  is  on  improving  signal  processing  techniques 
and  optimizing  system  performance  to  improve  target  resolution. 

Project  710816  Block  Funded  Support  for  Electromagnetic  Research 

This  is  research  in  the  area  of  electromagnetic  radiation  and 
scattering  including:  (1)  extension  of  the  Geometric  Theory  of 
Diffraction  (GTD)  for  convex  surfaces,  edges,  vertices  and  time  domain 
solutions;  (2)  the  GTD  combined  with  the  Method  of  Moments  (MM);  (3) 
extension  of  the  MM  codes  utilizing  polyqon  surface  current  patches  and 
wire-patch  attachment  modes;  (4)  transient  electromagnetic  phenomena 
including  target  identification,  radar  imagery,  K-pulse  techniques  and 
scattering  from  a  thin,  circular  disk;  and  (5)  improving  the  free  space 
scattering  range  at  The  Ohio  State  University  ElectroScience  Laboratory 
to  obtain  accurate  polarization  matrix  information  on  selected  reference 
targets. 
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Project  710964  Analysis  of  Ariborne  Antenna  Patterns 

The  objectives  of  this  program  are  to:  (1)  improve  the  aircraft 
model  for  far  field  pattern  computations  by  considering  a  more  realistic 
vertical  stabilizer;  (2)  study  various  ways  to  model  more  general  antenna 
types  such  as  a  monopole  in  the  presence  of  directors;  (3)  examine 
various  flat  plate  simulation  codes;  and  (4)  compare  various  calculated 
results  with  measurements  supplied  by  NASA  (Langley). 

Project  711930  Radar  Cross  Section  Studies 

The  objective  is  to  establish  the  GTD  techniques  required  to  treat 
the  scattering  of  missiles  and  to  develop  experimental  techniques  to 
obtain  the  scattered  fields. 

Project  712242  Formulate  Quasi -Optical  Techniques  for  Antennas  at  UHF 

The  goal  of  this  program  is  to  increase  the  electromagnetic 
effectiveness  of  Navy  ships  by  developing  low  cost,  integrated, 
systematic  EM  design  procedures. 


Project  712257  Application  of  Optical  Computing  Techniques  to  Jet 
Engine  Control 

This  program  involves  the  following  tasks:  (1)  survey  and  document 
control  requirements  for  jet  engines  using  information  supplied  by 
sponsor;  (2)  survey  and  document  the  field  of  optical  computing  as 
applied  to  jet  engine  controls;  (3)  generate  a  report  listing  the  various 
schemes  and  comparing  them  for  speed,  information  processing  capability, 
and  ability  to  withstand  the  necessary  environmental  conditions;  and 
(4)  make  recommendations  as  to  the  scheme  most  likely  to  satisfy  the 
requirments  stated  in  (1). 

Project  712861  Coal  Pile  Electromagnetic  Sensing  Research 

This  project  involves  a  research  program  in  electromagnetic 
subsurface  remote  sensing  as  applied  to  accurate  estimation  of  the 
quantity  of  coal  in  a  large  coal  pile.  This  problem  is  important  to 
inventory  control  at  coal  fired  generating  stations  in  the  electric  power 
industry.  Coal  quantity  depends  on  both  density  and  volume  of  the  pile, 
and  means  for  remote  sensing  of  both  of  these  parameters  is  sought. 

Project  713169  EO  Device  Signature  Reduction 

This  is  a  classified  program. 


Project  713176  Xenon  Probe  Laser  for  Atmospheric  Studies 

The  objective  is  to  construct  and  test  a  versatile  laser,  rugged 
enough  for  field  operation  and  tunable  to  various  lines  in  the  2-11.3  pm 
range,  for  use  as  an  atmospheric  probe  laser. 

Project  713206  Advanced  RCS  Reduction 
This  is  a  classified  program. 

Project  713302  Design  of  Dual  Band  Antennas 

The  Navy  frequently  has  need  for  antennas  which  will  operate  at  more 
than  one  band  of  frequencies.  This  project  addresses  the  design  of  a 
dual  band  reflector  antenna  utilizing  a  dichroic  surface  design  that  is 
based  on  extensive  experience  at  the  ElectroScience  Laboratory  with 
transparent  metallic  surfaces. 

Project  713402  Microstrip  Analysis  Techniques 

This  study  extends  the  state  of  the  art  in  microstrip  analysis 
techniques  by  considering  mutual  coupling  between  microstrip  antennas  on 
a  flat  strip  and  on  a  cylinder. 


Project  713533  Engineering  Calculations  for  Communications  Satellite 


Systems  Planning 


The  objective  of  this  program  is  to  enhance  the  computational 
capability  of  NASA/Lewis  Research  Center  for  calculating  the  performance 
of  proposed  communications  satellite  systems.  In  particular,  the  effort 
addresses  the  potential  interference  between  services  sharing  a  common 
frequency  band. 


Project  529623  Terrascan  Improvement 


The  objective  is  the  development  of  an  improved  Terrascan  unit  for 
the  location  of  buried  pipes.  Two  areas  are  being  addressed: 

1)  improvement  of  the  sensor  head  (antenna)  and  2)  improvements  in  the 
ease  of  operation  of  the  unit  and  its  adaptability  to  a  wide  range  of 
soil  conditions. 


Project  713656  Study  of  Meteorological  Factors  in  Earth  Satellite 
Propagation 


The  objectives  of  this  program  are  to  1)  conduct  gain  degradation 
measurements  on  Comstar  D/4  28.56  Ghz  synchronous  satellite  beacon  and 
2)  conduct  path  diversity  measurements  using  two  28  Ghz  radiometers 
separated  by  approximately  12  Km. 
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Project  713712  Antenna  and  Radome  Analysis  and  Design 


The  objective  of  this  study  is  to  analyze  the  performance  of  a 
design  candidate  antenna  with  a  tunnel.  Computer  models  are  to  be 
developed  for  the  scan  plane  so  that  worst  case  performance  can  be 
predicted. 

Project  312688  Geophysical  Research 

The  objective  of  this  program  is  to  determine  the  state-of-the-art 
in  underground  exploration  using  electromagnetic  techniques  and  define  a 
program  for  an  electromagnetic  exploration  system  (EMES ) . 

Project  713768  A  Theoretical  and  Experimental  Study  of  a  Potential 

Mobile  System  for  Evaluating  Electrical  Parameters  of  the 
Earth 

The  objective  of  this  study  is  to  evaluate  a  novel  approach  for 
rapidly  estimating  the  electrical  parameters  of  the  earth  over  a  broad 
frequency  band. 
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Project  713774  The  Effects  of  Atmospheric  Water  Vapor  on  Infrared 
Propagation 

The  objective  of  this  program  is  to  make  measurements  of  the 
absorption  of  infrared  radiation  by  pure  and  pressure-broadened  water 
vapor  samples  at  several  different  frequencies  using  a  frequency-doubled 
CO2  laser.  Results  of  these  measurements  are  to  be  incorporated  into  an 
analysis  of  continuum  absorption  as  wings  of  strong  absorption  lines. 

Project  714100  Navigation  System  Simulation 

The  objective  of  this  program  is  to  evaluate  the  simulation  and 
performance  of  the  GPS/JTIDS  Adaptive  Multifunction  Array  antenna  design. 

Project  714119  Support  Services  of  Electromagnetic  Design  of  RAM/RAS 
and  Non-Specular  Scattering 

This  is  a  classified  program. 

Project  714190  Resonant  Structure  NCTR 

The  objective  of  this  research  project  is  to  measure  the  complex 
polarization  scattering  matrix  of  resonant  structures. 
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Project  713774  The  Effects  of  Atmospheric  Water  Vapor  on  Infrared 
Propagation 

The  objective  of  this  program  Is  to  make  measurements  of  the 
absorption  of  Infrared  radiation  by  pure  and  pressure-broadened  water 
vapor  samples  at  several  different  frequencies  using  a  frequency-doubled 
CO2  laser.  Results  of  these  measurements  are  to  be  Incorporated  Into  an 
analysis  of  continuum  absorption  as  wings  of  strong  absorption  lines. 

Project  714100  Navigation  System  Simulation 

The  objective  of  this  program  Is  to  evaluate  the  simulation  and 
performance  of  the  GPS/JTIDS  Adaptive  Multifunction  Array  antenna  design. 

Project  714119  Support  Services  of  Electromagnetic  Design  of  RAM/RAS 
and  Non-Specular  Scattering 

This  Is  a  classified  program. 

Project  714190  Resonant  Structure  NCTR 

The  objective  of  this  research  project  Is  to  measure  the  complex 
polarization  scattering  matrix  of  resonant  structures. 
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Project  714258  Advanced  SATCOM  ECCM  Techniques 
This  Is  a  classified  project. 

Project  714349  Radar  Cross  Section  of  Aircraft/Avionics  Systems 
This  Is  a  classified  project. 

Project  714381  Electromagnetic  Reradiation  Theoretical  Prediction  for 
Complex  Building  Structures 

The  purpose  of  this  research  project  Is  to  develop  building  models 
with  sufficient  structural  detail  to  allow  prediction  of  phase  and  gain 
patterns  for  rooftop-mounted  antenna  systems  against  any  specified  HF 
SI gnal . 

Project  714414  Construction  of  Vehicle  Tracking  -  Phase  II 

This  project  provides  for  the  construction  of  a  ground-based  optical 
tracking  system.  This  system  will  use  three  Retlcon  digital  line 
scanning  cameras  to  track  a  number  of  light  sources  moving  over  a  portion 
of  the  Vehicle  Dynamics  Area  of  the  Transportation  Research  Center  of 
Ohio. 
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Project  714428  Radar  Cross-Section  (RCS)  Controls 

This  project  seeks  to  develop  design  parameters  for  slotted  metal 11 
radomes  and  RCS  control  of  OSR  antenna  assembly. 

Project  714505  Advanced  Adaptive  Antenna  Techniques 

The  objective  of  this  research  program  is  to  study  the  effect  of 
frequency  hopping  on  the  performance  of  an  adaptive  array  and  the 
performance  of  the  Frost  algorithm  with  multiple  constraints  for  use  In 
communications  systems. 

Project  312145  DBS  Antenna  Measurements 

The  purpose  of  this  program  Is  to  measure  the  performance  of 
candidate  antennas  for  DBS  application. 

Project  714614  Radar  Cross  Section  Studies 

A  bl static  near-  and  far-fleld  solution  will  be  developed  using  the 
Uniform  Geometrical  Theory  of  Dllffractlon  (UTD )  method  for  the 
scattering  from  airborne  targets.  Specifically,  basic  shapes  such  as 
finite  flat  plates,  cylinders,  ellipsoids,  and  cone  frustrums  will  be 
analyzed  separately  and  in  combination  with  one  another  to  provide  an 
efficient  approximate  model  for  scattering  bodies. 


Project  714861  On-Aircraft  Antennas 

The  objective  Is  to  develop  the  theory  and  the  practical  analysis 
techniques  necessary  to  predict  the  performance  of  line  array  of 
stacked  yagl  antennas  embedded  In  the  leading  edge  of  the  wing  of  an 
aircraft.  The  method  used  to  predict  these  antenna  patterns  shall  be  to 
adapt  the  aircraft  antenna  pattern  prediction  programs  developed  by  the 
Ohio  State  University  together  with  the  geometrical  theory  of  diffraction 
and  related  work  also  developed  at  the  Ohio  State  University. 

Project  714921  Computation  of  the  Performance  of  an  Element  Antenna 
Array  Mounted  on  an  Aircraft 

The  potential  communication  capability  provided  by  an  array  of 
antenna  elements  (5  monopoles)  mounted  on  an  aircraft  In  the  presence  of 
Interference  was  assessed.  The  study  Involved  the  computation  of  the 
radiation  pattern  of  each  of  the  monopoles  In  situ,  i.e.,  mounted  as 
specified  and  the  combined  performance  of  the  array  operating  In  an 
adaptive  mode.  The  Geometrical  Theory  of  Diffraction  was  used  to 
compute  the  radiation  patterns  of  Individual  monopoles.  The  Least  Mean 
Square  and  Power  Inversion  algorithms  were  used  to  study  the  performance 
of  the  array  In  adaptive  modes. 

Project  714946  Radar  Cross  Section  Studies 


This  Is  a  classified  project 


Project  715010  Stereologlcal  Computer-Assisted  Fluorometry:  Development 
of  Techniques  for  Analysis  of  Biological  Material 

This  project  Involved  the  design  of  a  prototype  device  to  use 
optical  fibers  to  transmit  data  on  cell  samples  Into  a  digital  computer 
for  statistical  processing. 

Project  715147  Basic  Electromagnetic  Scattering  Studies 

The  research  Is  directed  toward  the  further  development  of  the 
analytic  and  numerical  techniques  for  generating  the  scattered  fields  of 
prototype  aircraft  configurations. 

Project  715267  Mathematical  Models  to  Predict  the  Near  Field  and  Far 
Field  RF  Characteristics  of  Realistic  Jet  Intakes 

This  research  project  Is  concerned  with  an  analysis  of  the 
scattering  by  circular  and  rectangular  Inlet  shapes  In  the  near  field 
when  the  excitation  Is  a  plane  wave.  These  inlets  are  terminated  by  an 
Impedance  surface  within  at  some  distance  from  the  open  end. 

Project  715305  Design  of  Support  Structure  for  Low  RCS  Targets 

This  program  developed  and  constructed  a  low  cross-section  target 
support  pedestal  for  radar  cross-section  measurements. 


Project  71S346  Design  and  Calculate  Four  Panels 
This  Is  a  classified  project. 

Project  715354  Metallic  Radome  Design  and  Test 
This  Is  a  classified  project. 

Project  715370  Research  In  Transient  Response  Scattering  of  Antennas 

This  program  Is  Investigating  the  Impulse  time  domain  response 
scattering  of  several  antennas.  Compact  range  measurements  were  obtained 
covering  frequencies  from  1  to  18  GHz.  The  process  Is  being  developed 
for  use  as  a  standard  antenna  diagnostic  tool. 

Project  715498  Support  of  the  NASA  Propagation  Studies  and  Measurement 
Program  -  Meteorological  Factors  in  Earth-Satellite 
Propagation 

Two  experiments  are  being  performed  on  the  effects  of  precipitation 
on  satellite  communications  near  30  GHz.  The  first  evaluates  the 
Improvement  due  to  diversity  reception  on  the  ground  by  means  of  two 
radiometers  which  have  been  calibrated  against  the  Comstar  D4  satellite 
beacon  at  28.6  GHz.  The  second  Is  evaluating  the  gain  degradation  of  a 
large  paraboloid  due  to  angle-of-arrival  scintillations  of  the  beacon 
signal.  Future  experiments  are  being  planned  for  available  and  planned 
communications  satellites. 
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Project  715559  Research  on  Analysis  of  Large  Aperture  Reflector 
Antennas 

The  objective  Is  to  aid  NASA's  effort  to  advance  the  technology  of 
large  space  antenna  systems.  Available  electromagnetic  analyses  and 
computer  codes  will  be  extended  to  provide  the  capabilities  to  predict 
the  performance  of  large  deployable  reflector  antennas  with  mesh 
surfaces. 

Project  71S582  Advanced  Radar  Cross  Section  Reduction 
This  Is  a  classified  project. 

Project  715600  Laser  Radar  Cross  Section  and  Optical  Augmentation 
Descrlml nants 

The  laser  radar  cross  section  of  a  group  of  targets  was  measured  on 
the  ElectroSclence  Laboratory  optical  range  at  wavelengths  In  the  near-lR 
and  visible. 
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Project  715616  Development  of  a  Theoretical  Model  for  Ported  Coaxial 

Cable  Embedded  at  a  Uniform  Depth  In  a  Lossy  Half-Space 

The  guided  radar  uses  a  long  coaxial  cable  with  a  periodic  array  of 

slots  In  the  outer  conductor,  burled  In  the  earth  at  a  depth  of  a  few 

Inches.  When  an  intruder  walks  over  the  cable,  this  Induces  a  reflected 

wave  traveling  In  the  cable  toward  the  detector.  We  are  using  moment 
methods  and  computer  simulation  to  Investigate  the  target  coupling  and 
the  effects  of  planar  and  cylindrical  layering  In  the  earth. 

Project  715655  Advanced  Topics  In  Adaptive  Arrays 

Professor  Compton,  on  sabbatical  to  URL,  Is  exploring  mutual 
Interests  between  OSU  and  NRL  In  the  technical  area  of  Adaptive  Antenna 
Processing  Systems. 

Project  715723  Analysis  of  the  Scattering  from  an  Inlet  Geometry  with 
Lossy  Walls 

This  project  Is  concerned  with  a  study  of  the  effects  on  the 
scattering  by  an  Inlet  geometry  when  the  Inner  walls  of  the  Inlet  duct 
are  made  lossy. 
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Project  715728  Small  Array  Phase  Study 

The  purpose  of  this  study  Is  to  Investigate  the  relative  phase 
between  elements  of  a  small  array.  Theoretical  computations,  based  upon 
the  method  of  moments,  are  performed  to  determine  the  effects  of  mutual 
coupling  and  Interactions  with  the  array  support  structure  on  the 
elements  phases.  Also,  an  experimental  facility  Is  being  developed  to 
permit  the  accurate  measurement  of  the  relative  phase  of  the  array 
elements. 

Project  715731  Radar  Cross  Section  (RCS)  Prediction  Code  Modification 

The  analytical  methodology  will  be  developed  necessary  to  predict 
the  contributions  of  certain  electromagnetic  phenomena  to  the  RCS  of 
vehicles.  Specifically,  the  contribution  to  the  RCS  of  components 
partially  shadowed  by  other  components  on  the  vehicle  will  be  analyzed 
and  measurements  will  be  performed  to  validate  the  model. 


APPENDIX  III 


REPORTS  PUBLISHED  BY  ESL  OCTOBER  1982  TO  OCTOBER  1983 


529622-4  DEVELOPMENT  OF  TECHNIQUES  TO  PINPOINT  NATURAL  GAS  LEAKS  FROM 
BURIED  PIPES  WITH  AN  UNDERGROUND  RADAR  SYSTEM,  J.D.  Young, 
A.K.  Dominek,  J.  Izadlan  and  B.R.  Chen,  June  1983. 

710816-13  TRANSIENT  CURRENT  DENSITY  WAVEFORMS  ON  A  PERFECTLY  CONDUCTING 
SPHERE,  E.M.  Kennaugb  and  D.L.  Moffatt,  January  1983. 

710816-12  JOINT  SERVICES  ELECTRONICS  PROGRAM  -  FIFTH  ANNUAL, 

December  1982. 

711930-7  SCATTERING  BY  THIN  DIELECTRIC  STRIPS,  J.H.  Richmond, 

August  1983. 

710964-12  MOMENT  METHOD  ANALYSIS  OF  THE  CESSNA  150  AIRCRAFT, 

P.A.  Alexandropoulos,  E.H.  Newman  and  W.D.  Burnside,  October 
1982. 

711679-10  GENERAL  30  AIRBORNE  ANTENNA  RADIATION  PATTERN  CODE  USER'S 
MANUAL,  H.H.  Chung  and  W.D.  Burnside,  July  1982. 

711964-7  ELECTRICALLY  SMALL  ANTENNAS,  J.H.  Richmond  and  C.H.  Walter, 
December  1982. 

712242-11  QUASI -OPT I CAL  TECHNIQUES  FOR  ANTENNAS  AT  UHF  AND  ABOVE, 

R.C.  Rudduck  and  R.J.  Marhefka,  August-October  1981. 

712242-12  QUASI -OPTICAL  TECHNIQUES  FOR  ANTENNAS  AT  UHF  AND  ABOVE, 

R.C.  Rudduck  and  R.J.  Marhefka,  November  1,  1981  -  January 
31,  1982. 

712242-13  QUASI-OPTICAL  TECHNIQUES  FOR  ANTENNAS  AT  UHF  AND  ABOVE, 

R.C.  Rudduck  and  R.J.  Marhefka,  February  1,  1982  -  April  30, 
1982. 

712242-14  NUMERICAL  ELECTROMAGNETIC  CODE  -  BASIC  SCATTERING  CODE  NEC  - 
BSC  (VERSION  2)  PART  I:  USER'S  MANUAL,  R.J.  Marhefka  and 
W.D.  Burnside,  December  1982. 


712242-15  NUMERICAL  ELECTROMAGNETIC  CODE  -  BASIC  SCATTERING  CODE  NEC  - 
BSC  (VERSION  2)  PART  II:  CODE  MANUAL,  R.J.  Marhefka, 

December  1982. 

712242-16  NUMERICAL  ELECTROMAGNETIC  CODE  -  REFLECTOR  ANTENNA  CODE  NEC  - 
REF  (VERSION  2)  PART  I:  USER'S  MANUAL,  R.C.  Rudduck  and  Y.C. 
Chang,  December  1982. 

712242-17  NUMERICAL  ELECTROMAGNETIC  CODE  -  REFLECTOR  ANTENNA  CODE  NEC  - 
REF  (VERSION  2)  PART  II:  CODE  MANUAL,  Y.C.  Chang  and 
R.C.  Rudduck,  December  1982. 

712242-18  QUASI-OPTICAL  TECHNIQUES  FOR  ANTENNAS  AT  UHF  AND  ABOVE, 

R.C.  Rudduck  and  R.J.  Marhefka,  August  1,  1979,  -  December 
31,  .’982. 

712661-4  ANALYSIS  OF  ELECTROMAGNETIC  BACKSCATTERING  FROM  AN  INLET 

CAVITY  CONFIGURATION  -  FINAL  REPORT,  C.C.  Huang,  P.H.  Pathak, 
E.Y.  Lai  and  D.L.  Moffatt,  October  1982. 

712692-3  A  USERS  MANUAL  FOR  ELECTROMAGNETIC  SURFACE  PATCH  (ESP)  CODE: 
VERSION  II  -  POLYGONAL  PLATES  AND  WIRES,  E.H.  Newman  and 
P.  Alexandropoulos,  September  1983. 

712949-6  THEORY  OF  FIELD  STRUCTURE  FOR  PORTED  COAXIAL  CABLE  EMBEDDED 
IN  LOSSY  MEDIA,  J.H.  Richmond  and  R.J.  Garbacz,  March  1983. 

713206-3  USERS  MANUAL  FOR  COMPUTER  PROGRAM  'CADLAC',  E.K.  English  and 
B.A.  Munk,  January  1983. 

713206-4  SCATTERING  FROM  PERIODIC  ARRAYS  OF  ARBITRARILY  ORIENTED 
DIPOLE  ELEMENTS  IMBEDDED  IN  A  GENERAL  STRATIFIED  MEDIUM, 

E.K.  English  and  B.A.  Munk,  January  1983. 

713206-5  THE  SCATTERING  OF  PLANAR  ARRAYS  OF  ARBITRARILY  SHAPED  SLOT 
AND/OR  WIRE  ELEMENTS  IN  A  STRATIFIED  DIELECTRIC  MEDIUM, 

B.A.  Munk  and  L.W.  Henderson,  July  1983. 

713206-6  SCATTERING  FROM  N  ARRAYS  OF  FLAT,  LOSSY  ELEMENTS  IN  A 
STRATIFIED  MEDIUM  WITH  ARBITRARY  y  AND  e,  B.A.  Munk  and 
E.K.  English,  July  1983. 

713206-7  USER'S  MANUAL  FOR  COMPUTER  PROGRAM  'INDIAN',  E.K.  English  and 
B.A.  Munk,  September  1983. 

713206-8  USER'S  MANUAL  FOR  COMPUTER  PROGRAM  'FLATBK' ,  E.K.  English  and 
B.A.  Munk,  September  1983. 


713774-2  THE  EFFECTS  OF  ATMOSPHERIC  WATER  VAPOR  ABSORPTION  ON  INFRARED 
LASER  PROPAGATION  IN  THE  5  MICROMETER  BAND,  L.G.  Walter, 

May  1983. 

713402-5  MICROSTRIP  ANALYSIS  TECHNIQUES,  E.H.  Newman,  June  1982  - 
June  1983. 

713533-3  ENGINEERING  CALCULATIONS  FOR  COMMUNICATIONS  SATELLITE  SYSTEMS 
PLANNING,  C.A.  Levis  C.H.  Martin,  D.  Gonsalvez  and  C.W.  Wang, 
June  1983. 

713581-1  USER'S  MANUAL  FOR  REFLECTOR  ANTENNA  CODE  #2,  R.C.  Rudduck, 
and  R.J.  Paknys,  September  1982. 

713581-2  COMPUTER  CODE  FOR  SPHERICAL  AND  PARABOLIC  TORUS  REFLECTOR 
ANTENNAS,  R.J.  Paknys  and  R.C.  Rudduck,  September  1982. 

713603-8  THE  PERFORMANCE  OF  A  TRIPOLE  ADAPTIVE  ARRAY  AGAINST 

CROSS-POLARIZED  JAMMING,  R.T.  Compton,  Jr.,  October  1982. 

713603-9  SOPHISTICATED  JAMMERS  AND  ADAPTIVE  ARRAYS  -  FINAL  REPORT, 

R.T.  Compton,  Jr.,  October  1982. 

713656-3  RADIOMETRICALLY  INFERRED  SPACE -DIVERSITY  ATTENUATION  AT  28.6 
GHz:  CALIBRATION  ANO  INITIAL  RESULTS,  B.A.  Pigon  and  C.A. 
Levis,  February  1983. 

713671-2  THE  CALCULATED  PERFORMANCE  OF  SOME  RADAR  AND  TRANSMISSOMETER 
SYSTEMS  IN  RAIN  AND  FOG  AT  35,  94,  140  AND  240  GHz,  VOLUMES  I 
AND  II  -  COMPUTER  CODES,  E.L.  Utt  and  C.A.  Levis,  December 
1982. 

713711-1  MEASUREMENTS  OF  SHIP  BACKSCATTER  FOR  HF  SHIP  IDENTIFICATION 
STUDIES  -  FINAL  REPORT,  E.K.  Walton,  October  1982. 

713731-1  DATA  PROCESSING  OF  RADAR  MEASUREMENTS,  E.K.  Walton,  June 
1982. 

713768-1  A  THEORETICAL  AND  EXPERIMENTAL  STUDY  OF  A  POTENTIAL  MOBILE 
SYSTEM  FOR  EVALUATING  ELECTRICAL  PARAMETERS  OF  THE  EARTH, 

P.K.  Hayes  and  L.  Peters,  Jr.,  August  1982. 

714190-1  THE  OHIO  STATE  UNIVERSITY  NCTR  DATA  BASE  FILE  STRUCTURE, 

F.S.  Chen  and  E.K.  Walton,  October  1982. 


714921-1  PERFORMANCE  OF  AIRCRAFT  MOUNTED  ADAPTIVE  ARRAYS,  I.J.  Gupta, 
January  1983. 

714505-1  ADAPTIVE  ARRAY  BEHAVIOR  WITH  PERIODIC  ENVELOPE  MODULATED 

INTERFERENCE,  A.S.  Al-Ruwais  and  R.T.  Compton,  Jr.,  December 

1982. 

714505-2  ADVANCED  ADAPTIVE  ANTENNA  TECHNIQUES,  R.T.  Compton,  Jr., 
January  1983, 

714119-2  USERS  MANUAL  FOR  COMPUTER  PROGRAM  'CADFIT' ,  E.K.  English  and 
B.A.  Munk,  November  1982. 

714215-4  NEAR  FIELD  ANALYSIS  OF  AIRBORNE  ANTENNAS,  N.  Wang  and  W.D. 
Burnside,  March  1983. 

714258-1  EFFECT  OF  MUTUAL  COUPLING  ON  THE  PERFORMANCE  OF  ADAPTIVE 
ARRAYS,  I.J.  Gupta  and  A. A.  Ksienski,  October  1982. 

714258-2  ADAPTIVE  ARRAYS  FOR  MULTIPLE  SIMULTANEOUS  DESIRED  SIGNALS, 
I.J.  Gupta,  November  1982. 

714921-2  PERFORMANCE  OF  AIRCRAFT  MOUNTED  ADAPTIVE  ARRAYS,  I.J.  Gupta, 
May  1983. 

714119-1  DISCUSSION  OF  RCS  OF  MCDONNEL  DOUGLAS  VEHICLE  -  FINAL  REPORT, 
R.J.  Marhefka,  F.L.  Volakis,  L.  Peters,  Jr.,  and  W.D. 
Burnside,  August  1982. 

714190-2  RESONANT  STRUCTURE  NCTR,  D.L.  Moffatt,  J.D.  Young,  E.K. 

Walton  and  W.  Leeper,  January  1983. 

714215-1  GEODESIC  PATHS  FOR  THE  SIDE -MOUNTED  ANTENNA  ON  AN  ELLIPSOID 
MODEL,  J.G.  Kim  and  W.D.  Burnside,  October  1982. 

714215-3  RADIATION  PATTERNS  OF  AN  ANTENNA  MOUNTED  ON  THE  OFF-MID 
SECTION  OF  AN  ELLIPSOID,  J.G.  Kim  and  W.D.  Burnside,  June 

1983. 

714215-2  RADIATION  PATTERNS  OF  AN  ANTENNA  MOUNTED  ON  THE  MIDSECTION  OF 
AN  ELLIPSOID,  J.G.  Kim  and  W.D.  Burnside,  July  1983. 

714349-5  A  MODEL  STUDY  OF  THE  RCS  OF  AN  INLET,  R.J.  Marhefka  and  L. 
Peters,  Jr.,  August  1982. 
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and  W.O.  Burnside,  September  1983. 


INTERFERENCE,  A.S.  Al-Ruwals  and  R.T.  Compton,  Jr.,  July 
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THE  PERFORMANCE  OF  AN  LMS  ADAPTIVE  ARRAY  WITH  FREQUENCY 
HOPPED  SIGNALS,  L.  Acar  and  R.T.  Compton,  Jr.,  June  1983. 

ADVANCED  ADAPTIVE  ANTENNA  TECHNIQUES,  R.T.  Compton,  Jr., 
July  1983. 

ADVANCED  ADAPTIVE  ANTENNA  TECHNIQUES,  R.T.  Compton,  Jr.,  , 
1983. 


TECHNIQUES  FOR  ANALYSIS  OF  BIOLOGICAL  MATERIAL,  S.A.  Collins, 
Jr.,  October  1983. 
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ELECTROMAGNETIC  SCATTERING  BY  OPEN  CIRCULAR  WAVEGUIDES,  D.L.  Moffatt  and 
T.W.  Johnson,  reprinted  In  Radio  Science,  Vol.  17,  No.  6,  pp.  1547-1556, 
November-Dee ember,  1982. 

REGGE  POLES,  NATURAL  FREQUENCIES,  AND  SURFACE  WAVE  RESONANCES  OF  A 
CIRCULAR  CYLINDER  WITH  A  CONSTANT  SURFACE  IMPEDANCE,  N.  Wang,  reprinted 
In  IEEE  Transactions  on  Antennas  and  Propagation,  Vol.  AP-30,  No.  6, 
November  1982. 

AN  APERTURE-MATCHED  HORN  DESIGN,  W.D.  Burnside,  reprinted  IEEE 
Transactions  on  Antennas  and  Propagation,  Vol.  AP-30,  No.  4,  July  1982. 

A  SINGLE-PROBE  ON-SITE  METHOD  OF  MEASURING  THE  DIELECTRIC  CONSTANT  AND 
CONDUCTIVITY  OF  SOFT  EARTH  MEDIA  OVER  A  1-GHz  BANDWIDTH,  P.K.  Hayes, 
reprinted  In  IEEE  Transactions  on  Geoscience  and  Remote  Sensing,  Vol. 
GE-20,  No.  4,  October  1982. 

A  HYBRID  MM-GEOMETRICAL  OPTICS  TECHNIQUE  FOR  THE  TREATMENT  OF  WIRE 
ANTENNAS  MOUNTED  ON  A  CURVED  SURFACE,  L.W.  Henderson  and  G.A.  Thiele, 
reprinted  in  IEEE  Transactions  on  Antennas  and  Propagation,  Vol.  AP-30, 
No.  6,  November  1982. 

A  WIDE-BAND  ELECTRICALLY  SMALL  SUPERDIRECTIVE  ARRAY,  E.H.  Newman  and 
M.R.  Schrote,  reprinted  In  IEEE  Transactions  on  Antennas  and 
Propagation,  Vol.  AP-30,  No.  6,  November  1982. 

THE  EFFECT  OF  RANDOM  STEERING  VECTOR  ERRORS  IN  THE  APPLEBAUM  ADAPTABLE 
ARRAY,  R.T.  Compton,  Jr.,  reprinted  in  IEEE  Transactions  on  Aerospace 
and  Electronic  Systems,  Vol.  AES-18,  No.  5,  September  1982. 

HIGH  FREQUENCY  SCATTERING  BY  A  THIN  LOSSLESS  DIELECTRIC  SLAB, 

W.D.  Burnside  and  K.W.  Burgener,  reprinted  In  IEEE  Transactions  on 
Antennas  and  Propagation,  Vol.  AP-31,  No.  1,  January  1983. 

IMPROVED  IDENTIFICATION  OF  UNDERGROUND  TARGETS  USING  VIDEO  PULSE  RADARS 
BY  ELIMINATION  OF  UNDESIRED  NATURAL  RESONANCES,  J.L.  Volakls  and  L. 
Peters,  Jr.,  reprinted  In  IEEE  Transactions  on  Antennas  and  Propagation, 
Vol.  AP-31,  No.  2,  March  1983. 


ON  THE  CURRENT  DISTRIBUTION  FOR  OPEN  SURFACES,  E.H.  Newman  and  M.R. 
Schrote,  reprinted  In  IEEE  Transactions  on  Antennas  and  Propagation, 

Vol .  AP-31,  No.  3,  May  1983. 

PREDICTION  OF  ADAPTIVE  ARRAY  PERFORMANCE,  I.J.  Gupta  and  A. A.  Kslenskl, 
reprinted  In  IEEE  Transactions  on  Aerospace  and  Electronic  Systems,  Vol. 
AE-19,  No.  3,  May  1983. 

EFFECT  OF  MUTUAL  COUPLING  ON  THE  PERFORMANCE  OF  ADAPTIVE  ARRAYS,  I.J. 
Gupta  and  A. A.  Kslenskl,  reprinted  In  IEEE  Transactions  on  Antennas  and 
Propagation,  Vol.  AP-31,  No.  5,  September  1983. 

EIGENFUNCTIONS  OF  COMPOSITE  HERMITIAN  OPERATORS  WITH  APPLICATION  TO 
DISCRETE  AND  CONTINUOUS  RADIATING  SYSTEMS,  N.  Inagakl  and  R.J.  Garbacz, 
reprinted  In  IEEE  Transactions  on  Antennas  and  Propagation,  Vol.  AP-30, 
No.  4,  pp.  571-575,  July  1982. 

FIRST-ORDER  EQUIVALENT  CURRENT  AND  CORNER  DIFFRACTION  SCATTERING  FROM 
FLAT  PLATE  STRUCTURES,  F.A,  Slkta,  W.O.  Burnside,  T.T.  Chu  and  L. 

Peters,  Jr.,  reprinted  In  IEEE  Transactions  on  Antennas  and  Propagation, 
Vol.  AP-31,  No.  4,  pp.  584-589,  July  1983. 

THE  PERFORMANCE  OF  A  TRIPOLE  ADAPTIVE  ARRAY  AGAINST  CROSS-POLARIZED 
JAMMING,  R.T.  Compton,  Jr.  reprinted  In  IEEE  Transactions  on  Antennas 
and  Propagation,  Vol.  AP-31,  No.  4,  pp.  682-685,  July,  1983. 

SCATTERING  FROM  CYLINDRICAL  INHOMOGENEITIES  IN  A  LOSSY  MEDIUM, 

L.  Peters,  Jr.  and  J.H.  Richmond,  reprinted  In  Radio  Science,  Vol.  17, 
No.  5,  pp.  975-987,  September-October  1982. 

HIGH  FREQUENCY  SCATTERING  FROM  THE  EOTES  OF  IMPEDANCE  DISCONTINUITIES  ON 
A  FLAT  PLANE,  R.  Tlberlo  and  G.  Pelosl,  IEEE  Transactions  on  Antennas 
and  Propagation,  Vol.  AP-31,  No.  4,  pp.  590-596,  July  1983. 

AN  EXTENSION  OF  THE  UNIFORM  GTD  FOR  THE  DIFFRACTOIN  BY  A  WEDGE 
ILLULMINATED  BY  A  DIPOLE  CLOSER  TO  IT'S  EDGE,  R.G.  Kouyoumjlan,  G. 

Pelosl  and  R.  Tlberlo,  Alta  Frequenza,  Vol.  52,  No.  1,  pp.  24-33,  1983. 
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JSEP  RELATED  ORAL  PRESENTATIONS 
OCTOBER  1982  TO  OCTOBER  1983 


RESONANCE  REGION  BACKSCATTERED  SIGNATURES  FROM  A  CONDUCTING  ELLIPSOID, 
J.D.  Young  and  E.K.  Walton,  1983  International  IEEE/AP-S  Symposium  and 
National  Radio  Science  Meeting,  Houston,  Texas,  May  1983. 

IDENTIFICATION  OF  SCATTERING  MECHANISMS  FROM  MEASURED  IMPULSE  RESPONSE 
SIGNATURES  OS  SIMPLE  CONDUCTING  OBJECTS,  W.J.  Leeper  and  J.D.  Young, 

1983  International  IEEE/AP-S  Symposium  and  National  Radio  Science 
Meeting,  Houston,  Texas,  May  1983. 

MOMENT  METHOD  MODELING  OF  REALISTIC  STRUCTURES,  E.H.  Newman,  52nd 
Interservice  Antenna  Group  (I SAG)  Meeting,  Annapolis,  Maryland,  April 
1983. 

A  UNIFORM  GTD  ANALYSIS  OF  THE  EM  DIFFRACTION  BY  A  THIN-DIELECTRIC  HALF 
PLANE,  P.H.  Pathak  and  R.  Rojas-Teran,  1983  International  IEEE/AP-S 
Symposium  and  National  Radio  Science  Meeting,  Houston,  Texas,  May  1983. 

CONTRIBUTIONS  BY  EDWARD  M.  KENNAUGH  TO  HIGH  RESOLUTION  BROADBAND  RADAR 
POLARIMETRY ,  D.L.  Moffatt,  Second  Workshop  on  Polarlmetrlc  Radar 
Technology  Huntsville,  Alabama,  May  1983.  To  be  published  in  the 
Proceedings  of  the  Workshop. 

PARAXIAL  EDGE  DIFFRACTION,  R.G.  Kouyoumjlan  and  O.M.  Buyukdura,  1983 
URSI  International  Symposium  In  Santiago  de  Compostela,  Spain,  August 
1983. 

THE  MODERN  GEOMETRICAL  THEORY  OF  DIFFRACTION  AND  ITS  APPLICATION,  R.G. 
Kouyoumjlan,  Istanbul  Technical  University,  Istanbul,  Turkey,  September 
1983. 

TRANSIENT  CURRENT  DENSITY  WAVEFORMS  ON  A  PERFECTLY  CONDUCTING  SPHERE, 
E.M.  Kennaugh  and  D.L.  Moffatt,  NATO  Advanced  Research  Workshop  on 
Inverse  Methods  In  Electromagnetic  Imaging,  Bad  Windsheim,  Federal 
Republic  of  Germany,  September  1983.  To  be  published  in  the  Proceedings 
of  the  Workshop. 


TIME-DEPENDENT  RESPONSE  WAVEFORMS  AND  RECOGNITION  OF  IMAGING 
AUTHENTICITY  FEATURES,  D.L.  Moffatt  and  T.C.  Lee,  NATO  Advanced  Research 
Workshop  on  Inverse  Methods  In  Electromagnetic  Imaging,  Bad  Wlndshelm, 
Federal  Republic  of  Germany,  September  1983.  To  be  published  In  the 
Proceedings  of  the  Workshop. 

THE  OHIO  STATE  UNIVERSITY  COMPACT  RADAR  CROSS-SECTION  MEASUREMENT  RANGE, 
E.K.  Walton  and  J.D.  Young,  1983  AMTA  Annual  Meeting,  Annapolis, 
Maryland,  September  1983. 


August  1982 

UNIFORM  GEOMETRICAL  THEORY  OF  DIFFRACTION,  P.H.  Pathak,  The 
International  Society  for  Optical  Engineering,  San  Diego,  California, 
August  1982. 


JSEP  RELATED  DISSERTATIONS 
OCTOBER  1982  TO  OCTOBER  1983 

AN  OPEN  SURFACE  INTEGRAL  FORMULATION  FOR  ELECTROMAGNETIC  SCATTERING  BY  A 
MATERIAL  PLATE,  M.R.  Schrote. 

JSEP  RELATED  THESES 
OCTOBER  1982  TO  OCTOBER  1983 

ELECTROMAGNETIC  MODELING  OF  ARBITRARY  SURFACES  BY  POLYGONAL  PLATES 
AND/OR  WIRES,  P.  Alexandropoulos. 

IDENTIFICATION  OF  SCATTERING  MECHANISMS  FROM  MEASURED  IMPULSE  RESPONSE 
SIGNATURES  OF  SEVERAL  CONDUCTING  OBJECTS,  W.J.  Leeper. 
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APPENDIX  VI 


JSEP  RELATED  PAPERS  SUBMITTED  FOR  PUBLICATION 


OCTOBER  1982  TO  OCTOBER  1983 


CALCULATION  OF  THE  HIGH  FREQUENCY  DIFFRACTION  BY  TWO  NEARBY  EDGES 
ILLUMINATED  AT  GRAZING  INCIDENCE,  R.  Tlberlo,  and  R.G.  Kouyoumjlan, 
accepted  for  publication  In  the  IEEE  Transactions  on  Antennas  and 
Propagation. 

SCATTERING  BY  A  STRIP  WITH  TWO  FACE  IMPEDANCES  AT  EDGE -ON  INCIDENCE,  R. 
Tlberlo,  F.  Bessl,  G.  Manara  and  G.  Pelosl,  accepted  for  publication  In 
the  Journal  of  Radio  Science. 

ON  THE  EIGENFUNCTION  EXPANSION  OF  ELECTROMAGNETIC  DYADIC  GREEN'S 
FUNCTIONS,  P.H.  Pathak,  accepted  for  publication  In  the  IEEE 
Transactions  on  Antennas  and  Propagation. 

A  UNIFORM  GTD  APPROACH  TO  EM  SCATTERING  AND  RADIATION,  R.G.  Kouyoumjlan 
and  P.H.  Pathak,  to  appear  In  Handbook  on  Acoustic,  Electromagnetic  and 
Elastic  Wave  Scattering  -  Theory  and  Experiment,  Yol.  1,  edited  by  V.K. 
Varadan  and  V.V.  varadan,  North  Holland. 

TECHNIQUES  FOR  HIGH  FREQUENCY  PROBLEMS,  P.H.  Pathak,  to  appear  In  the 
Antenna  Handbook,  edited  by  Y.T.  Lo  and  S.W.  Lee,  ITT  Howard  W.  Sams  and 
Co.,  Inc. 

AN  OPEN  SURFACE  INTEGRAL  FORMU!  ATION  FOR  ELECTROMAGNETIC  SCATTERING  BY 
MATERIAL  PLATES,  E.H.  Newman  and  M.R.  Schrote,  accepted  for  publication 
In  the  IEEE  Transactions  on  Antennas  and  Propagation. 

POLYGONAL  PLATE  MODELING  OF  REALISTIC  STRUCTURES,  E.H.  Newman,  P. 

A1 exandropoul os  and  E.K.  Walton,  submitted  for  publication  In  the  IEEE 
Transactions  on  Antennas  and  Propagation. 

TIME  DOMAIN  ELECTROMAGNETIC  SCATTERING  BY  OPEN  ENDED  CIRCULAR  WAVEGUIDE 
AND  RELATED  STRUCTURE,  D.L.  Moffatt,  C.Y.  Lai  and  T.C.  Lee,  accepted  for 
publication  In  WAVE  MOTION. 

TRANSIENT  RESPONSE  CHARACTERISTICS  IN  IDENTIFICATION  AND  IMAGING,  D.L. 
Moffatt,  J.D.  Young  and  A. A.  Kslenskl,  to  appear  In  Microwave  Imaging 
Techniques  a  series  of  books  on  "Benchmark  Papers  in  fclectrlcal 
Engineering  and  Computer  Science",  edited  by  J.B.  Thomas,  Hutchinson 
Ross  Publishing  Company,  Stroudsburg,  Pennsylvania,  1983. 
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